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ABSTRACT 


This paper documents the development of a statistical 
analysis package for the TRS-80 microcomputer. The package 
is comprised of six interactive programs which are generally 
divided into topical areas. The major emphasis is on ex¬ 
ploratory data analysis and statistical inference, however, 
probability and inverse probability distributions are also 
included. 

The programming language is TRS-80 Level II BASIC en¬ 
hanced by the input/output commands available through the 
ESF-80 (Exatron Stringy Floppy) mass storage subsystem. 
With the modification of these few commands, the package is 
compatible with most floppy disk operating systems designed 
for the TRS-80 Model I or Model III microcomputers. 

vk. This statistical analysis capability implemented on a 
relatively inexpensive system provides a useful tool to the 
student or the trained analyst without ready access to a 
mainframe computer system. 
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I. INTRODUCTION 


Over the past decade computers have played an ever 
increasing role in the day to day operation and adminis¬ 
tration of the military services.' In the Marine Corps, for 
example, computer-based systems such as MMS (Manpower Man¬ 
agement System), JUMPS (Joint Uniform Military Pay System), 
MIMMS (Marine Corps Integrated Maintenance Management Sys¬ 
tem), and many others handle virtually all personnel 
administration and personnel pay matters, equipment main¬ 
tenance and operational status reporting, supply requisi¬ 
tioning, fiscal record-keeping, etc. The computers which 
support these systems are large mainframe devices usually 
shared by two or more major commands. While the small unit 
commander certainly benefits from the overall effects of 
these large computer systems, he has little or no access to 
a computer system for projects tailored to his individual 
unit. 

The small, relatively inexpensive microcomputer sys¬ 
tems which abound today offer the small unit commander the 
possibilities of a customized computer based system for 
unit level information processing, record keeping and 
problem solving. The hardware is available and within the 
budget of these lower levels of command but will be of 
little use unless appropriate software is also available. 
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f 

The purpose of this thesis is to provide a software package 
consisting of a number of statistical analysis procedures 
which can be effectively implemented on a microcomputer, 
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II. OVERVIEW 


A. PACKAGE CONTENT 

McNeil [Ref. 1] states that an effective data analyst 
requires as his basic tools a computer and a set of routine 
procedures. By his definition, the computer can range from 
the statistican with pencil and paper to a wide variety of 
electronic computational devices. Whatever form it takes, 
the function of the computer is to faithfully execute the 
collection of routine procedures selected by the analyst. 
By implication, the key to an effective statistical analy¬ 
sis package is the proper choice of procedures to be 
included. 

The collection of "routine procedures” contained in 
this package are intended to provide the potential user 
with a basic set of widely utilized statistical analysis 
tools. These procedures have broad application and are not 
designed for any one specific scenario. The package was, 
however, designed with a specific user in mind, i.e. the 
military officer with a working knowledge of statistical 
analysis techniques and a desire to apply these techniques 
at the small unit level. The potential applications remain 
broad but typically include such items as forecasting next 
quarter’s training support requirements, temporary addi¬ 
tional duty budget, spare parts demand, etc., without 


3 



complete knowledge of the number of personnel who will be 
assigned, what changes might be made in exercise schedules, 
equipment usage rate and many other pertinent facts. In 
other words, he must make plans, decisions, and forecasts 
based upon his best estimation of the values of random 
variables. A common sense approach to this problem is to 
collect data on previous observations of the random var¬ 
iable in question and use this data as a indicator of what 
can be expected in the future. 

In order to make knowledgeable decisions based on the 
available data it is first necessary to gain insight into 
the nature and structure of the data. Procedures utilized 
for this purpose are generally classified as exploratory 
data analysis techniques. A number of the procedures 
included in the package fall under this general classifi¬ 
cation. These include subroutines which calculate sample 
statistics such as sample mean, sample variance, skewness, 
kurtosis and the seven number summary. Also included are 
data plots such as histogram, stem and leaf, normal prob¬ 
ability plot and exponential probability plot. These pro¬ 
cedures were selected because the manner in which they 
represent the data is straightforward and easily inter¬ 
preted by the trained analyst. 

After the analyst has studied the nature and structure 
of the data he is usually interested in knowing what 
inferences concerning various population parameters can be 
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drawn from the data. Statistical analysis procedures which 
are concerned with this question fall under the heading or 
statistical inference. Statistical inference can be fur¬ 
ther subdivided. If, through exploratory data analysis of 
by other means, the analyst feels that it is reasonable to 
assume knowledge of the distributional family from which 
the data sample is drawn, then, according to DeGroot 
[Ref. 2], he is dealing in the area of parametric statis¬ 
tical inference. If, on the other hand, inferences are made 
about population parameters without assuming an underlying 
distributional form, the term nonparametric statistical 
inference is applicable. Procedures which rest on implicit 
population assumptions but which are tolerant of departures 
from these assumptions are called robust procedures. 

Parametric statistical inference procedures contained 
in this package include tests of hypotheses and confidence 
intervals for population parameters under the normal, 
exponential and Bernoulli distributional assumptions. While 
these three distributions are not the only ones of potential 
interest, they are probably the most common distributions 
assumed. 

Procedures which provide robust estimates for popula¬ 
tion parameters are also included in the package. The 
robust-t (Winsorized-t and trimmed-t) procedures are 
utilized for establishing a point estimate and confidence 
interval for the population mean. The Jackknife technique 
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is used for a point estimate and interval estimate for the 
population variance. The one true nonparametric statistical 
inference procedure included in the package is a subroutine 
which computes a confidence interval for the population 
median. 

As implemented in this package, the procedures thus 
far mentioned deal with univariate data, i.e. data which 
represents a single random variable. Also included is a 
multiple linear regression procedure which, using the 
method of ordinary least squares, tests for a statistically 
significant linear relationship between a dependent 
(response) variable and one or more independent (carrier) 
variables. The regression equation coefficients which 
estimate the linear relationship and a number of statistics 
which measure the ''goodness" of the fit are also calcula¬ 
ted. Linear regression is one of the most widely utilized 
statistical analysis procedures and is a powerful tool when 
properly used. A typical application for a military unit 
might be to study the relationship between the number of 
Not Operationally Ready (NOR) pieces of equipment reported 
for each month over some period of time (the response var¬ 
iable) and such carriers as the average number of techni¬ 
cians assigned for the month, number of days spent in the 
field during the month, average response time for supply 
demands, etc. 
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ANalysis Of VAriance (ANOVA) procedures, which are 
used to test for statistically significant differences 
between two or more population means, are also included. 
This program allows for one-way classification of popula¬ 
tions, e.g. scores on the GMS test (General Military Sub¬ 
jects) made by individuals classified according to MOS 
(Military Ocupational Speciality); two-way classification, 
e.g. GMS scores made by individuals classified according to 
MOS and sex; or three-way classification, e.g. GMS scores 
made by individuals classified according to MOS, sex and 
highschool completion status. As in the case of multiple 
linear regression, the ANOVA procedure is a very powerful 
and widely utilized analysis tool. 

Finally, the package contains several of the more 
common discrete and continuous probability distributions and 
inverse distributions and three common data transformations. 
The probability distributions and inverse probability 
distributions are of interest to the analyst in their own 
right and many also serve as subroutines for other statis¬ 
tical analysis procedures contained in the package. The 
data transformations (square root, cube root and natural 
logarithm) can be useful during the exploratory data 
analysis phase in helping to uncover structure in the data. 

Except in the more trivial cases, e.g. sample mean, 
the procedures mentioned above and the algorithms utilized 
to implement the procedures will be discussed in more 
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detail in later sections of this paper. A summary of the 
procedures contained in this package is given in Figure 1. 


B. TRS-80 MICROCOMPUTER SYSTEM 

As previously stated, the analyst can accomplish the 
calculations required to perform the statistical analysis 
procedures outlined above in a variety of ways. These 
range from pencil and paper to large mainframe computer 
systems. The computations required for very small data 
sets can easily be handled with pencil and paper and per¬ 
haps a handheld calculator while the practical considera¬ 
tions of time alone dictates that very large data sets 
require a mainframe computer system for effective analysis. 
It is analysis of the moderate size data set (perhaps 
twenty to several hundred data points) which would become 
quite tedious to accomplish by hand but which might not 
merit the expense of a large computer system, assuming one 
were available, that the microcomputer handles nicely. 

A number of microcomputers are currently available and 
new makes and models are appearing rapidly. This paper 
makes no attempt to draw qualitative comparisons between 
the competing brands. No such study was undertaken. The 
TRS-80 microcomputer was selected for this project simply 
because it appeared to offer capabilities comparable to the 
other popular makes at a slightly lower price. In addition, 
at least at the time the computer was purchased, the TRS-80 
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appeared to have the widest market acceptance. A more 
detailed description of the system follows. 

This package was developed on a TRS-80, Model I, Level 
II microcomputer. The basic system consists of a 16K RAM 
(Random Access Memory) CPU/Keyboard (Central Processing 
Unit), a 12 inch diagonal video display, and a cassette 
recorder for loading and recording programs and data. This 
basic system was enhanced with a Microtek 32K RAM memory 
expansion unit and a ESF-80 (Exatron Stringy Floppy) mass 
storage subsystem which was used in lieu of the cassette 
recorder. The entire system cost less than $1,400. 

This system proved to be reliable and provided the 
capabilities necessary for the successful completion of this 
project. A major enhancement, however, would have been 
at least one floppy disk system for loading programs and 
for input/output of data files. This addition would have 
significantly increased the speed of program loading and 
data file handling and provided the potential of alternate 
programming languages. 

C. PROGRAMMING LANGUAGE 

The programming language for this package is TRS-80 
Level II BASIC developed by Microsoft Inc. A Level II 
BASIC interpreter is resident in the ROM (Read Only Memory) 
contained in the TRS-80 CPU. This "built in" language was 
used as an expediency; the comparatively slow execution 
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speed of BASIC is not a significant factor for problems of 
the size this package is designed to handle. 

The native Level II BASIC was enhanced with several 
system commands made available through the ROM contained in 
the ESF-80 mass storage subsystem. These commands accom¬ 
plish several fundamental operations. The QSAVE and @LOAD 
commands are used to record and load program code. The 
0OPEN, 0CLOSE, SPRINT and ©INPUT commands are for data file 
handling. Finally, inclusion of the ©LOAD command in pro¬ 
gram code enables the chaining of programs (the ability to 
load a program while preserving the variables created by 
the previously loaded programs). 

D. PACKAGE PORTABILITY 

Without at least some minor modifications this package 
is portable only to an identical system. Virtually all 
DOS's (Disk Operating Systems) designed for the TRS-80 
Model I and Model III microcomputers, however, contain as a 
subset system commands which perform the same functions as 
those provided by the ESF-80 subsystem. Modification of 
these commands, which are contained in perhaps six lines of 
code in each of the six programs in the package, would make 
it portable to any TRS-80 Model I or Model III 48K RAM 
system equipped with at least one disk drive system. 

The most significant differences between Level II 
BASIC and the versions of BASIC implemented on other 
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microcomputers are keywords relating to display output. 
Because of the large amount of code directly relating to 
display output, it would require a sizeable effort to mod¬ 
ify this package for another microcomputer. The program 
logic, however, could serve as a model for programs written 
in other versions of BASIC or even FORTRAN. A partial 
translation of this package into FORTRAN-IV was accom¬ 
plished for implementation on a PDP-11-60 minicomputer 
system. 

E. LIMITATIONS 

Three interrelated constraints of this package are 
accuracy, speed and capacity. TRS-80 Level II BASIC 
includes the capability of defining double precision var¬ 
iables, thus significantly increasing the accuracy of 
computations. Double precision variables, however, require 
increased memory for storage and computation, with a resul¬ 
tant decrease in the capacity of the computer to handle 
large data arrays, and an increase in computation time. On 
the other hand, when variables are defined as single pre¬ 
cision some algorithms which use approximations requiring 
increased precision must be rejected in favor of others 
which are often slower. The goal then was to strike a 
proper balance so that a good degree of accuracy could be 
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attained and, at the some time, allow for reasonable speed 
of execution and capacity. 

The majority of the variables in this package are 
single precision variables, but integer variables were used 
where possible and a few procedures require the use of 
string, i.e. character or literal, variables. No double 
precision variables are used. Over the range of problems 
tested, the accuracy was found to be good (four or more 
decimal place accuracy) in most instances and fair (two or 
three decimal place accuracy) in a few cases. 

The capacity of the package is limited to a vector of 
200 observations for those procedures which apply to uni¬ 
variate data, i.e. the first four programs. The multiple 
linear regression program is limited to a data array of 
10X100 (response variable, 9 carrier variables and 100 
observations). The ANOVA program is limited to 10X20, 
10X10X10, and 10X10X10 arrays respectively for the one, two 
and three-way classification cases. The size of the data 
arrays handled by each program could be significantly 
increased if each program were loaded and executed as a 
single program. In order to chain the programs together 
and therefore pass variables back and forth between pro¬ 
grams, the array sizes are limited to those cited above. 

The package was designed for analysis of relatively 
small data arrays. Even so, due to the inherent nature of 
an interpreted language such as TRS-80 Level II BASIC, 
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subroutines which require many iterative calculations, e.g. 
several of the inverse probability distributions and some of 
the robust procedures, can require up to several minutes of 
execution time. Each algorithm was designed to execute as 
fast as possible, but whenever a choice had to be made 
between speed of execution and accuracy, the more accurate 
algorithm was selected. 

It is perhaps indicative of the current high degree of 
computational power available that an algorithm is consid¬ 
ered slow when it requires several minutes of execution 
time. A decade ago analysis of a moderate size sample 
which did not merit very expensive computer time could have 
consumed hours of tedious hand calculations. 

F. USER FRIENDLY FEATURES 

This package was designed to be used by the individual 
with a general understanding of statistical analysis pro¬ 
cedures and terminology. At the same time, it was recog¬ 
nized that there does not exist a completly standardized 
set of conventions for symbols and terms. While no attempt 
was made to define all procedures, parameters or terms, a 
brief explanation was provided where, in the subjective 
opinion of the author, it was felt confusion might exist. 
For example terms such as mean, variance, confidence level, 
sum of squares, F-statistic, etc. were assumed to be well 
understood by any potential user. When the user is 
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initially introduced to the term "exponential parameter 
lambda", however, he is advised that the term refers to the 
reciprocal of the exponential mean. 

Even the most knowledgeable user can inadvertently 
enter a value or character which would cause the program to 
terminate due to error or at least return an incorrect 
answer. Whenever the user is prompted for an input the 
value or character entered is first tested to see that it 
has some logical meaning within the context of the opera¬ 
tion or calculation to be executed. If it is recognized 
that it does not then the user is so advised and is allowed 
to reenter the value or character. If, for example, a data 
transformation or a data plot procedure is selected after 
summary statistics have been entered, the user is advised 
that the selected subroutine is not applicable since data 
have not been entered and he is prompted for another 
selection. 
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III. DISCUSSION OF PROGRAMS/ALGORITHMS 


A. STRUCTURE AND ORGANIZATION 

The package is divided into six chained programs. The 
control module is integrated into program 1 and must always 
be loaded first. The control module serves to define, 
dimension, and initialize variables for the first four 
programs and positions the variable pointer so that subse¬ 
quent programs can be loaded without loss of variables. 
The control module also contains the program instructions 
and program selection menu. 

After initial execution of the control module any of 
the six programs may be selected by the user. A video 
prompt tells the user which program tape to load into the 
ESF-80; the program is then automatically loaded and begins 
execution (auto run). The user is then free to move from 
one program to the next with the following limitation: when 
either program 5 or program 6 is loaded all variables are 
cleared and redimensioned in order to accomodate the mul¬ 
tidimensional arrays required by these programs. If the 
user desires to access one of the first four programs after 
using program 5 or program 6, it is necessary to first load 
program 1 and thus the control module in order to reini¬ 
tialize variables. 
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B. PROGRAM 1 


Program 1 contains data input/output procedures, data 
transformations and a number of common probability distribu¬ 
tions and inverse probability distributions. A detailed 
block diagram of program 1 is given in figure 2. 

1. Data Input/Output 

Data Input/Output (I/O) for the first four pro¬ 
grams is accomplished in program 1. Programs 5 and 6 con¬ 
tain their own data I/O subroutines. A one dimensional array 
of up to 200 data points may be manually entered via the 
keyboard, or a previously saved (recorded) data array may be 
loaded via the ESF-80. As a third alternative, the user may 
select one of four random number generators (normal, 
exponential, uniform and Bernoulli) to create data. Finally, 
the user is allowed to enter summary statistics only. This 
last capability is duplicated in program 3 to facilitate 
testing of various hypotheses. 

Any time data or summary statistics are entered, all 
previously entered data is cleared. The user is therefore 
given the option to save data on tape in the event that 
further analysis is desired. 

All data are stored in two single dimensional 
arrays, X(i) and Y(i). Both are defined as single preci¬ 
sion real arrays which are stored with seven digits of 
precision but are printed with a maximum of six digits of 
precision (if the seventh digit is 5 or greater then it is 
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BLOCK DIAGRAM OF PROGRAM #1 
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rounded up). All seven digits of precision are used for 
computations. The X array is used to save the raw data while 
the Y array may contain transformed data. All subroutines in 
the first four programs operate on the Y array. 

After the data has been entered either manually via 
the keyboard, loaded via a prerecorded data tape or 
generated by one of the four random number generators, two 
operations on the arrays are automatically executed. First, 
unless the data are already ordered from low value to high 
value, the Shell Sort subroutine is called to order the 
data. After the data are sorted the following sample 
statistics are calculated and saved for later use: sample 
mean, sample variance (unbiased estimate of population 
variance), maximum likelihood estimate of variance, skewness 
and kurtosis. 

2. Data Transformations 

McNeil [Ref. 3] claims that data transformations are 
the most powerful tools avialable to the data analyst. Such 
adamant statements inevitably invite arguments but it is at 
least true that, when properly applied, data trans¬ 
formations can be very useful to the analyst. They can 
reveal structure in the data that is otherwise obscured and 
can reexpress the data so that assumptions such as normal¬ 
ity are more nearly satisfied. The following data 
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transformations are certainly not exhaustive but are three 
common transformations which have proved useful in data 
analysis work. 

The square root, cube root and natural logarithm 
data transformations are optional after the raw data has 
been entered and sorted. The original data may be restored 
to the Y array from the X array by returning to the data 
transformations subroutine and selecting the 'No Transfor¬ 
mation* option. After each data transformation sample sta¬ 
tistics are recalculated. 

3. Probability Distributions 

The probability distributions menu contains five 
univariate continuous probability distributions, i.e. the 
normal, exponential, chi-square, Student*s-t and the F 
distributions and three discrete probability distributions, 
the binomial, negative binomial and Poisson distributions. 
The discrete distributions calculate both point probability 
and cumulative probability while the continuous distribu¬ 
tions calculate the left and right tail cumulative prob¬ 
abilities. The specific algorithms utilized are contained in 
the program listings at the end of this paper. 

Kennedy and Gentle [Ref. 4] point out that 'rough* 
approximations which provide two or three figure accuracy 
over certain ranges of parameter values would be adequate if 
one were only interested in probabilities. There is a need 
for greater accuracy in statistical analysis packages. 
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however, since these values are used in evaluating other 
expressions for which accurate results are required. A 
number of the probability distribution algorithms contained 
in program 1 are used as subroutines by other probability 
distributions, inverse probability distributions and, in 
later programs, used for evaluations such as confidence 
intervals and tests of hypotheses. Some of these applica¬ 
tions will be discussed in following paragraphs. 

In order to achieve the highest degree of accuracy 
possible, several of the probability distributions combine 
two or more algorithms. It is common that the faster 
algorithms work well over certain parameter ranges but lose 
accuracy quickly outside these ranges. It thus became 
necessary to utilize slower, often iterative algorithms in 
order to maintain an acceptable level of accuracy. Appendix 
A shows comparisons of probabilities calculated by the 
algorithms contained in this package and tabled values given 
in Handbook of Tables for Probability and Statistics [Ref. 
5]. 

a. Normal Distribution 

The central limit theorem simply stated says 
that the sample mean of a random sample drawn from any 
distribution with finite variance is, in the limit, norm¬ 
ally distributed about the true population mean. This 
theorem makes the normal distribution one of the most 
important distributions in statistical analysis. A number of 
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the statistical inference procedures contained in this 
package have the normal distribution as their foundation and 
several other probability distribution subroutines utilize 
relationships with the normal distribution. 

The algorithm for the normal distribution is 
based on a polynomial expansion of the right tail cumu¬ 
lative probability given in Abramowitz and Stegun [Ref.6]. 
This algorithm gives results which are exact when compared 
to tabled values (see Appendix A). The nominal computation 
time is one second. 

b. Exponential Distribution 

The exponential distribution is commonly used to 
represent the distribution of time between the occurrences 
of certain events. One of the most frequent applications is 
in relaibilty theory. 

The exponential cumulative probability is 
particularly easy and fast to calculate. The closed-form 
expression for the right tail cumulative probability is 
evaluated and the resulting probabilities are exact to 
single precision accuracy. 

c. Chi-square Distribution 

The importance of the chi-square distribution in 
statistical analysis is based on the theorem that the sum of 
squares of k independent and identically distributed 
standard normal random variables has a chi-square 
distribution with k degrees of freedom. Extensions of this 
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theorem serve to test hypotheses concerning and establish 
confidence intervals for the population variance under the 
assumption of normality. 

Two different algorithms were utilized to 
approximate chi-square cumulative probabilities. A rela¬ 
tionship between the chi-square distribution and the normal 
distribution which is given in Abramowitz and Stegun [Ref. 
7] was used for the large degrees of freedom case (degrees 
of freedom greater than thirty). This algorithm is fast 
(approximately one second) and is accurate to a minimum of 
three decimal places. For smaller degrees of freedom, 
however, the accuracy of this approximation falls off rap¬ 
idly. Poole and Brochers [Ref. 8] give a BASIC program based 
on a series expansion of the right tail cumulative 
probability. This algorithm is slower (six seconds are 
required for thirty degrees of freedom) but maintains at 
least three decimal place accuracy down to one degree of 
freedom. 

d. Student’s - t Distribution 

The t distribution, like the chi-square distri¬ 
bution, has wide application in statistical analysis work. 
By definition, if a standard normal random variable is 
divided by the square root of a chi-square random variable 
over its degrees of freedom then the resulting random 
variable has a t distribution with the same degrees of 
freedom. It can be shown that if a random sample of size n 
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is drawn from a normal distribution then the square root of 
the sample size times the difference between the sample mean 
and true population mean all divided by the square root of 
the sample variance has a t distribution with n-1 degrees of 
freedom. An important point is that this t statistic does 
not depend on the true population variance. 

Two algorithms were utilized to obtain 
t-distribution cumulative probabilities. Abramowitz and 
Stegun [Ref. 9] give a distributional relationship with the 
normal distribution which produces good results (three to 
four decimal place accuracy in one second) when degrees of 
freedom exceed thirty and a series expansion of the 
cumulative distribution function for smaller degrees of 
freedom. The latter algorithm returns four or more decimal 
place accuracy in one to three seconds, 
e. F Distribution 

The F distribution is still another sampling 
distribution of great importance in statistical analysis 
work. A random variable which is the ratio of two inde¬ 
pendent chi-square random variables each divided by its 
degrees of freedom has an F distribution with the respective 
numerator and denominator degrees of freedom. This statistic 
is often referred to as the F ratio and serves, for example, 
to test hypotheses concerning differences between variances 
of independent normal populations. 
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Kennedy and Gentle [Ref.10] point out that 
existing algorithms for F distribution probabilities do not 
give a high degree of accuracy. This reference does, how¬ 
ever, give a FORTRAN program based on a series expansion of 
a basic integral form of the left tail cumulative prob¬ 
ability which is among the best available algorithms. This 
package utilizes a direct translation of the FORTRAN pro¬ 
gram given in the above reference. The algorithm generally 
produces three decimal place accuracy and occasionally four 
decimal place accuracy. Computation time is a function of 
degrees of freedom and ranges from a nominal one second for 
small degrees of freedom to appoximately seven seconds for 
120 numerator and denominator degrees of freedom, 
f. Discrete Distributions 

The binomial, negative binomial and Poisson 
distribution algorithms perform direct calculations of their 
respective probability functions to produce point 
probabilities. Cumulative probabilities are calculated by 
summation of point probabilities with one exception. The 
binomial cumulative probability utilizes the Poisson 
approximation if the number of trials is greater than 45 and 
the probability of occurrence on a single trial is either 
less than or equal to .1 or greater than or equal to .9. The 
normal approximation is used if the number of trials is 
greater than 45 and the probability of occurrence on a 
single trial times the number of trials is greater than 5. 
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All point probabilities are accurate to four or 
more decimal places. The normal and Poisson approximations 
to the binomial cumulative probability are accurate to three 
or more decimal places; the other cumulative probability 
procedures produce a minimum of four decimal place accuracy. 

4. Inverse Probability Distributions 

Exponential percentiles are calculated with the same 
ease as exponential probabilities and good closed form 
approximations for normal and chi-square percentiles were 
found. Abramowitz and Stegun [Ref. 11] was the source of the 
normal percentiles approximation while the algorithm for 
chi-square percentiles is a direct translation of FORTRAN 
program listed in "Collected Algorithms from ACM" [Ref. 12], 
These algorithms give accurate results when compared to 
tabled values (see Appendix 2) and all execute in one to 
three seconds. 

It was more difficult to obtain accurate results for 
t-distribution and F-distribution percentiles. Abramowitz 
and Stegun [Ref. 13] give closed form approximations for 
both distributions which perform well for larger degrees of 


freedom (ten 

or greater 

for 

the 

t-distribution 

and 

twenty/twenty 

or greater 

for 

the 

F-distribution). 

These 


approximations, however, lose accuracy rapidly for smaller 
degrees of freedom. In these cases it was necessary to 
resort to an algorithm which starts with a rough 
approximation of the percentile and, calling the respective 
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probability distribution subroutine, iteratively refines the 
approximation until the returned probability converges to 
the specified value. This algorithm produces accurate 
percentiles but computation time can extend up to several 
minutes. 

C. PROGRAM 2 

Program 2 contains the six subroutines which are 
classified as exploratory data analysis procedures. Their 
primary purpose is to provide information about the data 
(raw or transformed) which reveals patterns, trends, spread, 
symmetry, how closely the data follows some classical 
distribution, etc. The information about the data is 
presented in two basic formats. The first two subroutines in 
this program present a numerical summary of the data while 
the last four subroutines give graphical representations 
(plots) of the data. 

1. Numerical Summaries 

The first numerical summary of the data displays 
several classical sample statistics, each of which is a 
measure of some facet of the data's structure. A measure of 
the central location and/or typical value of the data is 
given by the sample mean. Two measures of the spread of the 
data are given. The first is the sample variance, i.e. the 
unbiased estimate of population variance under the 
assumption of normality, and the second is the maximum 
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likelihood estimate of population variance under the same 
assumption. The sample statistic called skewness provides a 
measure of the symmetry of the data. A value of zero 
indicates a symmetric distribution while positive and 
negative values indicate right and left skewed distribu¬ 
tions respectively. Finally, kurtosis provides a measure of 
how much of the data is distributed into the tail of the 
distribution. For reference, a normal distribution has 
kurtosis of zero while the exponential distribution, which 
has a longer tail than the normal distribution, has kur¬ 
tosis of six. 

The second numerical summary of the data is the 
seven number summary. The "seven" refers to the basic 
elements of the display. These consist of the median, upper 
and lower quartiles (75 and 25 percentiles), the upper and 
lower eighths (87.5 and 12.5 percentiles ) and the extremes 
(maximum and minimum). 

Like the mean, the median gives a measure of the 
central or typical value of the data. The three remaining 
pairs of numbers give a measure of the spread of the data. 
Two other numbers associated with each of these pairs give a 
measure of skewness of the data. The first is simply the 
midspread or average of each of the pairs which can be 
compared with the median and with each other. The second 
measure of skewness makes this comparison a little more 
automatic. If M represents the median, Q1 the lower 
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quartile, Qu the upper quartile then (((Qu-Ql)/2)-M)/(Qu-Ql) 
will be zero if the quartiles are symmetric about the median 
while a positive or negative number will suggest the 
corresponding direction of skewness. A similar number is 
calculated for eighths and for extremes. 

2. Data Plots 

The histogram is a graphical representation of the 
data in which each point is assigned to a class interval 
depending upon its value. The number of data points in a 
class interval determine the height of the graphical bar 
representing that cell. The resulting display is then a 
picture of the distribution which generally approximates the 
shape of the underlying probability distribution function. 

An important parameter in determining the usefulness 
of the histogram is the width of the class interval. The 
extreme cases illustrate this point. If the class interval 
is so wide that all data points fall into a single cell then 
about the only information provided is the number of data 
points in the sample. On the other hand, if the class 
interval in so narrow that every data point belongs to a 
separate cell then the picture is very difficult to 
interpret. The histogram subroutine contained in this 
program can have up to fifteen cells. The initial class 
interval is automatically determined so that there will be a 
range of at least five cells from the minimum to the maximum 
data point. Some of these cells may be empty. After the 
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initial histogram is drawn the user is given the option to 
increase the class interval or, subject to the fifteen cell 
limitation, decrease the class interval. 

A graphical representation of the data similar to 
but containing more information than the histogram is the 
stem and leaf plot. The significant difference between the 
two plots is that the incremental elements in each cell of 
the stem and leaf, i.e. the leaves, are numbers instead of 
an incremental portion of a solid bar. The most significant 
digit of each data point is represented by the class 
interval labels running along the stem (vertical axis) of 
the plot and the leaves are the next most significant digit 
of each data point. 

McNeil [Ref. 14] suggests several advantages the 
stem and leaf enjoys over the histogram. As well as pro¬ 
viding the same information about the shape of the distri¬ 
bution it reveals intracell distribution. It is also pos¬ 
sible to obtain reasonable estimates of percentiles simply 
by counting leaves. 

As with the histogram, the initial class interval is 
automatically determined by the computer. After the first 
plot the user is again allowed to redraw the plot with a 
different class interval, this time subject to a limitation 
of sixteen cells. 

A fairly simple procedure for testing the conformity 
of the data to some specified theoretical probability 
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distribution is to plot the ordered sample against the 
corresponding theoretical percentiles. A one to one 
correspondence would, of course, plot as a straight line. 
Systematic departure from a straight line plot indicates 
that the sample does not conform well to the theoretical 
distribution. 

This procedure was implemented for the normal and 
exponential distributions and these two plots are the final 
graphical representations of the data contained in this 
program. These plots perform well for fairly significant 
departures from the theoretical distributions. Due to the 
relatively low resolution graphics capability of the TRS-80 
video display, however, slight to moderate departures from 
the theoretical distributions are difficult to detect. 

D. PROGRAM 3 

The statistical analysis procedures contained in program 
3 are classified as parametric statistical inference proce¬ 
dures. They include classical tests of hypotheses and cal¬ 
culation of confidence intervals for normal, exponential and 
Bernoulli population parameters. The specifics of these 
classical procedures are given in almost all first course 
statistics textbooks and will not be reiterated here. The 
discussion will be limited to the basic format of the 
procedures. Major portions of each of the algorithms 
contained in this program have already been discussed, i.e. 


36 



the probability distributions and inverse distributions 
which were developed for program 1. A detailed block dia¬ 
gram of the hypotheses tested and confidence intervals 
calculated by this program is given in Figure 3. 

1. Tests of Hypotheses 

The basic format for the various test of hypotheses 
is as follows. Once the user has selected the assumed 
distribution and the specific null hypothesis to be tested 
then the appropriate test statistic (usually a combination 
of sample statistics already calculated in program 1) is 
calculated and displayed. Using this test statistic, the 
appropriate probability distribution subroutine is called 
and the probability at which this statistic will become 
statistically significant (P-value) is returned and dis¬ 
played . The user is then prompted to enter the probability 
of type I error (alpha value) he desires to establish for 
the test. If alpha is less than the P-value then the null 
hypothesis is accepted, otherwise it is rejected. 

A more classical approach to hypothesis testing is 
to calculate a value for a statistically significant sta¬ 
tistic given the alpha value and then compare the actual 
test statistic calculated from the sample to this value to 
determine acceptance or rejection. This approach was re¬ 
jected in favor of the one given above for two reasons. 
First the accuracy and speed of the probability distribu¬ 
tion subroutines are generally better than those for the 


37 



BLOCK DIAGRAM OF PROGRAM It 3 


I PROGRAM *31 




7~1 

1 EXPONENTIAL SWLEl 

BERNOULLI SAWLEl 

1) Hypotheses Concerning Mean 

a) MU = HUo -vs- HU =/= HUo 

b) HU <= HUo -vs- HU > HUo 

c) MU >= HUo -vs- HU < HUo 

2) Confidence Interval far Mean 

1) Hypotheses Concerning Proportion 

a) P = Po -vs- P =/= Po 

b) P <= Po -vs- P > Po 

c) P >= Po -vs- P < Po 

2) Confidence Int. for Proportion 


NORMAL SAfflfl 

1) Hypotheses Concerning Mean 

3) HU = HUo -vs- MU =/= HUo 

b) HU <= HUo -vs- HU > HUo 

c) HU >= HUo -vs- HU < HUo 

2) Hypotheses Concerning Variance 

a) VAR = mo -vs- m =/= mo 

b) VAR <= VARo -VS- VAR > VARo 

c) m >= VARo -vs- VAR < VARo 

3) Confidence Interval for Mean 

1) Confidence Interval for Variance 


FIGURE 3 
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inverse distribution subroutines. 


Calculation of the 


P-value uses the probability distribution subroutines while 
calculation of a statistically significant value for the 
statistic uses the inverse probalility subroutines. Second, 
and perhaps most importantly, the displayed P-value provides 
more information than a simple accept/reject criteria. 

2. Confidence Intervals 

The format for the confidence interval subroutines 
is similar to that utilized for tests of hypotheses. The 
user selects the assumed population and the desired 
confidence level. After this information has been entered 
the appropriate inverse probability distribution is called 
and returns the value at which a corresponding sample sta¬ 
tistic would become significant. A confidence bound for the 
parameter of interest is calculated from this value and 
displayed. This algorithm is first utilized to calculate 
two sided confidence limits and then a upper and lower one 
sided confidence limit. 

3. PROGRAM 4 

This program contains two robust and one nonparametric 
statistical inference procedures. 

1. Robust-t 

When it is suspected that the random sample contains 
outliers, then it is often desirable to obtain an estimate 
for the population mean that places less weight on the 
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extreme points. Two techniques for obtaining such an 
estimate are the Winsorized-t and trimmed-t procedures. 
Winsorization is accomplished by symmetrically setting one 
or more of the extreme values of the ordered sample equal 
to its nearest neighbor. Trimming, on the other hand, omits 
the symmetrically placed extreme values. It should be noted 
that in both cases there is an implicit assumption that the 
true distributional form is symmetric. The depth to which 
the extreme values are Winsorized/trimmed is referred to as 
the level of Winsorization/trimming. After the data are 
Winsorized/trimmed the point estimate of the population mean 
is calculated in the usual manner. Dixon and Massey [Ref. 
15] point out that there is little loss of efficiency if the 
data are Winsorized (at mild levels) when the extreme points 
are actually valid and that much improved estimates result 
when the extremes are not valid observations. Trimmed 
estimates are somewhat less efficient. 

It has also been shown that the t-ratio which is 
calculated using the sample variance of the Winsorized data, 
multiplied by an appropriate constant (see Dixon and Massey 
[Ref. 16] for details) closely follows the t-distribution 
with h-1 degrees of freedom where h is the number of 
obervations not Winsorized/trimmed. This statistic is used 
to calculate confidence bounds for the true population mean. 

This program allows Winsorizing/trimming of the data 
up to the tenth level. Computation time is a function of 
the selected level but is not excessive in any case. 
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2. Jackknife 


If the distribution from which a random sample is 
drawn is normal then a confidence interval for the true 
population variance is properly based on the chi-square 
sampling distribution. Such confidence intervals are, 
however, very sensitive to the underlying distribution and 
do not perform well when the assumption of normality is 
violated. A robust method for handling departures from 
normality is the jackknife technique. The jackknife 
technique is discussed in some detail in Mosteller and Tukey 
[Ref. 17]. 

This procedure, more than any other in this package, 
shows up the speed disadvantage of an interpreted language 
like TRS-80 Level II BASIC. The crux of this procedure is 
to divide the data into g groups of size h and then 
iteratively calculate the sample variance g times, each time 
leaving out one group of h. The usual method is to let h=l 
and therefore g=n where n is the size of the sample. This 
eliminates any arbitrariness in the formation of the groups 
but means that a sample of one-hundred data points requires 
one-hundred calculations of the sample variance. This 
specific example required 25 minutes of calculation time. 
As a partial solution to this problem, the user is allowed 
to set h to any integer between one and five so long as the 
data size is evenly divisible by h. When h is greater than 
one the ordered data must be randomized to avoid 
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systematically leaving out the h smallest values first, then 
the next h smallest values, etc. Randomizing the data does 
introduce some of the arbitrariness mentioned earlier and 
results in slightly different answers each time the 
procedure is repeated. A number of sample problems, however, 
displayed reasonable stability from one calculation to the 
next. Application of the jackknife technique to a sample of 
one hundred observations divided into twenty groups of size 
5 required approximately five minutes of calculation time. 

3. Confidence Interval for the Median 

By the definition of the median, the probability 
that any random observation drawn from any continuous 
population is less than or equal to the median equals 
one-half. Based on this fact it can be shown that the 
probability that the median lies between the ith and the jth 
order statistic of a random sample has a binomial 
distribution with parameters p=.5 and n equal to the sample 
size. Since no assumptions are made about the population 
from which the random sample is drawn other than that it is 
continuous, it is possible to obtain a truly nonparametric 
confidence interval for the population median. 

Confidence intervals for the median calculated by 
this program are in exact agreement with tabled values 
contained in Dixon and Massey [Ref.18]. Calculation time 
varies from a few seconds to several minutes depending on 
the sample size. 
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F. PROGRAM 5 


Program 5 contains a multiple linear regression 
procedure which, by the method of Ordinary Least Squares 
(OLS), estimates the coefficients of the basic linear model, 
Y = X B + e, where Y is the vector of observations of the 
dependent variable, X is the matrix of independent 
variables, B is the vector of coefficients estimated and e 
is the error vector. The program also calculates and 
displays a number of the standard statistics each of which 
measures some facet of the quality of the OLS solution. 

The computational method used by this program is 
straightforward matrix algebra applied to the basic linear 
model. The greatest difficulty with this method is in 
obtaining accurate results when the X'X matrix is inverted. 
According to Draper and Smith [Ref. 19] one of the prime 
causes of roundoff error during this operation is the 
presence of widely differing orders of magnitude in the 
numbers involved. In order to minimize this problem the 
data are normalized or centered as described in the 
aforementioned reference. Two additional advantages result 
from normalizing the data. If Z represents the normalized X 
matrix then Z'Z is the correlation matrix and can be 
displayed to provide the analyst with additional information 
about the data. The coefficients which are directly 
estimated using the normalized data are referred to as the 
beta coefficients. The beta coefficients have the 
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advantage of being independent of the units of measurement 
of the original data. This program displays both the 
standard regression equation coefficients and the beta 
coefficients. 

Other standard statistics calculated and displayed are 
the ANOVA table statistics, R-square, adjusted R-square, 
standard error of estimate, partial-F for the coefficients 
and standard error of the coefficients. Each of these 
statistics with the exception of the adjusted R-square 
statistic is discussed in detail in Draper and Smith [Ref. 
20]. The adjusted R-square statistic is discussed in 
Intriligator [Ref. 21]. 

A perhaps not so standard feature of this program is a 
display of the diagonal elements of the hat matrix. The hat 
matrix is defined as X(X*X)“ 1 X’, where these symbols have 
the same meaning as given earlier. The diagonal elements 
of this matrix give a measure of the influence or leverage 
of the ith observation of the response variable on the 
corresponding fitted value. Large values are indicative of 
high leverage observations and may indicate that these 
observations are outliers. Hoaglin and Welsch [Ref. 22] 
give as a rule of thumb for ’large* values, any value in 
excess of twice the rank of the hat matrix divided by the 
number of observations. All diagonal elements of the hat 
matrix are displayed and large values are highlighted. 
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This program allows the user to regress the dependent 
variable on all or any subset of the independent variables. 
It is therefore possible for the user to use various 
algorithms for selecting the 'best subset' of carriers. 

The final feature of this program is a plot of the 
estimates of the response variable against the residuals. 
According to McNeil [Ref. 23] and other recognized 
authorities in the field of data analysis, such a display of 
the residuals is one of the best tools in existence for 
obtaining a measure of the overall fit of the data to the 
model. 

The regression procedure contained in this program was 
tested against the Portland Cement example contained in 
Draper and Smith [Ref. 24]. The data are thirteen 
observations of a response variable and four carrier 
variables. This regression was completed in 48 seconds. 
The estimated regression equation coefficients and all 
statistics agreed with the reference to a minimum of three 
decimal places and in most cases four or more decimal 
places. 

G. PROGRAM 6 

Program 6 is the last program in this package and 
contains ANalysis Of VAriance (ANOVA) procedures. The 
algorithms for the ANOVA procedures were taken from Duncan 
[Ref. 25] who makes the claim that ANOVA is one of the most 
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powerful tools in statistical analysis. He further states 
that the basic procedure of ANOVA consists of classifying 
and cross-classifying some statistical result of interest to 
determine if a specified classification is statistically 
important in affecting the result. A potential application 
of ANOVA was given in an earlier section of this paper. 

It was attempted to include as much flexibility in the 
ANOVA procedures contained in this program as feasible. 
One, two and three-way classification of the statistical 
result of interest is allowed. The one-way classification 
procedure allows up to ten levels of the classification 
(ten ranks in the example presented earlier) and up to 
twenty observations at each level. The program also 
accomodates an unequal number of observations at each level. 
The two-way classification ANOVA procedure allows ten levels 
of each classification and up to ten replications. The 
three-way classification again allows up to ten levels of 
each of the three classifications but due to computer memory 
limitations does not allow for replications. The two- and 
three-way classification procedures calculate interaction 
statistics as well as the main classification statistics. 
All three procedures are accurate to a minimum of three 
decimal places. 
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IV. CONCLUSION 


When microcomputers were first introduced only a few 
years ago, they were widely perceived as interesting toys 
used primarily for video games. This perception is rapidly 
changing and today microcomputers are increasingly accepted 
as serious tools for many small business applications. Just 
as the military services have adopted electric typewriters, 
photo-copy machines and other modern office machinery, it is 
not hard to imagine that relatively inexpensive microcomput¬ 
ers will also be adopted and find their way to all levels of 
command. A popular commercial application of microcomputers 
is in the area of wordprocessing. It may well be that this 
will be the initial application of microcomputers in the 
military services. Whatever the intended application might 
be, the end result will be ready access to a high degree 
of computational power. 

If this potential ready access to computational power is 
to be put to use by the military officer trained in 
specialized areas such as operations research, then it is up 
to these officers to develop the specialized software. 
Specialized software for microcomputers has been the subject 
of this thesis. It has presented a software package 
containing some of the basic statistical analysis tools 
required for effective analysis work. Many areas of interest 
to the trained analyst have not been addressed. Other recent 
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efforts in this direction are addressed by Morgeson [Ref. 
26] who developed a similar package for the Apple II Plus 
microcomputer and Duff [Ref. 27] who has developed a network 
optimization package for the same microcomputer. Additional 
efforts in the development of statistical analysis software 
are highly recommended and encouraged. 
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APPENDIX A 


ACCURACY COMPARISONS FOR PROBABILITIES 
NORMAL DISTRIBUTION 

PROBABILITY 


STATISTIC 

APPROXIMATION 

TABLED VALUE 

0.00 

.5000 

.5000 

0.25 

.5987 

.5987 

0.50 

.6915 

.6915 

0.75 

.7734 

.7734 

1.00 

.8413 

.8413 

1.25 

.8944 

.8944 

1.50 

.9332 

.9332 

1.75 

.9599 

.9599 

2.00 

.9772 

.9772 

2.25 

.9878 

.9878 

2.50 

.9938 

.9938 

2.75 

.9970 

.9970 

3.00 

.9987 

.9987 

3.25 

.9994 

.9994 

3.50 

.9998 

.9998 

3.75 

.9999 

.9999 

4.00 

1.0000 

1.0000 
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CHI-SQUARE DISTRIBUTION 


PROBABILITY 


STATISTIC 

D.F. 

APPROXIMATION 

TABLED VALUE 

7.88 

1 

.9950 

.9950 

5.99 

2 

.9500 

.9500 

1.21 

3 

.2495 

.2500 

3.36 

4 

.5005 

.5000 

9.24 

5 

.9001 

.9000 

2.18 

8 

.0250 

.0250 

3.57 

12 

.0100 

.0100 

27.5 

15 

.9751 

.9750 

34.8 

18 

.9900 

.9900 

10.1 

19 

.0496 

.0500 

9.59 

20 

.0250 

.0250 

33.9 

22 

.9497 

.9500 

34.4 

25 

.9004 

.9000 

27.3 

28 

.4990 

.5000 

49.6 

29 

.9900 

.9900 

47.0 

30 

.9751 

.9750 
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STUDENT'S - t DISTRIBUTION 


PROBABILITY 


STATISTIC 

D.F. 

APPROXIMATION 

TABLED VALUE 

31.821 

1 

.9900 

.9900 

1.886 

2 

.9000 

.9000 

.271 

4 

.6001 

.6000 

3.365 

5 

.9900 

.9900 

1.943 

6 

.9500 

.9500 

.703 

9 

.7501 

.7500 

2.228 

10 

.9750 

.9750 

4.140 

14 

.9995 

.9995 

1.341 

15 

.9001 

.9000 

.257 

19 

.6000 

.6000 

2.086 

20 

.9750 

.9750 

2.492 

24 

.9900 

.9900 

3.725 

25 

.9995 

.9995 

1.311 

29 

.8999 

.9000 

2.750 

30 

.9949 

.9950 

1.684 

40 

.9500 

.9500 

1.296 

60 

.9000 

.9000 

3.373 

120 

.9995 

.9995 
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F - DISTRIBUTION 


PROBABILITY 


STATISTIC 

D.F./D.F. 

APPROXIMATION 

TABLED VALUE 

647.80 

1/1 

.9750 

.9750 

198.50 

1/2 

.9950 

.9950 

53.59 

3/1 

.9000 

.9000 

5.79 

2/5 

.9501 

.9500 

21.92 

4/6 

.9990 

.9990 

4.24 

5/10 

.9751 

.9750 

5.76 

6/12 

.9950 

.9950 

1.93 

8/25 

.9001 

.9000 

2.16 

10/30 

.9495 

.9500 

3.64 

12/40 

.9990 

.9990 

2.06 

15/60 

.9749 

.9750 

2.09 

24/120 

.9950 

.9950 

5.23 

10/3 

.9000 

.9000 

46.76 

15/4 

.9990 

.9990 

5.17 

20/6 

.9750 

.9750 

2.63 

30/30 

.9950 

.9950 

63.06 

120/1 

.9000 

.9000 
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APPENDIX 3 


ACCURACY COMPARISONS FOR PERCENTILES 
INVERSE NORMAL DISTRIBUTION 

PERCENTILE 


PROBABILITY 

APPROXIMATION 

TABLED VALUE 

.0005 

-3.29 

-3.29 

.0102 

-2.32 

-2.32 

.2514 

-0.67 

-0.67 

.4960 

-0.01 

-0.01 

.5199 

0.05 

0.05 

.6554 

0.40 

0.40 

.7549 

0.69 

0.69 

.9015 

1.29 

1.29 

.9251 

1.44 

1.44 

.9750 

1.96 

1.96 

.9901 

2.33 

2.33 

.9951 

2.58 

2.58 

.9990 

3.09 

3.09 

.9995 

3.29 

3.29 
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INVERSE CHI-SQUARE DISTRIBUTION 

PERCENTILE 

PROBABILITY D.F. APPROXIMATION TABLED VALUE 


500 

1 

.455 

.455 

950 

2 

5.99 

5.99 

025 

3 

.216 

.216 

990 

5 

15.1 

15.1 

005 

8 

1.34 

1.34 

975 

10 

20.5 

20.5 

010 

12 

3.57 

3.57 

900 

15 

22.3 

22.3 

995 

20 

40.0 

40.0 

750 

22 

26.0 

26.0 

025 

25 

13.1 

13.1 

975 

28 

44.5 

44.5 

100 

30 

20.6 

20.6 
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INVERSE STUDENT'S 


t DISTRIBUTION 


PERCENTILE 


PROBABILITY 

D.F. 

APPROXIMATION 

TABLED VALUE 

.9900 

1 

31.840 

31.821 

.9500 

2 

2.920 

2.920 

.0250 

3 

-3.182 

-3.182 

.7500 

5 

.726 

.727 

.9750 

10 

2.229 

2.228 

.9995 

15 

4.073 

4.073 

.0500 

18 

-1.734 

-1.734 

.9000 

20 

1.326 

1.325 

.6000 

25 

.256 

.256 

.2500 

28 

-.683 

-.683 

.9995 

30 

3.646 

3.646 

.0100 

40 

-2.424 

-2.423 

.0005 

60 

-3.460 

-3.460 

.9950 

120 

2.618 

2.617 
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INVERSE F - DISTRIBUTION 


PERCENTILE 


PROBABILITY 

D.F./D.F. 

APPROXIMATION 

TABLED VALUE 

.975 

1/1 

647.50 

647.80 

.995 

1/2 

198.10 

198.50 

.900 

3/1 

53.61 

53.59 

.950 

2/5 

5.79 

5.79 

.990 

4/6 

9.15 

9.15 

.975 

5/10 

4.24 

4.24 

.995 

6/12 

5.76 

5.76 

.900 

8/25 

1.93 

1.93 

.950 

10/30 

2.17 

2.16 

.990 

12/40 

2.67 

2.66 

.975 

15/60 

2.06 

2.06 

.995 

24/120 

2.09 

2.09 

.900 

10/3 

5.23 

5.23 

.975 

60/20 

2.23 

2.22 

.995 

30/30 

2.63 

2.63 

.900 

120/1 

63.34 

63.06 

.950 

120/120 

1.35 

1.35 


56 



COMPUTER LISTINGS 



10 REM INSTRUCTIONS/PROGRAM SELECTION SUBROUTINE 

20 IF K9O0 GOTO 1000ELSE CLEAR 
30 RANDOM 
40 CLEAR 330 

SO DIM X<201 ),Y<201),NI<15) 

60 DEFINT I-K 5 DEFSTR I 

69 LBS-CHHSU23) : RB*=CHRS(129> ! MB*»CHRS<176>: BLS-CHRK292) 

70 REM INSTRUCTIONS •** 

80 CLS S PRINT TAB(20>rSTATISTICAL ANALYSIS* : PRINT TAB(23>!*by R. P. Isbell* ! PRINT TABdlJrNaval Po 
stgraduate School, Monterey, Ca." 

90 IF VC9O0 GOTO 290 

100 PRINT 5 PRINT* This program ie the first of six ‘chained ' programs in* 

110 FRINT*this package. It contains the first set of statistical analy-* 

120 FHINT*sis subroutines and also serves as the control module for the* 

130 PRINT*package. This program must always be loaded first as it accom-’ 

140 PRINT'plishes variable definition, dimensions arrays and, most unpor-" 

190 FRINT*tantly, positions the variable pointer to allow subsequent * 

160 FRINT*programs to be Ipaded without loss of variables. If at any * 

170 PRINT'time it becomes necessary tD interupt program execution and re-* 

180 PRINT’tssue the RUN command, the action must be taken with Program #1* 

190 PRINT*U.e, this program) loaded,* 

200 PRINT : PRINT* Press the SPACE BAR to display additional information.*; 

210 IF INKEYS O * * GOTO 210 

220 CLS i PRINT* Data input for programs 61 thru 64 is accomplished in pro- gram #1 and is limited to a data 
vector of 200 observations. All data and variables are saved as you switch back and forth between the first fo 
ur programs,* 

230 PRINT SPRINT" When either program #9 or program 66 is loaded all ajrrent variables are cleared in order t 

o provide suffiaent memory far the multidimensional data arrays used by these programs. Any* 

240 FRINT*data which was entered in program 61 that you want to retain should be saved to tape prior to load! 
ng either of the last two programs. 

290 PRINT 5896,“PRESS SPACE BAR TO DISPLAT PROGRAM SELECTION MENUS* 

260 IF INKEYSO* * GOTO 260 

270 K9-1 S REM **• PROGRAM SELECTION MENU *«* 

280 CLS S PRINT TAB(20)5*PROGRAM SELECTION MENU* 

290 PRINT 5128,‘61 Data I/O, Transformations, Proto. & Inv. Proto. Distributions"SPRINT 

300 PHINT*62 Sample Stat's, 7-Num-Sum, Data Plots*:PRINT 

310 FRINT*63 Parametric Tests of Hypotheses & Confidence Intervals'SFRINT 

320 PRINT**4 Robust Estimates and Confidence Intervals’SFRINT 

330 PRINT*63 Multiple Linear Regression*SPRINT 

340 PRINT*6i Analysis of Variance* 

390 KF-896 S KS-6 ! GG5UB 900 

420 IF IO-l GOTO 1000 

430 IF 10-2 THEN IT-"l* 5 GOTO 490 

440 IF 10-3 OR 10-4 THEN ZT-*2* ELSE ZT»“3* 

490 IF 10-3 OR 10-9 THEN 10=1 ELSE 10-2 

460 CLS I PRINT-INSERT TAPE 6*3TJ* INTO ESF-80 AND PRESS SPACE BAR- 
470 IF INKEYSO* * GOTO 470 
480 0LOAD 10 

900 REM ♦ * * OPTION SELECTION SUBROUTINE •** 

910 PRINT 5KP*9,"OPTION DESIRED -SLBtSMBSSRBSS 
920 n-INXEYS S IF ZI-** THEN 520ELSE IO-VALOI) 
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530 IF I0<1 OR IO>KS THEN PRINT 0KP, BL*i! PRINT 0KP,-«** NOT A VALID OPTIOM ****!! FOR 1-1 TO 100 
0 ! NEXT ! PRINT OKP, BL»| ! GOTO 510 
540 RETURN 

600 PRINT 0KP*5,LB»!" *!L*r or “IRV,* “IRBSi* *;LB*,*MB*iRB«5 

610 20«INKEY* ! IF ZQ='* THEN 610ELSE PRINT 0KP+18.Z0}: FOR 1=1 TO 100 ! NEXT 

620 IF ZQOLS AND IQOCHR4(ASCO.*)-32) AND ZQOR* AND ZQOCHRKASCUWH32) THEN 600 

630 RETURN 


ESSCRA M 

1000 REM ••• PROGRAM t 1 
1010 REM MAIN MENU * * * 

1020 CLS ! PRINT TAB(20);*»** MENU W 

1030 PRINT 0128,-<1) Data I/O 

1040 PHINT"<2) Data Transformations 

1050 PRINT*(3> Probability Distributions 

1060 PHINT"(4) Inverse Probability Distributions 

1070 FKINT*(5) Return to Program Selection Menu 

1080 KP=704 ! KS*5 : GOSUB 500 

1150 IF 14=1 AND 10=2 OR 14=0 AND IN-0 AND 10=2 PRIMT07O4,CHR*(253);:PRINT 0704,THERE IS NO DATA T 
O TRANSFORMMFOR 1=1 TO 1500:NEXT'.PRINT 07O4,CHR»(253)iGOTO 1080 
1160 ON 10 GOSUB 2000, 3500, 5000, 5200, 270 
1170 GOTO 1020 


BA1A OfEST L a m m 


2000 REM DATA INPUT/OUTPUT SUBROUTINE 

2010 CLS ! PRINT 011,'*** DATA INPUT/OUTPUT#** 

2020 PRINTIFRINTIPRINT-MENU:* 


2030 PRINT- 
2040 PRINT- 
2030 PRINT- 
2060 PRINT* 

o*»s) 


(1) Input Data via Keyboard 

(2) Input Data via Taoe 

(31 Input Summary Statistics ONLY 

(4) Create a random sample using random number generator* 


PRINT- 


2070 PRINT- (5) SAVE data to tape 

2080 KP-768 ! KS=3 ! GOSUB 500 

2100 IF 10=1 OR 10=2 THEN 15=1 ELSE 15=0 

2120 ON 10 GOTO 2200, 2700, 2900, 4000, 3300 


(for demonstration purp 


mm DATA 

2200 REM INPUT DATA VIA KEYBOARD ••• 

2210 CLS 114-0 

2220 PRINTTo end data input just press ENTER in response to prompt- 
2230 IN-0 
2240 PRINT 

2250 X*«" 1 PRINT*Enter XrON-UlTJONPUT X* 

2260 IF X*=“ GOTO 2310 

2270 T»=LEFT»(X*,1> 5 IF TM*0* OR T*»T" OR T*«-2* OR T$=*3* OR T*=*4* OR T*=*5* OR TS=*6* ORT4-T* OR 
T»«‘Q* OR T*=-9- OR TS=-.- OR T*=“-- OR T*=-*- THEN GOTO 2290 
2280 PRINT"Reenter last value* 1 GOTO 2250 
2290 IN-IN+t ! X(IN)=VAL(X*) 1 YttNl=X(IN) 

2295 IF IN=1 THEN 2250 ELSE IF XONKXUN-1) THEN 15=0 
2298 IF IK>199 THEN 2310 
2300 GOTO 2250 

2310 PRINTS IF IN-0 THEN PRINT-NO DATA ENTERED*; ELSE IF IN-1 PRINT*ANALYSIS CAN NOT BE PERFORM 

ED ON A SAMPLE OF 1*} ELSE PHINT-Number of data points -"OH : FOR 1=1 TO 500 ! NEXT 

2320 FOR l-IN+l TO 200 S X(I)«0 S Yai=X(I>! NEXT 

2330 IF IN=0 THEN RETURN ELSE IF IN-1 THEN IN«0 : RETURN 

2340 IF 15=1 GOTO 3800 
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SHELL SORT 


2400 REM * * * SHELL SORT ROUTINE **• 
2410 CLS ! PRINT Q438,-*** ORDERING DATA 
2420 IF IN>3 GOTO 2510 
2430 IF IN-3 GOTO 2460 

2440 IF X(1»X<2> THEN XT-Xtll : XU1-XC2) ! X(2)-XT 
2430 GOTO 2300 

2460 IF X(i»X(2) THEN XT-X12) ! X<2)*X(1> ! X(1)«XT 
2470 IF X(2»X(3) THEN XT-X13) ! X<3)-X(2) ! X(2)-XT 
2400 IF X(1»X(2) THEN XT=X<2> ! X12)-X<1) S X(1)-XT 
2300 FOR I«1 TO IN : Y(I)-X(I> 1 NEXT I ! GOTO 2660 
2310 IM-INT(L0G(IN)/L0G<2)>-1 
2520 H ■ 2CIM - 1 
2330 K ■ H 

2340 FOR I » K ♦ 1 TO IN 
2330 XT*X<I) 

2360 J - I - K 
2370 IF J < 1 GOTO 2620 
2300 IF XT >- X(J> GOTO 2620 
2390 XU*K) » X(J) 1 Y(J>K) * Y1J) 

2600 J ■ J - K 

2610 GOTO 2370 

2620 XU»K) - XT : YU«-K) - JCT 

2630 NEXT I 

2640 H ■ (H - l)/2 

2630 IF H >■ 1 GOTO 2330 

2660 GOSUB 3800 

2670 RETURN 


INPUT DATA VIA TAPE 

2700 REM INPUT DATA VIA TAPE **• 

2710 CLS i 14-0 ! PRIHT’Insert DATA TAPE into ESF-80* 

2720 PRINT 1 PRINT'Pres* SPACE BAR after inserting the tape- 

2730 IF INKEY* O ■ * GOTO 2730 

2740 PRINT : INPUT "Enter Data File •‘JFI 

2730 GOPEN FI 

27601-1 

2770 0INFUT Xa> 

2780 Ya>-xa) 

2783 IF 1-1 THEN 28p0 ELSE IF XdKXU-l) THEN 13-0 
2790 IF XID-999999 GOTO 2810 
2800 I-I-M ! GOTO 2770 
2810 0CLOSE 

2820 IN*1-1 ! PRINT 5 PRINT-NUMBER OF DATA POINTS -*3-1 1 FOR 1-1 TO 1300 S NEXT I 1 IF 13-0 THEN GOT 
O 2410ELSE GOTO 3800 


mm SVttMART STATISTICS flBLX 

2900 REM *** INPUT SUMMARY STATISTICS ONLY *** 

2910 14-1 : CLS ; PRINT TAB(20)“SUMMARY STATISTICS INPUT- 

2920 PRINT ! INPUT-Simole Siie“;iN 1 IF IN-1 OR IN-0 GOTO 2310 

2930 PRINT (?192,“Is sample assumed to be from an Exponential Distribution?" 

2940 KP-256 S L»*-y" 1 R*»*n* 1 GOSUB 600 

2930 IF IQ-"n-OR IQ--N-GOTO 2990ELSE PRINT 0192,-WUl you enter-5LBS{-M-nffl*I*ean or Exponential Para 
meter -JLBtrL'JRBSJ-ambda-; 

2960 KP-256 1 L»-*m* 1 R*-“l- ! GOSUB 600 

2970 IF ICK-m- OR IQ-’M* THEN PRINT 1 INPUT'Enter Mean'lXB ! V1-XB*XB J V-Vl 5 GOTO 3170 
2980 PRINT 1 INPUT'Enter Lambda-JXB ! XB-l/XB ! V1-XB«XB ! V-Vl 1 GOTO 3170 
2990 PRINT 3192, *Is the sample assumed to be from a Bernoulli Population?" 

3000 KP-256 ! L*-"y- ! RS-'n" ! GOSUB 600 
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3010 IF ZQ=“n" OR Z0»"N" THEN GOTO 3030ELSE PRINT ! PRINT'Enter number of occurrences out of“;iN!*trials" 
J1INPUT XB : IF XB>IN OR XBCO GOTO 3020ELSE IF XBOINT(XB) GOTO 3020ELSE XB-XB/IN ! V1 = XB«(1-XB)! 
V*Vl 1 GOTO 3170 

3020 PRINT G192.CHR*<250) ! PRINT Q192,*YOU MOST ENTER AN INTEGER FROM 0 TO*?IN 5 PRINT: FOR 1*1 TO 
2000 : NEXT : GOTO 3010 

3030 PRINT 0192,BL* ! PRINT 1 PRINT 0192.**:SINFUT"Sample Mean'IXB 
3040 PRINT I PRINT"WiU you enter variance or std. deviation ?* 

3050 KP=339 : L*»*v* ! Rs*"d" ! GOSUB 600 

3060 IF ZQ=*v" OR ZQ*“V" THEN V**"Sample Variance* ELSE VS**Sample Std. Deviation* 

3070 PRINT ! PRINT V*!!INPUT SV 

3080 IF V**“Sample Variance* THEN GOTO 3100ELSE V1*SV«SV ! V*((IN-11/IN»*V1 
3090 GOTO 3170 

3100 PRINT"Is the sample variance! *5CHRf(123i:*l"lCHR4(125);* the M.L.E. (i.e. i/n) 

3110 PRINT* *iCHR«(123>;*2*:CHR»(125n* unbiased estimator (i.e. 1 'n-l> 

3120 PRINT 0640,"SELECT APPROPRIATE ESTIMATOR ";CHR*<123>I* "iCHRV125> 

3130 ZI=INKEY* ! IF ZI=“* GOTO 3130ELSE PRINT 067O.ZI!FOR 1*1 TO 100! NEXT 
3140 IF ZIO*l* AND ZIO"2“ THEN PRINT 0670," *I!GOTO 3130 
3150 IF ZI=*1* V*SV ! V1*<IN/<IN-1»*SV 
3160 IF ZI="2* THEN Vt=SV ! V*«IN-11/IN)»SV 

3170 CLS ! PRINT"When summary statistics are entered in lieu of data, several subroutines are not available. 
These include!" 

3180 Z1*CHR*(123):Z2»CHR*(1251!PRINT!PRINT* *!Z1!" Data Transformations *!Z2 ! PRINT" "iZli* Five Numb 

er Summary *JZ2 ! PRINT" *;Z1I* All Data Plots *!Z2 I PRINT" "5Z1I* All of the Subroutines in Prog. #4 *;Z2 
3190 FOR 1*1 TO 200 ! X(I)*0 ! Y(I1*0 ! NEXT 
3200 PRINT 0396,"PRESS SPACE BAR TO CONTINUE* 

3210 IF INKEYSO* * GOTO 3210 
3220 RETURN 


SME PAT A IS-IAFS 

3300 REM *•* SAVE DATA TO TAPE •** 

3310 CLS ! PRINT“INSTRUCTIQNS:* 

3320 PRINT" (1) Insert DATA TAPE 

3330 PRINT* (2) If you use a blank tape then FILE » will be 1 

3340 PRINT" (3) If the tape is not blank then the FILE • must be 

3350 PRINT* the number of the last data file recorded PLUS 1 

3360 PRINT ! PRINT "When the tape is ready press SPACE BAR 

3370 IF INKEY* O " * GOTO 3370 

3380 PRINT ! INPUT“Enter FILE »"iIO 

3390 YGN+1)*999999 

3400 GOFEN IQ 

3410 FOR 1*1 TO IN+1 

3420 SPRINT Y(I) 

3430 NEXT I 
3440 0CLOSE 

3450 PRINT IN’,*data points recorded as DATA FILE **!IO ! FOR 1-1 TO 2000 ! NEXT ! RETURN 


DATA TRANSFORMATIONS 

3500 REM ♦** DATA TRANSFORMATION SUBROUTINE 
3510 CLS ! PRINT TAB(22)"DATA TRANSFORMATIONS* 

3520 PRINT ! PRINT "MENU!" ! PRINT*<1> NO Transformation 
3530 PRINT"(2> Square Root Transformation 
3540 PRINT"(2) Cube Root Transformation 
3550 PRINT"(4) Logarithmic Transformation 

3560 PRINT ! PRINT* NOTE II Other transformations must be made before data entry 
3570 PRINT" NOTE 21 Raw data may be recovered by returning to this sub- 
3580 PRINT* routine and selecting Option l 

3590 PRINT" NOTE 31 Option 2 will not work with negative data and option 4 requires all data great 

er than zero 
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3400 KP-896 : KS-4 : GOSUB 500 

3610 CLS ! PRINT*E X EC UTING TRANSFORMATION OPTION*!IO 
3620 FOR 1=1 TO IN 

3630 ON 10 GOTO 3640, 3650, 3670, 3690 
3640 Y(I)-XU) ! GOTO 3710 

3630 IF XU) < 0 CLS ! PRINSCAN NOT TAKE THE SQUARE ROOT OF NEGATIVE DATA!* 5 FOR J-l TO 2000 1 KE 

XT J : GOTO 3300 

3660 Y<I>-S0R<X(I)) i GOTO 3710 

3670 YU)-ABS(XU))C.3333333 I IF XU) < 0 YU)-Ya)*-t 

3680 GOTO 3710 

3690 IF XU) <* 0 CLS ! PRINT’CAN NOT TAKE THE LOG OF DATA LESS THAN OR EQUAL TO 2ERO* t FOR J-l T 
O 2000 ! NEXT J ! GOTO 3500 
3700 TU)-LOG(XU» 

3710 NEXT I 

CALCULATE SAMPLE STATISTICS 

3800 REM • • • CALCULATE SAMPLE STATISTICS • *• 

3810 CLS : PRINT 0454,*» INITIALIZING DATA BASE** 

3820 SM-0 ! SUM-0 ! 52-0 : FOR 1-1 TO IN ! SM-SM+YU) 1 DV-YU1-5UM ! SUM-SUM+DV/11 S2-S2+DV • fYUHSUF 
> t NEXT 

3830 XB-SM/IN ! V1-S2/UN-1) ! V-S2/IN 

3840 Gl-0 ! G2-0 i FOR 1-1 TO IN ! Gl*Gl+(Ya>-XB>l3 5 G2-G2*(Ya>-XB>£4 ! NEXT I ! Gl-Gl/IN/Vtl.5 I G2-1G2 

/IN/Vt2>-3 

3830 RETURN 

RAMPQM NVMB5B QgNgHATOT 

4000 REM RANDOM NUMBER GENERATOR SUBROUTINE 

4010 CLS : 14-0 : print0i8,-random number generator* 

4020 PRINT ! PRINT ! PRINT'Distribution Options!* 

4030 PRINT* (1) Normal 
4040 PRINT* (2) Uniform 
4050 PRINT* (3) Exponential 
4060 PRINT* (4) Bernoulli 
4070 KP-643 ! KS*4 ! GOSUB 300 

4080 PRINTQ704,"*} i INPUT'How many random numbers do you want" JIN 

4083 IF IN<-1 THEN PRINT 0640, CHRS(250> ! PRINT 0640,'WHAT? ... Please request a sample size of at least 2 
* 1 GOTO 4080 

4090 ON IO GOSUB 4300, 4400, 4600, 4700 

4100 CLS • PRINT'Do you want to see the random tangle *!! L*-*y* ! R*-*n* ! XP-33 ! GOSUB 600 
4110 IF ZQ-*n* OR ZQ-*N* THEN 2410ELSE PRINT 
4120 J-0 

4130 FOR 1-1 TO IN 

4140 J-J+l ! IF J > 3 PRINT* * 

4130 IF J > 3 THEN J»1 

4160 PRINT USING*#*#e*e*.###“5XU>!:PRINT* *! 

4170 IF INTO/30W/50 GOTO 4190 
4180 GOTO 4210 

4190 PRINT ! PRINT*PRESS SPACE BAR TO DISPLAY NEXT SET OF NUMBERS* 

4200 IF INKEY* O * " GOTO 4200 
4210 NEXT I 

4220 PRINT ! PRINT*Press SPACE BAR to continue* 

4230 IF IN KEY* O * * GOTO 4230 
4240 GOTO 2410 

NORMAL RNG 

4300 REM ••• NORMAL RNG ••* 

4310 CLS ! PRINT 010,’NORMAL RANDOM NUMBERS BY BOX-MULLER METHOD- 

4320 PRINT ! PRINT 

4330 INPUT’What is the mean“!MU 

4340 INPUT'What is the vanance*?VAR ! SIGMA-SQRCVAR) 

4330 FOR 1-1 TO IN STEP 2 
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4360 Ul=RND<0) 5 U2=RND<0> 

4370 X<I>*SOR(-2»LOGIU1»*COS<6.28315»U2>:X(I>»X(I>«SIGMA ♦ MU 5 Y(I)-X(I) 

4380 X<I+l>*SQR<-2*LOG(Ul>J*SIN<6.283l5*U2>:Xa-H>-X(Irl>»SIGMA ♦ MU S Y(Itl)*X<H-l) 
4390 NEXT I 
4395 RETURN 


HMESR M RN -Q 

4400 REM ••• UNIFORM RNG ••• 

4410 CLS 

4420 PRINT02O,“UNIFORM RANDOM NUMBERS- 

4430 PRINT ! PRINT ! PRINT- (1) Continuous 

4440 PRINT* (2) Uniform Integers 

4450 KP*384 ! KS=2 S GOSUB 500 

4460 PRINT : INPUT“Enter Lower LimifJLL 

4470 INPUT’Enter Upper Linut-JUL 

4480 ON 10 GOTO 4490, 4540 

4490 FOR I*! TO IN 

4500 X1I1»RND(0)*<UL-LL) ♦ LL 

4510 Y(I)»X<I> 

4520 NEXT I 
4530 RETURN 
4540 FOR 1*1 TO IN 

4550 X(I)*INT(RND(0>e<UL-LL) + LL ) ! IF Xtl) O LL THEN X(I)«X<I)*1 
4560 Ya>*xa» 

4570 NEXT I 
4580 RETURN 


EXPONENTIAL RNG 

4600 REM ••• EXPONENTIAL RNG 

4605 CLS ! PRINT 015,"EXPONENTIAL RANDOM NUMBERS- 

4610 PRINT 0128,-Will you enter the Exponential Parameter Lambda or the Mean* : KP* 192 : Lf**]* 1 R<-*m* 1 GO 
SUB 600 

4620 IF 7Q--1- OR ZO*L* THEN ZP»‘Lambda* ELSE ZP*“Mean* 

4630 PRINT 0256,-Enter -JZPJIINPUT MU 
4640 IF ZP»*Lambda- THEN MU-l/MU 
4650 FOR >1 TO IN 
4660 U1«RND(0> 

4670 Xa>» -MU ♦ LOG(Ul) 

4680 rai-xa) 

4690 NEXT I 
4695 RETURN 


BBUWVLU RNQ 

4700 REM ••• BERNOULLI RNG eee 
4710 CLS 

4720 PRINT015,-BERNOULLI RANDOM NUMBERS 

4730 PRINT 1 INPUT'What is the probability of success on a single trial"IPS 
4740 FOR I-t TO IN 
4750 X(I)“RND(0> 

4760 IF Xdl <» PS THEN XdM ELSE XlD-0 
4770 NEXT I 
4780 RETURN 
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PROBABILITY DISTRIBUTIONS MENU 


5000 

5010 

5020 

5030 

5040 

5050 

5060 

5070 

5080 

5090 

5100 

5110 

5120 

5130 

5140 

5150 

5160 

5170 

5180 

5190 


REM *•♦ PROBABILITY D I 5 T R I B U T I O M S SUBROUTINE 
CLS ! PRINT TAB(18)“PROBABILITY DISTRIBUTIONS 
PRINTIPRINT! 12*0 
PRINT-MENU!" 

PRINT“tl) Normal Distribution 
FRINT"(2) Binomial Distribution 
PRINT-13) Negative Binomial Distribution 
FRINT"t4> Poisson Distribution 
PRINT"(5) Exponential Distribution 
FRINT"(6) Chi-Square Distribution 
PRINT"(7) Student's - t Distribution 
PRINT"(8) F-Oistnbution 
KP*896 ! KS*8 ! GOSUB 500 

ON IO GOSUB 6000,7000,8000,9000,10000,11000.12000,13000 
IF 10*2 OR 10*3 OR 10*4 OR 10=5 THEN RETURN 

PRINT ! PRINT ! PRINT "Prob. 1 X <*-;PX;*) * *! ! PRINT USING**.*###*!F0X 
PRINT "Prob. < X > “IPX:*) * “J! PRINT USING-#.i*#*“iQOX 
PRINT 9960," Run Again or Quit?*!! KP*978 i L$*-r" ! Rs-*q- ! GOSUB 600 
CLS ! 12*0 ! IF Z0*“r“ OR ZQ-*R* GOTO 5130ELSE RETURN 
GOTO 1020 


inverse aasa giiaMS 

5200 REM ••• INVERSE DISTRIBUTIONS * * * 

5210 CLS ! PRINT TAB(15)“INVERSE PROBABILITY DISTRIBUTIONS 

5220 PRINT 9192,"MENU!-!FRINT*(1) Inverse Normal 

5230 PHINT"(2) Inverse Chi-Square 

5240 PRINT-O) Inverse Student's - t 

5250 FRINT"<4) Inverse F 

5260 PRINT"(5) Inverse Exponential- 

5270 K P-960 ! KS-5 ! GOSUB 500 

5280 12*0 

5290 ON IO GOSUB 14000, 15000, 16000, 17000, 18000 
5300 IF 10*5 GOTO 5350 

5310 PRINT ! PRINT ‘For F(x) » Prob. ( X O x ) «-{FR 

5320 IF XQ-9E36 THEN PRINT-x - infinity- ELSE IF XQ*9E-36 PRINT *x ■ - infinity- 
5330 IF XQ-VE36 OR XQ-9E-36 THEN GOTO 5350 

5340 PRINT -x - -J! IF ABSCXQD1000 THEN PRINT INT(XQ) ELSE IF ABS(XQ»100 PRINT USING-**#*.»-;XQ EL 
SE IF ABS1XQD10 PRINT USING" •##.*# a 5XQ ELSE PRINT USING"«*.**»*;XQ 
5350 PHINT 9960,* Run Again or Quit"!! LS*“r* i RS*-q- 1 KF-976 ! GOSUB 600 
5360 IF ZQ-'R* OR ZQ-'r* GOTO 5280ELSE RETURN 


NORMAL DISTRIBUTION 

6000 REM ••• NORMAL DISTRIBUTION SUBROUTINE 
6010 IF 12-1 GOTO 6070 

6020 CLS {PRINT 920,‘NORMAL DISTRIBUTION* 

6030 PRINT!PRINT!INPUT-Enter Mean'iMU 
6040 INPUT-Enter Vanance“{VAR 

6050 IF VAR<*0 PRINT*** ERROR ** "iVAR'IS NOT A VALID VARIANCE-!GOTO 6040 
6060 INPUT-Enter Normal Statistic x *{NX {IF 12-0 THEN PX-NX 
6070 SD«SQRtVAR) 

6080 XN*(NX-MU)/SD!AX>ABS(XN) 

6090 XZ-4989423*EXP((-AX*AX)/2) 

6100 NT*1/<1*.2316419*AX) 

6110 0OX*XZ*«.3193813*NTM.3363638*NT£2>*(1.781478*NTC3>-<l.82t256*NTt4»U.330274*NTC3H 

6120 IF XNCO GOTO 6150 

6130 FOX-I-QOX 

6140 RETURN 

6150 FOX-QOX 
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4160 QOX- 1-FOX 
6170 RETURN 


BIHQMIAfaPCT OTV .TOM 

7000 REM * * « BINOMIAL DISTRIBUTION SUBROUTINE * • • 

7010 IF 12-1 GOTO 7330 

7020 CLS S PRINT TAB(20)"BINOMIAL DISTRIBUTION 
7030 PRINT 

7040 PRINT" The normal scenario for the Binomial Distribution is K" 

7050 PRIHT"oaurrences out of N trials where P is the probability " 

7040 PRINT'of occurrence on a single trial. 

7070 PRINT ! INPUT'Enter N"!N 

7080 IF N<0 OR NOINT(N) THEN PRINT"** ERROR ** N MUST BE A POSITIVE INTEGER"tGOTO 7070 
7090 INPUT"Enter P"iPB 

7100 IF PB-l THEN PRINTTHE PROBABILITY OF*!NJ"OUT OF*iN‘," a 1.0“SGOTO 5170 

7110 IF PB-0 THEN PRINTTHE PROBABILITY OF 0 OUT 0F"SN}"- 1.0"!GOTO 5170 

7120 IF PB<0 OR PB>1 THEN PRINT"** ERROR «« "SPBSTS NOT A VALID PROBABIUTY"!GOTO 70*0 

7130 INPUT'Enter K"iXX 

7140 K-XK 

7150 IF K<0 OR KOXK THEN PRINT"** ERROR ** K MUST BE A POSITIVE INTEGER" iGOTO 7130 

7140 IF K>N THEN PRINT"THE PROBABILITY OF “IKS" OUT OF "INI" ■ O'SGOTO 5170 

7170 IF K-0 OR N-K THEN BC-1 i GOTO 7240 

7180 IF K<0 THEN PR-0 i GOTO 7250 

7190 TB-N 1 BB-K 

7200 IF N<120 GOTO 7230 

7210 GOSUB 7740 

7220 GOTO 7250 

7230 GOSUB 7430 

7240 PR - BC • PBCK * (1-PBHtN-K) 

7250 PRINT S PRINT"Prob. of exactiyiKl'occurrencBS out oF*;Ni"trials « "{SPRINT USING"#.####‘iPR 
7240 PRINT 096O,*Do you want Cumulative Probability "ll KP-992 S L*-"y" J RM'n* S GOSUB 400 
7270 IF ZQ-"n* OR ZQ«“N* PRINT Q940,BL*{1 GOTO 5170 
7280 PRINT S INPUT"Enter K'iXK 

7290 K-XK S IF K>N THEN GOTO 7160ELSE IF K<0 OR KOXK PRINT** ERROR •* K MUST BE A POSITIVE IN 
TEGER'SGOTO 7280 
7300 IF N<-45 GOTO 7330 

7310 IFPB3.1 AND FB<.9 AND PB*N>5 THEN PRINT"NORMAL APPROXIMATION HILL BE USED" SGOSUB 75801 
GOTO 7410 

7320 IF PBC-.l OR FB3-.9 THEN PRINT’PQISSON AFROXIKATION WILL BE USED" ! GOSUB 7640t GOTO 7410 
7330 PR-0 

7340 IF K<0 GOTO 7410 

7350 IF K >» N THEN PR-1 1 GOTO 7410 

7340 FOR BB-0 TO K 

7370 TB-N 1 GOSUB 7430 

7380 PR-PR* BC • PBCBB • (l-PB)t(N-BB) 

7390 IF PR>-1 GOTO 7410 
7400 NEXT BB 

7410 IF 12-1 RETURN ELSE PRINT 1 PRINT’Prob. of'SKPor Ins occurrences out of*!Ni"trials > "ISPRINT USING" 

•.•*«*tPR 

7420 PRINT l GOTO 7240 

7430 BC-0 I B«0 S BT-0 1 SB-BB I QB-TB-BB 

7440 FOR BI-1 TO TB 

7450 IF SB>-QB GOTO 7510 

7440 IF SB >• BI THEN B-B ♦ LOG(BI) 

7470 IF OB > BI GOTO 7550 

7480 IF OB - BI THEN BC—B 5 GOTO 7550 • 

7490 BC-BC * LQGlBI) 

7500 GOTO 7550 

7510 IF QB >- BI THEM B-B ♦ LOGCBI) 

7520 IF SB > BI THEN GOTO 7550 
7530 IF SB - BI THEN BC—B 1 GOTO 7550 
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7540 BC=BC ♦ LQG(BI) 

7550 NEXT BI 
7540 BC-EXP1BC) 

7570 RETURN 
7530 MU=PB 

7590 VAR-1PB • <1-PB> )/N 

7600 NX-UC ♦ »5)/N 

7610 12-1 

7620 GOSUB 6000 

7630 12*0 

7640 PR-FOX 

7650 RETURN 

7660 IF PBC-.l THEN LAM-N*PB ELSE LAM»N*U-PB1 
7470 TI-l 

7400 IF PB<«.1 THEN NC-K ELSE NC-N-K-l 

7690 12-1 

7700 GOSUB 9230 

7710 12-0 

7720 IF PB>=.9 THEN PR-l-PR 
7730 RETURN 

7740 REM * « * NORMAL APPROXIMATION 
7750 IF PB>*.9 AND (1-PB)*N<5 THEN GOTO 7080 
7760 IF PB<.1 AND PB«N<5 THEN GOTO 7880 
7770 MU-PB 

7780 VAR-(PB * (1-PBH/N 
7790 NX-<K*.5)/N 
7800 12-1 
7810 GOSUB 6000 
7820 PR-FOX 
7830 NX-0C-.5J/N 
7840 GOSUB 6000 
7850 PR-FR-FOX 
7860 12-0 
7870 RETURN 

7880 REM •** POISSON APPROXIMATION 

7890 IF PBC-.l THEM LAM*N*PB ELSE LAM-N*<1-PB> 

7900 TI-l 

7910 IF PB<-.1 THEN NC-K ELSE NC-N-K 

7920 12-1 

7930 GOSUB 9000 

7940 12-0 

7950 IF PR}-.9 THEN PR-l-PR 
7960 RETURN 


HSQAim BINPWftt, Pi?TRfflVTtgN 

8000 REM * • * NEGATIVE BINOMIAL SUBROUTINE **♦ 

0010 IF 12-1 GOTO 8170 

8020 CLS 1 PRINT TABU5PNEGATIVE BINOMIAL DISTRIBUTION 

8030 PRINT 1 PRINT* The normal scenario for the Negative Binomial Distribution is* 

8040 PRINT'probability of exactly K failures prior ta the Rth success* 

8050 PRINT "when the probability of success on a single trial is P 
8060 PRINT 5 INPUT'Enter K*?XK ! K-XK 

8070 IF K<0 OR KOXK THEN PRINT*** ERROR •* K MUST BE A POSITIVE INTEGER*'GOTO 8060 
8080 INPUT*Enter R-JRA 

8090 IF RA<0 OR RAOINT(RA) THEN PRINT*** ERROR •« R MUST BE A POSITIVE INTEGER*!GOTO 8080 
8100 INPUT"Enter P*;PB 

8110 IF PB-0 THEN PRINT'THE PROBABILITY OF*;K’,*FAILURES PRIOR TO*iRA!*SUCCESSES - 1.0'IGOTO 5170 

8120 IF PB-1 THEN PRINTTHE PROBABILITY OF*JKi*FAILURES PRIOR TO*!RA:*SUCCESSES * O.OVGOTO 5170 

8130 IF PB<0 OR FB> 1 THEN PRINT*** ERROR ** 'iPBJ'IS NOT A VALID PROBABIUTY"U30TO 8100 

8140 IF K-0 AND RA-0 THEN PRINT 1 PRINT**** Reenter parametersUGOTO 8060 

8150 IF K-0 THEN PR - PBCRA 1 GOTO 8200 

8160 IF RA-0 THEN PR - (l-PB)CK ! GOTO 8200 
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6170 TB-RA+K-1 J BB=K 

8160 GOSUB 7430 

8190 PR-BC * PBCRA ♦ (l-PB)CK 

6200 PRINT I PRINT"Prod, of ezactly'IKI’failures prior to'IRA;" successes" 

8210 PRINT’is equal to ‘SPRINT USING"».####"JPR 

8220 PRINT 0960,‘Do you want Cumulative Probability "J*. KP-992 i LS="y“ 5 R*»*n" 1 GOSUB 600 
8230 IF ZQ="N“ OR ZQ='n" THEN PRINT 0960, BL*J: GOTO 5170 
8240 PR-0 ! PRINT ! INPUT’Enter IC'JXK 

8250 K-XK 5 IF K<0 OP. KOXK THEN PRINT"** ERROR *« tC MUST BE A POSITIVE INTEGER-.GOTO 8240 
8260 IF K=0 AND RA-0 THEN PRINT : PRINT"*** Reenter parameter’lGOTO 8240 
8270 IF K-0 THEN PR-PBCRA 5 GOTO 8330 

8280 IF RA-0 THEN FOR 1=1 TO K ! PR=PR+(l-PB)tI ! NEXT ! GOTO 8330 

8290 FOR BB-0 TO K 

8300 TB-RA * BB - 1 1 GOSUB 7430 

8310 PH-PR ♦ (BC * PBCRA * (l-PB)CBB) 

8320 NEXT BB 

8330 IF RA«t IX="st* ! GOTO 8360 
8340 IF RA=2 IX-"nd" t GOTO 8360 
8350 IF RA-3 ZX-"rd* ELSE ZX="th" 

8360 PRINT ! PRINT "Prob. of"!K;"or fewer failures prior to the’JRAJ IX!" lucres*" 

8370 PRINT'is equal to "J1PRINT USING"*.*###"JPR 
8380 PRINT! GOTO 8220 


POISSON DISTRIBUTION , 

9000 REM «•• POISSON DISTRIBUTION SUBROUTINE 
9010 IF 12-1 GOTO 9090 

9020 CLS 5 PRINT TAB(20)"POISSON DISTRIBUTION 
9030 PRINT ! INFUT’Enter Rate <lambda)*JLAM 

9040 IF LAM<*0 PRINT"** ERROR •* ENTER A VALUE > O-.GOTO 9030 
9050 INPUT"Enter time period tt)"5TI 

9060 IF TK-0 PRINT"** ERROR •* ENTER A VALUE > O'tGOTO 9050 
9070 INPUT "Enter number of occurrences" iNC 

9080 IF NCCO OR NCOINTtNC) PRINT"** ERROR ♦* ENTER A POSITIVE INTEGER VALUE"tGOTO 9070 
9090 PD-0 

9100 IF NCOI GOTO 9140 
9110 FOR PI-2 TO NC 
9120 PD-PD ♦ LOG (FI) 

9130 NEXT PI 

9140 PR“(-LAM*TI) ♦ (NC»LOG(LAM*TD) - PD 
9150 PR-EXP(PR) 

9160 IF 12-1 RETURN 

9170 PRINT 1 PRINT "Prob. of exactly" iNCJ-occutrencas in" 015" time units ■ "JIPRINT USING"9.#*#*"JFR 
9180 PRINT 0960,"Do you want Cumulative Probability? "J! KP-992 1 Ls-*y“ ! RS-"n" S GOSUB 600 
9190 IF 2Q=*N* OR ZQ=“n" THEN PRINT 0960, BLSJJ GOTO 5170 
9200 PRINT 1 INFUT’Enter number of ocxurrences'JNC 

9210 IF NCCO OR NCOINTtNC) THEN PRINT"** ERROR •• ENTER A POSITIVE INTEGER VALUE"IGOTO 9200 

9220 IF NC-0 THEN FR*EXP<-LAM«TI>:GOTO 9330 

9230 PR-0 

9240 FOR I-l TO NC 

9250 PU-0 ! PD-0 

9260 FOR PI-1 TO I 

9270 FD-PD ♦ LOG (FI) 

9280 NEXT PI 

9290 PU»(-LAM*TI) * d*LOGlLAM*TI» - FD 
9300 FR-PR + EXPtPU) 

9310 NEXT I 

9320 PR-PR ♦ EXP(-LAMeTI) 

9330 IF 12-1 THEN RETURN 

9340 PRINT.; PRINT"Prob. of"JNCJ"or fewer ocrurrenaes" 

9350 PRINT"ln"rnj"time units ■ "JJPRINT USING "#.*#***JPR 
9360 PRINT : GOTO 9180 
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EXPONENTIAL DISTRIBUTION 


10000 REM EXPONENTIAL DISTRIBUTION SUBROUTINE 

10009 CLS PRINT TAB(20>*EXPONENTIAL DISTRIBUTION* _ 

10010 PRINT : PRINT* An exponential distribution is often used in a practical* i PRINT *proble«n to represent 
the distribution of the tune that elapses* 

10019 PHINT“before the occurrence of some event. Perhaps the moat common* 

10020 PRINT*example is the 'interarrival times * between renewals in a* 

10030 PRINT'Poisson Process. A particularly nice feature of the exponen-* 

10040 PRINT’tial distribution is its memoryless property. That is. it is* 

10090 PRINT*not necessary to consider past occurrences of an event in order* 

10060 PRINT’to calculate probabilities far future occurrences of the event.* 

10070 PRINT ! PRINT * Press the SPACE BAR to proceed.*; 

10080 IF INKEYS O * * GOTO 10080 
10090 CLS 

10100 PRINT*NOTE* While the parameters and answers are expressed in terms* 

10110 PRINT*af time, be aware that this is not the only context possible.* 

10120 PRINT ! PRINT'You may use either the Exponential Rate. LAMBDA (1/MEAN) or the MEAN* 

10130 KP-384 J L*-*l* ! R*-“m* ! GOSUB 600 
10140 ZEF-ZQ 

10190 PRINT S PRINT ! PRINT*Enter MIF ZE»*1* OR ZE-*L* INPUT*LAMBDA“!LAM ELSE INPUT"MEAN*[LAM 
10160 IF LAM<=0 THEN PRINT*** ERROR ** YOU MUST ENTER A VALUE > 0*!GOTO 10190 
10170 IF ZE«*M* OR ZE»“m* THEN LAM-1/LAM 
10180 INPUT’Enter number of time umts’STI 

10190 PRINT ! PRINT’Prob. of“JTIi*time units or less until the next* S PRINT"occurrence Is «*!(PRINT USING**.# 
•♦•*;i-EXP(-LAMBDA*TI) - 

10200 PRINT I PRINT’Prob. of more than*!TC5*time units until the next* ! PRINT’occurrence is -‘{SPRINT USING* 
#.##***5EXP<-LAMBDA*TI) 

10210 PRINT ! PRINT’Expected time until the next occurrence is*i 1 /LAMBD ANPRINT* time units* 

10220 PRINT 0960, BLS : GOTO 5170 


CHI - SQUARE DISTRIBUTION 

11000 REM ♦** CHI-SQUARE DISTRIBUTION SUBROUTINE 
11010 IF 12-1 GOTO 11080 

11020 CLS! PRINT 320,"CHI-SQUARE DISTRIBUTION* 

11030 PRINT 0192,**i:iNPUT*Enter Chi-Square Statistic x *JCX 

11040 IF CX<«0 PRINT*** ERROR »• CHI-SQUARE STATISTIC MUST BE > 0*SGOTO 11030 

11050 IF 12-0 THEN PX-CX 

11060 INPUT'Enter Degrees of Freedom* IDF 

11070 IF DFC-0 OR DFOINT(DF) THEN PRINT-** ERROR ** ENTER AN INTEGER > 0*K3OTO 11060 
11080 IF DF>30 GOTO 11290 
11090 DP-1 

11100 FOR I»INT(DF> TO 2 STEP -2 
11110 DP-DP*I 
11120 NEXT I 

11130 NF-C XI <INT((DF+1 )/2))*EXP(-C X/21/DP 
11140 IF INT(DF/2)-DF/2 GOTO 11170 
11150 PC-SQR12/CX/3.141593) 

11160 GOTO 11180 

11170 PC-1 

11180 LC-t 

11190 MC-1 

11200 D-DF 

11210 D-D+2 

11220 MC-MC*CX/D 

11230 IF MCC.000001 GOTO 11260 

11240 LC-LC*MC 

11290 GOTO 11210 

11260 FOXaPC » NP* LC 

11270 QOX-l-FOX 

11280 RETURN 
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11290 NX=«CX/DF>t(l/3> - (l-(2/(9*OFmi/SQR(2/(9«DF» 
11300 MU=0i VAR=1 > 12=1 
11310GO5OB 6000 
11320 RETURN 


STUDENT'S -1 DISTRIBUTION 

12000 REM * * * t-DISTRIBUTION SUBROUTINE 

12010 IF 12*1 GOTO 12060 

12020 CLSiPRINT 020,“T-DISTRIBUTION* 

12030 PRINT ! PRINT 1 IKPUT’Enter t-Statistic x *iTX 
12040 INPUT'Enter Degrees of Freedom"iDF 

12030 IF DF<=0 OR DFOINT1DF) THEN PRINT"** ERROR ♦* ENTER AN INTEGER > 0‘tGOTO 12040 

12060 IF TXCO 13=1 ELSE 13=0 

12070 IF DF>25 GOTO 12500 

12080 THETA*ATN(ABS(TX)/SQH(DFl> 

12090 IF DFOl GOTO 12140 
12100 AT=.63661<’8»THETA 
12110 FOX-(AT+l)/2 
12120 QOX-l-FOX 
12130 GOTO 12450 

12140 IF INT(DF/2)*DF/2 GOTO 12270 

12130 IF DF*3 GOTO 12420 

12160 TS=0 1 TP*1 S TC=COSITUETA) 

12170 FOR 1*1 TO INT«DF-3>/2> 

12180 TP»TP*2#I/(1*2*I> 

12190 TC*TC*<COS(THETA)C2) 

12200 TS*TS+(TP*TC) 

12210 NEXT I 

12220 TS-TS+C OS (THETA) 

12230 AT».6366198*(THETA*(SIN(THETA>*TSH 

12240 FOX«(AT>l)/2 

12230 QOX-l-FOX 

12260 GOTO 12450 

12270 IF DF02 GOTO 12320 

12280 AT-SIN(THETM) 

12290 FOX»(AT+l)/2 
12300 QOX-l-FOX 
12310 GOTO 12430 
12320 TS*0 ! TP-1 !TC*1 
12330 FOR I«1 TO INTHDF-2J/2) 

12340 TF-TF*«2*I>-1>/(2*I) 

12330 TC-TC*(CQS(THETA)C2> 

12360 TS»TS+(TF*TC> 

12370 NEXT I 

12380 AT-SIN (THETA HH1+TS) 

12390 FOX*(AT+l>/2 
12400 QOX-l-FOX 
12410 GOTO 12450 

12420 AT- >6366198*(THETA*(S1N(THETA>*COS (THETA))) 

12430 FOX-(AT*l)/2 
12440 QOX-l-FOX 
12430 PX-TX 
12460 IF 13-0 RETURN 
12470 QOX-FOX 
12480 FOX-l-QOX 
12490 RETURN 

12300 NX»( TX * (l - ,23/DF) l/SQRG ♦ ,3*TXC2/DF) 

12310 MU-0 1 VAR-1 112-1 1 GOSUB 6000 
12320 FX-TX 
12330 AT-2*FOX - 1 
12340 RETURN 
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F - DISTRIBUTION 

13000 REM *** F-DISTRIBUTION SUBROUTINE 

13010 IF 12*1 GOTO 13110 

13020 CLS I PRINT 020,'"F-DISTRIBUTION* 

13030 PRINT 0128,-?: INPUT"Enter F-Statistic x "IFX 

13040 IF FX<0 THEN PRINT0128, CHRS12S3),'SPRINT 0128,"F-STATISTIC MUST BE A POSITIVE NUMBER“;:FOR I 

-1 TO 13OO:NEXT:PRINT0128,CHR*(253)!:GOTO 13030 

13050 INPUT"Enter Degrees of Freedom in Numerator'iDN 

13060 IF DN>0 AND DN-INTtDH) GOTO 13080 

13070 PRINT"** ERROR •« ENTER AN INTEGER > O’lGOTO 13050 

13080 INFUT*Enter Degrees of Freedom in Denominator'iDD 

13090 IF DD>0 AND DD*INT(DD) GOTO 13110 

13100 PRINT"** ERROR •* ENTER AN INTEGER > 0“!GOTO 13080 

13110 PX*FX J IF FX*0 THEN FOX*0 ! GOX-1 5 RETURN 

13120 IF DN > 60 AND DD > 60 GOTO 13720 

13130 REM «*• APPROX. FOR SMALL TO MODERATE D. F. ♦*• 

13140 K1*DN ! K2-DD i PX*FX 

13130 A*DN*FX/(DN*FX+DD) ! Al-t-A 

13160 IF Al< IE-37 THEN A1-1E-37 

13170 D1-DN*,3 : D2*DD«.3 l D3-D1+D2-1 

13180 R-o: si*o : S2*o: dl»i : c*.23 

13190 XM«*1 1 XK*1 I FI-3.141393 ! KN-DD 

13200 KM-INTID2)*2 

13210 IF KKOKK GOTO 13390 

13220 KN-D2-1 

13230 IF KN-0 GOTO 13320 

13240 FOR I«1 TO ICN 

13250 SI-DL*S1*R 

13260 D2-D2-1 

13270 D3-D3-1 

13280 TEM-A1/D2 

13290 R*D3*TEM 

13300 S2-(R*TEM)*S2 

13310 NEXT I 

13320 Sl-DLtSl*R 

13330 DL-0 

13340 T—1 

13330 D3»-l 

13360 S2-A»S2 

13370 C*C*4 

13380 GOTO 13370 

13390 KM-IKT(D2> 

13400 IF KN-0 GOTO 13490 
13410 FOR 1*1 TO KN 
13420 S1-DL*S1*R 
13430 D2-D2-1 
13440 D3-D3-1 
13430 TEM-A1/D2 
13460 R*D3»TEM 
13470 S2-(R>TEM)«S2 
13480 NEXT I 
13490 S1*XK*S1 
13300 S2*XK«S2 
13310 ART*SQR1A 11 
13320 XM-XM*ART 
13330 T»(XM-ART>/A1 
13340 D3»-J 
13530 XK-2/FI 
13360 C-C *2 

13370 IF 0-873 GOTO 13660 
13S90 D2*Dl 
13390 D3-D2+D3 
13600 S2-SI 
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13610 Sl-0 
13620 Al-A 

13630 IF AK IE-37 THEN A1-1E-37 
13640 KN-DN 
13650 GOTO 13200 

13660 IF CO.125 THEN DL-4/PI ♦ ATN1T) 

13670 FQX-XM*(S2-S1)-DL 
13600 IF FOXCO THEN FOX-O 
13690 IF FOXM THEN FOX-1 
13700 QOX-l-FOX 
13710 RETURN 

13720 REM *** APPROX. FOR LARGE D. F. *** 

13730 NX-(«FXClt/3»*<l-12/(9*DD»H-<l.-i2/<9*DN»)>/SGR<(2/(9*DN»+<FXt<2/3»*2/<9»DD» 

13740 MU-0 1 VAR-1 112-1 

13750 GOSUB 6000 

13760 PX-FX 

13770 RETURN 


INVERSE NORMAL 

14000 REM ••• INVERSE NORMAL SUBROUTINE 

14010 IF 12*1 GOTO 14100 

14020 CLS l PRINT 320, "INVERSE NORMAL* 

14030 PRINT 5 PRINT 1 INFUT’Entar Maan’iMO 
14040 INPUT *Entar Varianca‘!VR 

14050 IF VRC-0 PRINT"** ERROR •• "!VR{"IS NOT A VALID VARIANCE’IGOTO 14040 
14060 INPUT'Entar Probability li.a. F(x) * Prob, 1 X <■ x ) - "JPR 
14070 IF PR-0 THEN XQ-9E-36 ELSE IF PR-1 XQ-9E+36 
14060 IF PR-0 OR PR-1 RETURN 

14090 IF PR<0 OR PR>1 THEN PRINT"** ERROR *« "JPRJ'IS NOT A VALID PROBABILITY" (GOTO 14060 
14100 IF PR>.5 THEN RL*SQR1L0GU/«1-PRH2>» ELSE RL-SQR(LOG( 1/1PRC2M 

14110 XQ-HL-«2.313517+l.a02853*RL>+<.010328»RL[2»/(i+«1.4327B8»RL>+I.189269*RLt2>*(.001308*RL£3»> 
14120 IF PR>».5 THEN XQ-MO*(SQR(VR)*XQ) ELSE XQ*M0-(SGR(VR1*XQ) 

14130 RETURN 


INVERSE CHI - SQUARE DISTRIBUTION 

15000 REM INVERSE CHI-SQUARE SUBROUTINE 

15010 IF 12-1 GOTO 15090 

15020 CLS 1 PRINT 320,‘INVERSE CHI-SQUARE" 

15030 PRINT 1 INPUT" Ent»r Degrea* of Freadom'JDF 

15040 IF DF<-0 OR DFOINT(DF) THEN PRINT"** ERROR ** ENTER AN INTEGER > O’lGOTO 15030 
15050 INPUT "Entar Probability i.a. P(x) - Prob. ( X <- x 1 • “,FR 
15060 IF PR-0 THEN XQ-0 1 RETURN 
15070 IF PR-1 THEN XQ-9E36 ! RETURN 

15080 IF PR<0 OR FR>1 THEN PRINT"** ERROR «« "IPRl'IS NOT A VALID PROBABIUTY'tGOTO 15050 
15090 MO-O I VR-t 112-1 J PI-PR 
15100 IF DlOl GOTO 15150 
13110 PR>.5*<1-PI» 

15120 GOSUB 14000 
15130 XQ-XQI2 1 PR-PI 
15140 RETURN 
15150 IF DF>2 GOTO 13180 
13160 XQ—2*L0G(1-P» 

13170 RETURN 
15180 GOSUB 14000 
13190 DQ-XQ 

13200 IF DF><2*INT(4*ABStDQ)») GOTO 15300 

13210 X9-1<0000886 - .2237368/DF - .01513904/DFI2 

15220 X9-X9+1DFI-.3 * DQ • 1.4713941 ♦ .02607083/DF - .008986007/DFI2H 

13230 X9-X9 ♦ U/DF * DQC2 • 1.0001348028 * ,01128186/DF * .02277679/DFI21) 

13240 X9-X9 ♦ (OFC-1.3 * DQI3 • (-.008553069 - .01133761/DF - .01323293/DF(2» 
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13250 X9-X9 ♦ (DF£-2 • DGC4 * (.00312338 ♦ ,003169654/DF - .006930336/DFC2)) 
13260 X9-X9 (DFC-2.5 ♦ D0t3 * (-.0008426812 + .00253001/DF ♦ .001060438/DF£2» 

13270 X9-X9 + (DFt-3 * DQC6 * (.00009780499 - .001430117/DF ♦ .001365326/DFC2H 
13280 XQ=DF*X9£3 
13290 RETURN 

13300 X9»l - 2/(9»DF) ♦ (4*DQ£4 ♦ 16*DOt2 -28)/(1215*DF£2) 

13310 X9-X9 ♦ <8»DQC 6 + 720*DQC4 ♦ 3126*DQ£2 + 2904)/(229633*DFE3l 
13320 X8-DQ/3 ♦ (-DGC3 ♦ 3*DQ)/(162»DF) 

13330 XB-XB - (3*DQC5 ♦ 40*DQC3 ♦ 45*DQ)/t5832*DFt2> 

15340 X8-X8 ♦ (301»DQC7 - 1517#DQC5 -32769*OQC3 - 79349*DQ)/(7873200*DF£3> 
13330 X8- SOR(2/DF) * X8 
15360 XQ-DF ♦ ((X8+X91C3) 

15370 RETURN 


IWYERSE -1 

16000 REM • *• IN VERS£-t SUBROUTINE 
16010 IF 12*1 GOTO 16090 

16020 CLS I PRINT G16, *INVER5E-t PROBABILITY DISTRIBUTION* 

16030 PRINT J INPUT“£nter Degrees of Fraedom“iDF 

16040 IF DFC-0 OR DFOINT(DF) THEN PRINT*** ERROR ** ENTER AN INTEGER > 0‘IGOTO 16030 
16030 INPUT*Enter Probability i.e. F<xl - Prob. ( X <« x ) » *JPR 
16060 IF FR-1 THEN XQ-9E36 ELSE IF PRO XQ-9E-36 
16070 IF PR-1 OR PRO RETURN 

16080 IF PR<0 OR PR>1 THEN PRINT”** ERROR ** *iPR;*IS NOT A VALID PROBABILITY"GOTO 16030 
16090 moo : vr-1 : 12*1 
16100 GOSUB 14000 
16110 XC-XQ 

16120 XQ-(.25 * (XCC3 * XCD/DF 

16130 XQ-XG) ♦ t(l/96) « (5«XCC3 ♦ 16*XC£3 ♦ 3*XC»/DFE2 

16140 XQ-XQ ♦ ((1/384) « (3*XCt7 * 19*XC£3 ♦ 17*XCE3 -I3*XC»/DF£3 

16130 XQ-XQ ♦ ((1/92160) • (79*XCE9 * 776*XCC7 ♦ 1482*XCE3 -1920*XCE3 -943«XC))/DF£4 

16160 XQ- XC * XQ 

16170 IF DF>3 RETURN 

16100 IF PR>.05 AND PR<.95 GOTO 16230 

16190 IF DF>1 GOTO 16210ELSE IF PR>.999 OR PR<.001 THEN DELTA-1000 ELSE DELTA-100 
16200 GOTO 16240 

16210 IF DF-2 AND PR>.99 THEN DELTA-100 ELSE DELTA -10 
16220 GOTO 16240 
16230 DELTA-1 

16240 IF PR>-.3 THEN TX-XQ ELSE TX—XQTX-XQ 
16230 IF PR>-J THEN PE-PR ELSE FE-l-PR 
16260 GOSUB 12000 

16270 IF ABS(PE-FOX) 5- .00001 THEN GOTO 16290ELSE IF PR-PE THEN XQ-TX ELSE XO—TX 
RFTtTPM 

16290 IF PE < FOX THEN TX-TX-DELTA 
16300 IF PE < FOX THEN DELTA-DELTA/10 
16310 IF PE < FOX THEN GOTO 16260 
16320 TX-TX*DELTA 
16330 GOTO 16260 


71 



INVERSE? 


17000 HEM **• IN VERSE-F SUBROUTINE 

17010 IF 12-1 GOTO 17130 

17020 CLS ! PRINT @ 22, “INVERSE - F" 

17030 PRINT ! PRINT 

17040 INPUT"Enter Degrees of Freedom in Numerator";DN 

17050 IF DN>0 AND DN-INT(DN) GOTO 17070 

17060 PRINT"** ERROR ** ENTER AN INTEGER } 0-JGOTO 17040 

17070 INPUT"Enter Degrees of Freedom in Denomimator'JDD 

17080 IF DD>0 AND DD-INT(DD) GOTO 17100 

17090 PRINT"** ERROR •« ENTER AN INTEGER > 0"!GOTO 17070 

17100 INPUT’Enter Probability i.e. F(x) » Prob. II(=>il« "JPR 

17110 IF PR<0 OR PRM THEN PRINT"** ERROR *♦ *;PRJ*IS NOT A VALID PROBABILITY" iGOTO 17100 

17120 IF PR}. 9 99 PRINT “THIS SUBROUTINE DOES NOT CALCULATE QUANTILES CORRESPONDING TO FROBA 

BIUTIES GREATER THAN ,999“’.GOTO 17100 

17130 IF PR® l THEN PR».9999 

17140 IF PR-0 THEN XQ*0 ! RETURN 

17150 IF DN<20 OR DDC20 GOTO 17230 

17160 M0®0 t VR®1 S 12*1 

17170 G0SU8 14000 

17180 H - 2 * < <1/(DD-1) ♦ 1/<DN-I»t-1 ) 

17190 LM - (XGI2 - 3>/6 

17200 WX-<XQ • SQRtH * LMH/H - ( (l/(DN-li - l/(DD-l) > • < LM * 5/6 - 2/<3*H> » 

17210 XQ * EXPtt*WXJ 
17220 RETURN 

17230 IF DNOl AND DDOl OR DN»1 AND DD-1 GOTO 17330 

17240 FH-FR ! IF DN-1 THEN PR-<l*PR>/2 ELSE PR-1-PR/2 

17250 IF DN*t THEN DF-DD ELSE DF-DN 

17260 12-1 ! GOSUB 16000 

17270 XQ-XQ*XQ i IF DD*1 THEN XQ*1/XQ 

17280 PR-PH 

17290 IF PR7-.99 THEN XQ*.994*XQ S GOTO 17320 

17300 IF FR}-.975 AND DN>3 THEN XQ-.993«XQ i GOTO 17320 

17310 IF PR7-.93 AND DN73 THEN XQ*.996*XQ 

17320 RETURN 

17330 IF PR}. 999 AND DD*1 THEN XQ*9E36 ! RETURN 
17340 IF PR}«.997 AND DD-1THEN DELTA-500000 5 GOTO 17410 
17350 IF PR>.993 AND DD-1 THEN DELTA-15000 i GOTO 17410 
17360 IF PR>-.99 AND DD-1 THEN DELTA-5000 i GOTO 17410 

17370 IF PR>->975 AND DD-1 OR FR}-.999 AND DD-2 THEN DELTA-1000 t GOTO 17410 

17380 IF PR}-.9 AND DD-1 OR PR}.99 AND DD02 OR PR}-.999 AND DDC-3 THEN DELTA-100 ! GOTO 17410 

17390 IF PR}-.99 AND DD<-3 THEN DELTA-30 1 GOTO 17410 

17400 IF PR}-.975 AND DD<*2 THEN DELTA-50 ELSE DELTA-10 

17410 LL-0 

17420 FX-DELTA 

17430 12-1 ! GOSUB 13000 

17440 IF ABS(FR-FOXK».00001 THEN XQ-FX 1 RETURN 

17450 IF PR < FOX THEN DELTA-DELTA/2 ! FX-DELTA S GOTO 17430 

17460 LL-FX S FX-FX*DELTA ! UL-FX 1 GOSUB 13000 

17470 IF ABS(FR-FOXX®.00001 THEN XQ-FX 5 RETURN 

17480 D ELTA-IUL-LL1/2 

17490 FX-LL+DELTA 

17500 GOSUB 13000 

17510 IF ABS<PR-FOXK-.00001 THEN XQ-FX S RETURN 
17520 IF PR } FOX THEN LL-FX ELSE UL-FX 
17530 GOTO 17480 
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ggaS ES StSPE HTEAU 


10000 REM * * • INVERSE EXPONENTIAL SUBROUTINE 
18010 CLS : PRINT 520,“INVERSE EXPONENTIAL* 

18020 PRINT : PRINT* You may either enter the exponential parameter (LAMBDA) or* 

18030 PRINT'the exponential MEAN (1 /lambda).* 

18040 KP=320 ! L*»“l* ! R*«*m* ! GOSUB 600 
18050 ZEP’ZQ 

18060 IF ZEP=’l* OR ZEP“*L* THEN ZEP»*LAMBDA" ELSE ZEP-*M£AN" 

18070 PRINT I PRINT* Enter *;ZEPi 
18080 INPUT LAMBDA 

18090 IF LAM <* 0 THEN PRINT*** ERROR ** ENTER A VALUE > OVGOTO 18070 
18100 PRINT ! INFUT’Enter Probability i.e. Prob.l T C« t ) «*5FR 
18110 IF PR*0 THEN TI=0 i FR»1 ! GOTO 18160 

18120 IF PR-1 THEN PRINT ! PRINT’For Prob. < T O t) ■ 1.0 t - infinity*iRETURN 

18130 IF PRCO OR PR>1 THEN PRINT*** ERROR ♦* *iPR}*IS NOT A VALID PROBABILITT*!GOTO 18100 

18140 PB=1-PR 

18150 IF ZEP»*LAMBDA* TI- - LOG(PB)/LAMBDA ELSE TI ■ - LOG(PB)»LAMBDA 
18160 PRINT ! PRINT’For Prob. I T (» t I **iPRi* t -’HI 
18170 RETURN 


END OF PROGRAM »1 
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PROGRAM 


MENU 

10 REM •••MAIN MENU*** 

20 CL5 ! PRINT TAB120) “*** MENU •*•* 

30 PRINT 

40 PRINT*!!) Sample Statistics (mu, sigma, gl, g2) 

SO FRINT"<2) Seven Number Summary 

AO PRINT"(3) Histogram, Stem & Leaf, Probability Plots 

70 PRINT*<4) Return to Program Selection Menu 

80 PRINT 0445, 1 'NUMBER OF OPTION DESIRED “iCHR*(123)iCHR*<176)!CHRS(123) 

90 ZI-INKEY* ! IF 21="* GOTO 90ELSE IO=VAL(2I) 

100 PRINT 0471,Zi;: FOR 1=1 TO 100? NEXT 

110 IF IO<l OR 1054 THEN PRINT 064O,CHR»<253) ! PRINT 0440,"**« NOT A VALID OPTION •***! FOR I«1 TO 

1000 ! NEXT ! PRINT 0440, CHRS(253) GOTO 80 

120 ON 10 GOSUB 500, 1000, 2000, 10000 

130 GOTO 10 

140 CLS 

150 PRINT*Option"iIO;*is not applicable since only sample statistics were entered.* 

140 PRINT094O,‘Press SPACE BAR to return to main menu*; 

170 PRINT 'Kurtosis »‘!G2 
180 IF INKEYtO* * GOTO 180 
190 CLS 
200 GOTO 10 


DISPLAY SAMPLE STATISTICS 

500 REM *** SAMPLE STATISTICS SUBROUTINE 
310 CLS 

320 PRINT TAB(23)"SAMPLE STATISTICS* 

530 PRINT 0192,"Sample Mean «*5XB 

540 PRINT ‘Sample Variance (unbiased estimate) =*!V1 

550 PRINT-M.L.E. for Variance =‘}V 

340 PRINT *Std< Deviation (sqr, root of sample var.) **1SQR(V1) 

570 PRINT 'Skewness »*JG1 
580 PRINT ‘Kurtosis -*;G2 

590 PRINT * PRINT-NOTE: The unbiased estimate of variance * 

400 PRINT* *;CHRS(123)S* i*. SUM (X(i> - Xbar)t2 / (n-1) ‘;CHRS(123>i* is used in- 
410 PRINT* all subsequent calculations throughout the package,* 

420 PRINT 0940,‘PRESS SPACE BAR TO RETURN TO MAIN MENU*? 

430 IF INKEY* O * * GOTO 430 
440 RETURN 


SEVEN NUMBER SUMMARY 

1000 REM ••• SEVEN NUMBER SUMMART SUBROUTINE 
1010 CLS 

1020 PRINT TABU D’SEVEN NUMBER SUMMARY- 
1030 MD=INT(J*<IN+1» ! QD=INT(,3*<MD+1)) : ED=INT(.3*«3D+1» 

1040 IF IN/2=INT(IN/2) THEN M=.5«(Y(MD)+Y(MD+1» ELSE M=Y(MD> 

1050 IF QD=J*(MD+1) THEN QL»Y(0D) ! QU=Y(IN-QD+l> 5 GOTO 1070 
1040 QL=J*(Y(QD)+Y<QD+1)) 5 QU«*5*<Y<IH-QD)+YCIN-QD+1)) 

1070 IF ED=.5#(QI>1> THEN EL=Y(ED> 5 EU-YGN-ED+l) 1 GOTO 1090 
1060 EL=.5*(Y(EDHYIED+1» t EU*.5*(Y<IN-EDM-YUN-ED+1)> 

1090 TV»M 1 GOSUB 1250 : PRINT 0144,‘MEDIAN « *il PRINT USING P*!M5 
1100 PRINT 0237,‘MEASURES OF*? 

1110 PRINT 0247,“LOWER MEAN U PP ER SKEWNESS- 

1115 PRINT 0331,*- - - -* 

1120 PRINT 0384,-QUARTER:* t TV*QL ! GOSUB 1230 : PRINT 0394,**::PRINT USING P*5QLJ 
1130 TV=J*<QLrQU) ! GOSUB 1250 : PRINT 04O3,**{:PRINT USING PS I TV[ 
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1140 TV-QU : GOSUS 1250 ! PRINT 04|6,-"i:PRINT USING P»;QUi 

1150 TV* (.5»(QU+0L> - M)/(QU-QU : GOSUS 1250 ! PRINT 0431,"“!:PRINT USING PSJTV! 

1153 PRINT 0512,'EIGTH!"!! TV=EL : GOSUB 1230 ! PRINT 0522,*"',:PRINT USING PSiEL! 

1160 TV*i5*(EL*EU)! GOSUB 1250 ! PRINT 0S33."*;:PRINT USING PSJTV! 

1165 TV*EU ! GOSUB 1250 : PRINT 0S44,**;:PRINT USING PSJEU! 

1170 TV*I.5«<EL»EU) - MJ/tEU-EL) ! GOSUB 1250 ! PRINT 0559,*";:PRINT USING P*?TV; 

1175 PRINT 0640,'EXTREMEi"! ! TV*Y<1> : GOSUB 1250 : PRINT 065O,"*i:PRINT USING PS! TUI! 

1180 TV-.5t<Y<ll*Y<IN» : GOSUB 1250 : PRINT 0661,"";:PRINT using ps;tv; 

1190 TV-YCIN) : GOSUB 1250 ! PRINT 0672,“"SPRINT USING PV,Y(IN)S 

1200 TV»<,3*<Y(INHYID) - M)/<Y(IN>-Ytl» ! GOSUB 1250 : PRINT 0687,"";:PRINT USING P*iTV 
1210 PRINT 0896,"Press SPACE BAR when ready to return to program menu*! 

1220 IF INKEY40* * THEN 1220 
1230 RETURN 
1250 HB*ABS(TV) 

1255 PI S<1 )»*♦.*•••" : F1S<2>*"##.###“ J Pl*<3>-"•*•.#*“ Pl*<4)-" •♦••>•* ! Pl*<5>*"«****" 1 P1»<6>-*M** 

cctc- 

1260 IF HB7-.0001 AND HB<1 OR HB*0 THEN PS-P1SU) 1 RETURN 
1270 IF HB>-1 AND HB<10 THEN PS=P1S<2> t RETURN 
1280 IF HB>*10 AND HBOOO THEN P4=P1»(3) ! RETURN 
1290 IF HB>»100 AND HBC1000 THEN P4«P1*(4) ! RETURN 
1300 IF HB3-1000 AND HBC10000 THEN P*-P1»(5)! RETURN 
1310 PS*P1S(6) 1 RETURN 


VSSU QF PATA.&Q.TS 

2000 REM • DATA "PLOTS SUBROUTINE 
2010 CLS 

2020 PRINT TAB(23)*D ATA PLOTS" 

2030 PRINT ! PRINT 

2040 PRINT "MENU!" ! PRINT* tl) Histogram 
2030 PRINT" Q) Stem & Leaf Plot 
2060 OC*CI 

2070 PRINT* (3) Normal Probability Plot 
2080 PRINT" (4) Exponential Probability Rot 

2090 PRINT 0643,"NUMBER OF OPTION DESIRED “;CHR»(123)!CHRS(176);CHRS(125> 

2100 a-INKEY* ! IF II*"" GOTO 2100ELSE IO*VAL(ZI) 

2110 PRINT 0671,a;: FOR >1 TO 100 : NEXT 

2120 IF 10<l OR IO>4 PRINT 0640, CHR*(230>i! PRINT 0640,"«** NOT A VALID OPTION aea'JSFOR 1*1 TO 1000! 
NEXT ! PRINT 064O,CHR*<253>! GOTO 2090 
2130 IF IO>*3 GOTO 2200 

2140 PRINT 0640,“Min Value **‘,Y(1)!* Max Value **rY(IN)!PRINT" Range *'"!YON>-YU > 

2150 PRINT t PRINT"AFTER VIEWING "j! IF 10*1 PRINT*HISTOGRAM"! ELSE PRINTTTEM & LEAF"! 

2160 PRINT* PRESS SPACE BAR TO CONTINUE"! * 

2170 IF INO30 OR 10-2 THEN FOR 1-1 TO 2000 5 NEXT 
2180 SC-lOtONT<LOG<Y<INKYU>>/LOG<10>>) 

2190 NB*0 ! CI*1 

2200 ON 10 GOSUB 3000, 4000, 5000, 6000 
2210 RETURN 


HISTQCRAM 

3000 REM ••• HISTOGRAM SUBROUTINE 
3010 NB-0 

3020 IF SC<1 THEN GF«INT<Y<1)/SC> ELSE GF-INT(Y<1)/INTISC)I 
3030 IF SC<1 GOTO 3060 

3040 IF GF ♦ N8*CI < INTlYttNl/INTfSCI) ♦ IE-6 AMD KB <- 13 THEN NB*MB*1 ! GOTO 3040 
3030 GOTO 3070 

3060 IF GF ♦ NB*CI C aHT<Y(IN)/SC) ♦ IE-6) AND NB C- IS THEN NB-NB-M ! GOTO 3060 

3070 PC-NB/15 5 IF PCM THEN CI*CI*1 ! KB-0 ! GOTO 3020 

3080 IF BF*O"0N" AND PCC1/3 THEN CI-CI/2 I NB*0 l GOTO 3020 

3090 n-i 

3100 FOR I-l TO NB ! NIUV-0 

3110 FOR J*n TO IN ! IF Y(J)/SC < IaCI ♦ INTrr<l>/SC) - IE-6 THEM NND-NKIW1 ELSE H-J ! GOTO 3130 
•3120 NEXT T 
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3130 IF II > IM GOTO 3150 
3140 NEXT I 

3150 MN*0 i FOR I-l TO NB i IF NI(I) > MN THEN MN-NKI) 

3160 NEXT I 

3170 IF MN < 11 THEN VS*. 25 : GOTO 3200 

3180 IF MN < 22 THEN VS-.S : GOTO 3200 

3190 IF MN <* 45 THEN VS-1 ELSE VS-INTIMN/43I+1 

3200 IF SC < I GOTO 3260 

3210 CLS ! CN*INTiY( 1 l/INTISC)) ! FOR 1*0 TO 4*<N8> STEP 4 
3220 TN-CN-INTICN) l IF TN*0 OR TN**5 THEN PRINT 0961+1, CHS 
3230 CN*CN+CI 

3240 IF CN/SC > 15»CI + INT(T(1)/SC> GOTO 3310 
3250 GOTO 3300 

3260 CLS ! CN=INTnrU)/5C> t FOR I * 0 TO 4*NB STEP 4 

3270 TN-CN-INT(CN) ! IF TN=0 OR TN=>.3 THEN PRINT 0961+1, CNJ 

3280 CN=CN+CI 

3290 REM IF CN > 15«CI ♦ INTtYdl/SC) GOTO 1882 
3300 NEXT I 

3310 IF VS>1 THEN LN-2 ELSE LN»0 
3320 FOR 1=0 TO 14 

3330 IF VS>1 THEN PRINT 0896-1*64, MID*<STR*<VS*LN>,2,3>+*-'ii LN-LN+3 ! GOTO 3370 
3340 IF 1-2 THEN K-2 ELSE IF 1*6 THEN K-5 ELSE IF 1-10 THEN K*8 ELSE IF 1*14 THEN K-U 
3350 IF VS-.5 THEN K=2*K 

3360 IF 1*2 OR 1*6 OR 1*10 OR 1-14 THEN PRINT 0896-1*64,MID*<STR*(K>.2,2>+'-*; 

3370 NEXT I 

3380 Xl*9 ! FOR 1-1 TO NB 1 IF NHD-0 GOTO 3440 

3390 JU-NKD/VS - 1 1 EC-NHD/VS-JU 1 IF EC<J5 THEN GOTO 3410 

3400 IF VS>2 THEN JU-JU* l 

3410 FOR J*0 TO JV ! Y1-44-J i IF Y1<0 THEN Y1=0 

3420 SETtXl.Yl) i SET(Xl-t,Yl> 1 SET(X1-2,Y1)! SET(XI+2,Y1) 5 SET(X1+1,Y1)! SET(X1-3,Y1) ! IF XK124 THEN S 
ET(X1+3,Y1> 

3430 NEXT J 

3440 Xl-Xl+8 ! NEXT I 

3450 IF SCOl THEN PRINT038,'Actual Values are*;VAL<STR*lSC)>;: PRINT 0102,'times horizintal values'! 

3460 IF INKEY* O " * GOTO 3460 

3470 CLS ! PRINT'OPTIONS!" ! PRINT* <1) Return to Menu* 1 PRINT* <2J Redraw Histogram with larger dass in 
terval* 1 IF MB<=7 PRINT* (3) Redraw with smaller class interval* 

3480 PRINT 0453,"NUMBER OF OPTION DESIRED -;CHR*U23>;CHR*<176>5CHR*U25) 

3490 XI-INKEY* ! IF II*" GOTO 3490ELSE IO-VALOI) 

3500 PRINT 0479,115! FOR 1*1 TO 100 ! NEXT 
3510 IF NB37 THEN IU*2 ELSE IU*3 

3520 IF IOC1 OR IO>IU PRINT 0448,CHRS<253> ! PRINT 0448,'•** NOT A VALID OPTION ***'! FOR 1-1 TO 100 
0 1 NEXT ! PRINT 0448,CHR*<253): GOTO 3480 
3530 IF 10-2 THEN BF*-*ON' ELSE SF*-“OFF* 

3540 IF 10-1 THEN RETURN ELSE IF 10*2 THEN CI*2*CI ELSE CI-CI/2 

3550 CLS ! PRINT*AFTER VIEWING HISTOGRAM, PRESS SPACE BAR TO CONTINUE'; 

3560 NB*0 ! GOTO 3020 


gSM.ftLSAE 

4000 REM ••• STEM t LEAF SUBROUTINE •** 

4010 CLS ! IF SCOt THEN PRINT'ACTUAL VALUES ARE'SCi'TIMES THOSE SHOWN* 
4020 SP*INTfY<l)/SC+lE-6) 

4030 IF SP<0 THEN SP-SP-CI 
4040 EP*INTCY(IN)/SC) + Cl 

4050 IF ABS(EP-SP)/CI<*3 THEN CI-CI/2 ! GOTO 4010 
4060 16-0 ! JC*l ! IP*'* 

4070 FOR CO-SP TO EP STEP Cl 

4080 IF CO < 0 THEN CJ*INT(C01+1 ELSE CJ-INT(CO) 

4090 IF CJ * 0 AND CO < 0 GOTO 4120 

4100 PRINT USING*+#**;CJ; 5 PRINT CHR*(124)} 

4110 GOTO 4130 

4120 PRINT* -0*;CHR*U24)! 

4130 FOR J-JC TO IN 

4140 IF YU) O INTI Y<1>) GOTO 4170 
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4130 IF Ylll/SC > CO»CI GOTO 4350 
4160 GOTO 4160 

4170 IF YtD/SC ♦ IE-6 >= CO +CI THEN GOTO 4330 

4180 IF YtJ) >= 0 GOTO 4230 

4190 IF Y(J) O INT (Y(J)i GOTO 4230 

4200 IF Y(J)/SC - .001 <= CO+1 AND Y(J)/SC ? CO - Cl ♦ 1.001 THEN SL=<10*<Y(J)/SC - INTCY<J>/SC*lE-4»> ELS 

EGOTO 4310 

4210 IF 16=0 THEN 16=1 

4220 GOTO 4250 

4230 IF Y<J)/SC+TE-6 >» CO AND Y(J)/SC+lE-6 < CO*CI THEN SL=(10*<Y(J)/SC -INT(YtJ)/SC*lE-6») ELSE GOT 
0 4310 

4240 IF 16-0 THEN 16=1 

4230 IF YtJKO AND SL O 0 THEN SL-10-SL 

4260 IS=MID*< STR*(SL) ,2,1) 

4270 IF ZS»“." THEN ZS«*0* 

4280 IF COCO THEN IF=IS + ZP ELSE ZP=ZP*ZS 

4290 JC=JC*I ! IF JC > IN THEN JC-IN 

4300 GOTO 4330 

4310 IF 16=1 THEN 16=2 

4320 IF 16=2 GOTO 4350 

4330 NEXT J 

434016-0 

4330 PRINT ZP 5 ZP=** ! 16=0 1 NEXT CO 

4360 PRINT096O,‘PRESS SPACE BAR TO CONTINUE*; 

4370 IF INKEY* O * * GOTO 4370 

4380 IF ABS«EP-SP)/CID8 THEN RETURN ELSE CLS! PRINT'Do you want to see the Stem 6 Leaf with a smaller 
Class Interval* 

4390 PRINT 0197,CHR*(123>!‘ Y or N *;CHR*(125);* *;CHR*<123);CHR*U76);CHR*<12S> 

4400 H-INKEY* ! IF ZI=‘* GOTO 440OELSE PRINT 0210,HI! FOR 1=1 TO 100 ! NEXT 
4410 IF ZI=“y* OR ZI=*Y* THEN CI=CI/2 ! GOTO 4000 
4420 IF ZlO'n* AND ZIO*N* THEN GOTO 4390 
4430 RETURN 


NORMAL PROBABILITY PLOT 

3000 REM *** NORMAL PROBABILITY PLOT SUBROUTINE 
3010 CLS 

3020 PRINT-NORMAL PROBABILITY PLOT OF!* 

3030 PRINT" (II Every data point 
3040 PRINT* (2) Every other data point 
3030 PRINT* <31 Every Nth data point 

3060 PRINT 0643,“NUMBER OF OPTION DESIRED *;CHR*(123)«HR*(176)X:HR*(123) 

3070 □"INKEY* ! IF ZI=“* GOTO 3070ELSE IO»VAL(2I> 

5080 PRINT 0671,□!*. FOR 1=1 TO 100 ! NEXT 

3090 IF IOC1 OR IO>3 THEN PRINT 0640, CHR*(230)!J PRINT 0640,**** NOT A VALID OPTION ****JJ FOR 1=1 TO 
10001 NEXT 1 PRINT 0640, CHR*(250);i GOTO 3060 
3100 IF 10=1 OR 10=2 J=IO 1 GOTO 3120 
3110 PRINT @709,**;:iNPUT-Enter N"JJ 
3120 CLS 1 PRINT 064,'‘+3 -*5 
3130 PRINT 0192,"*2 -*{ 

3140 PRINT 0320,'‘*1 -*5 
3130 PRINT 0448,* 0 -’} 

3160 PRINT 0576,*-l -*! 

3170 PRINT 0704,*-2 -*5 
3180 PRINT 0832,‘-3 -*; 

3190 FOR I-l TO 43 
3200 SET (3,1) 

3210 NEXT I 
3220 FOR 1=3 TO 123 
3230 SETd,43) 

3240 NEXT I 

3230 PRINT 020,"NORMAL PROBABILITY PLOT*? 

3260 SP=INT(Y(1)I 1 IF Y(IN> ■ INT (Y(IN)I THEN EP-YCIN) ELSE EP=INT(TdN)Hl 

3270 PRINT 0962, SP1 

3280 I« INT( LOG<ABS(EF»/LGG<10> 1 
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3290 PRINT QlOzl-I, SIH*(EP»; 

3300 RG=EP - SP 

3310 MP- SP ♦ RG/2 ! IF MP>10 THEN MF-INT(MP) 

3320 PRINT 0990 - INT(I/2), MPJ 

3330 12-1 :mo=o : VR-1 

3340 FOR I»1 TO IN+1 STEP J 

3350 IF I5IN THEN I-IN 

3360 GV- 125 - INT( (EP - YCI))/RG ♦ 120) 

3370 PR® I/tIN+l) 

3380 GOSOB 7000 

3390 IF XQ < 0 GOTO 5420 

3400 DW - 22 - INT( CXQ/3.5) » 22 ) 

3410 GOTO 5430 

3420 DW* 22 ♦ INT ( (ABS(XQ)/3.5) « 22) 

3430 SET(OV.DW) 

3440 IF I-IN GOTO 3460 \ 

3430 NEXT I 

5460 IF INKEY4 O " * GOTO 5460 
3470 RETURN 


EXPONENTIAL PROBABILITY PLOT 

6000 REM *•• EXPONENTIAL PROBABILITY PLOT • ** 
6010 CLS 

6020 PRINT 064,-6 -* 

6030 PRINT 0192,“5 -* 

6040 PRINT 0320,-4 -* 

6030 PRINT 0448,*3 -* 

6060 PRINT 0376,-2 -* 

6070 PRINT 0704,-1 -• 

6080 PRINT 0832,-0- 
6090 FOR I-l TO 39 
6100 SET<3,I) 

6110 NEXT I 

6120 PRINT 0835>STR1NGS<6O,*,-)! 

6130 PRINT 020,-EXPONENTIAL PROBABILITY PLOT*; 

6140 CI-YCIN)/20 

6130 HS - 10 C (INT CLOG(CI*13)/LOG<10» -1) 

6160 Cl- INT (CI/HS) ♦ 1 
6170 FOR I»0 TO 19 
6180 ZH-STR«<CI*(I+l» 

6190 IF CI*<I+1) >-100 AND INTa/2) • 1/2 GOTO 6210 

6200 PRINT 0898+1*3,ZHi 

6210 NEXT 1 

6220 PRINT 0990, -X*5 

6230 IF HS - 1 THEN GOTO 6240ELSE PRINT * ••JHSi 
6240IG-0 
6230 IX-IN 

6260 FOR 1-20 TO 1 STEP -1 

6270 IF YCTX) > I*CI*HS THEN IG-IG+1 5 IX-IX-l : GOTO 6270 
6280 IF IG-0 OR IG-IH GOTO 6320 
6290 07-1*6 

6300 DW- 40 - INTt-LOGdG/GN+D) * 6) 

6310 SET <OV,DW> 

6320 IH-IG 
6330 NEXT I 

6340 PRINT 0133,‘Prtsi SPACE BAR to cantlnu*-; 

6330 IF INKEY* O * * GOTO 6330 
6360 RETURN 
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INVERSE NORMAL 


7000 REM **• INVERSE NORMAL SUBROUTINE 
7010 IF PR-.3 THEN XQ=MO ! RETURN 

7020 Cl*2.513317 J C2*.S028S3 i C3*,010328 ! C4*l. J C3-1.432788 J C6*.l«9269 ! C7-.001308 
7030 IF PH>.5 THEN RL=SQR(LQG(C4/((C4-PR)*(C4-PR)))) ELSE RL*SCIR(LOG<C4/<FH*FR)» 
7040XQ“RL-«Cl+<C2*RLH<C3*RL*RL»/<C4+<C3*RL>+<C6*HL*RLH<C7«RL*RLaRLm 
7030 IF PR} .3 XQ*MO*<SQR<VR)*XC» ELSE XO»MO-(SOR(VRHXQ) 

7060 RETURN 


PROGRAM SELECTION MENU 


10000REM •** PROGRAM SELECTION MENU •** 

10010 CLS 

10020 PRINT TAB(20)“FROGRAM SELECTION MENU’ 

10030 PRINT 0128,**1 Data I/O, Transformations, Prob. & Inv. Prod. Distributions "PRINT 

10040 PRINT"*2 Sample Stat's, 7-Num-Sum, Plots”!PRINT 

10050 PRINT**3 Parametric Tests of Hypotheses & Confidence Intervals"SPRINT 

10060 PRINT*#4 Robust Estimates and Confidence Intervals‘5PRINT 

10070 PRINT**3 Multiple Linear Regressjon-iPRINT 

10060 PRINT" 46 Analysis of Variance- 

10090 PRINT 0900,"NUMBER OF PROGRAM DESIRED *;CHRSU23)iCHR*(176»;CHRS(123> 

10100 n-INKEY* : IF a-- GOTO IOIOOELSE IO*VALai) 

10110 PRINT 0927,ZI!1 FOR 1*1 TO 100! NEXT 

10120 IF I0<1 OR IO>6 THEN PRINT 0896,CHR*<233K:PR1NT 0696,-*** NOT A VALID OPTION ***’;; FOR 1-1 T 

O 1000 : NEXT ! PRINT 0896,CHRS(253)i; GOTO 10090 

10130 IF 10*2 GOTO 10 

10140 IF IO-l THEN ZT**1" J GOTO 10160 

10130 IF 10*3 OR IO»4 THEN ZT=“2" ELSE ZT*’3* 

10160 CLSlFRINT’INSEKT TAPE *“tZTr INTO ESF-80 AND PRESS SPACE BAR* 

10170 IF INKEY* O * * GOTO 10170 

10180 IF 10*1 OR 10*3 OR 10*3 THEN 10*1 ELSE 10*2 

10190 QLOAD 10 


END OF PROGRAM *2 
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PROGRAM 4*3 


MENU 

10 REM •••MAIM MENU*** 

20 LS-CHRSU23) : RS=CHRS(125> ! BLS-CHRH232) 

30 CLS ! PRINT TAB<20) “•** MENU •«*■ : PRINT 

40 PRINT'TESTS OF HYPOTHESES & CONFIDENCE INTERVALS FOR DATA ASSUMED!* 

30 PRINT @256," It) Norm*! (or Control Limit Theorem Applies) 

60 PRINT" (2) Exponential 

70 PRINT" <3> Bernoulli 

60 PRINT S PRINT'OTHER OPTIONS!" 

90 PRINT ! PRINT" <41 Input Summary Statistics 
100 PRINT" (3) Return to Program Selection Menu 
110 JP-896 ! KS-3 ! GOSUB 760 

120 IF 10*2 AND XB<=0 GOTO 130ELSE IF IO-3 AND XB<0 OR IO-3 AND XB>1 GOTO 130 
130 ON IO GOSUB 1000, 4000, 6000, 300, 22000 
140 GOTO 10 
ISO CLS 

160 IF IO-2 THEN PRINT "EXPONENTIAL PARAMETERS MUST BE GREATER THAN ZERO* ! PRINT ! PRINT”THE 
SAMPLE MEAN « *;XB!‘AND THEREFORE DOES NOT CORRESPOND TO A" ! PRINT-EXPONENTIAL POPULATION 
"SGOTO 180 

170 PRINT "THE SAMPLE MEAN FROM A BERNOULLI POPULATION MUST BE A NUMBER BETWEEN 0 6 1"! PR 
IHT ! PRINT"THE SAMPLE MEAN « *JXB!"AND THEREFORE DOES NOT CORRESPOND TO A" ! PRINT"BERNOULL 
I POPULATION" 

190 PRINT 0896,"PRESS SPACE BAR TO RETURN TO MAIN MENU" 

190 IF INKEYSO" * GOTO 190 
200 GOTO 10 


INPUT SUMMARY STATISTICS 

300 REM *** INPUT SUMMARY STATISTICS •** 

310 CLS: 14-1 

320 PRINT TAB(20)*INPUT SUMMARY STATISTICS* 

330 PRINT 0128," If summary statistics are entered then any data which was previously entered and all samp 
le statistics calculated For such data will be cleared." 

340 PRINT ! PRINT* Do you wish to proceed?" 

330 T1S-"T" ! T2S»"N" 1 TP-448 ! GOSUB 710 
360 IF 2Q-"n" OR ZQ»"N" THEN 14-0 ! RETURN 
370 CLS ! PRINT-CLEARING DATA..." 

380 FOR 1-1 TO 200 ! YUJ-0 5 X<I)«0 ! NEXT ! IN-0 ! Vl-0 ! V«0 ! XB-0 ! CLS 
390 PRINT TAB(20 ("INPUT SUMMARY STATISTICS* 

400 PRINT ! INPUT’Enter Sample Sixe"iIN ! PRINT 0128, BLS 

410 PRINT 0128,"Is data assumed to be from an Exponential Population ?* 

420 T1*-*Y“ ! T2S-"N* ! JP-256 ! GOSUB 710! IF ZQ-"n" OR ZQ-"N* GOTO 500 

430 PRINT @128, BLS ! PRINT 0128,"Will you enter the MEAN or the Exponential Parameter LAMBDA?" 

440 T1S-"M* ! T2S-"L" ! GOSUB 710 

430 PRINT 0128, BLS ! PRINT BLS ! PRINT 0128,"Enter "J! IF Z0-*m" OR IQ-"M" PRINT "Mean"! ELSE PRINT "La 
mbda‘1 

460 INPUT XB i IF XSOO THEN PRINT 064," *♦* ERROR *** EXPONENTIAL PARAMETERS MUST BE GREATER 

THAN ZERO"!GOTO 450 

470 IF ZQ»*1* OR ZQ-"L" THEN XB-l/XB 

480 V1-XB*XB ! V-Vl 
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490 RETURN 

S00 PRINT Q128. BL* i PRINT 3128,‘Is the sample assumed to be from a Bernoulli Population?* 

510 JP*256 ! T1*»“Y“ ! T2*»‘N* ! GOSUB 710 
520 IF ZQ-“n“ OR ZQ=“N‘ GOTO 560 

530 PRINT 3128, BL* 1 PRINT BL* 1 PRINT 3120,‘Enter the number of occurrence* out of*!INJ*trials*inNPUT XB 
540 IF XB<0 OR XBMN GOTO 550ELSE IF XBOINT(XB) GOTO 550ELSE XB»XB/IN ! V1*XB#(1-XBI S V=*V1 i RET 
URN 

550 PRINT 8128, BL* ! PRINT 364, '*** ERROR ENTER AN INTEGER FROM 0 TO *!IN 1 GOTO 530 
560 PRINT 3128, BL* 1 PRINT BL* 1 PRINT 3128,‘Enter Sample Mean*i!INPUT XB 
570 PRINT 3256,“Do you want to enter Sample VARIANCE or STD. DEVIATION?* 

580 JP*320 ! T1*«"V" ! T2*»‘D* 1 GOSUB 710 
590 PRINT 3256, BL* ! PRINT BL* 

600 IF ZQ*“v“ OR ZQs-V* PRINT 3256,‘Enter Sample VariancB-JlINPUT VI { V«<(IN-n/IN)#Vl S RETURN 
610 PRINT 3256,‘Enter Std. Deviation‘!!INPUT VI i V1»V1*V1 i V*<(IN-1)/IN>*V1 5 RETURN 


OPTON SRSqTI Qtt 

700 REM *** OPTION SELECTION #*# 

710 PRINT 8JP+5,L*5" "ITUi* or *1T2*;* *JR*J* *JL*;CHR*<176);R* 

720 ZQ-INKEY* S IF ZQ*“ GOTO 720 

730 PRINT 3JPM0, ZQ J FOR IJ*t TO 100 1 NEXT 

740 IF ZQOT1* AND ZQOCHR*(ASC(Tl*H32J AND ZOOT2* AND ZQOCHR*(ASC(T2*)*32) THEN GOTO 710 
750 RETURN 

760 PRINT 3JP+5,‘NUMBER OF OPTION DESIRED ‘;L*5CHR*<176)!R* 

770 H-INKEY* i IF Z>** GOTO 770ELSE IO«VAL(ZI) 

780 PRINT QJP+31, ZIJ ! FOR I=»l TO 100 1 NEXT 

790 IF IO<l OR IO>KS THEN PRINT 8JP* BL* 1 PRINT 3JF,‘«*e NOT A VALID OPTION ****!! FOR I-l TO 1000 1 
NEXT i PRINT 3JP, BL* S GOTO 760 
800 RETURN 


MENU OF TESTS FOR NORMAL DATA 

1000 REM MENU OF TESTS FOR DATA FROM NORMAL DISTRIBUTION ♦ 

e * 

1010 CLS : FRINT*DATA ENTERED IS ASSUMED TO BE NORMAL 
1020 PRINT 1 PRINT-MENUS* 

1030 PRINT* 111 Test Hypothesis Conaming the Mean 
1040 PRINT* (2) Test Hypothesis Conaming the Variance 
1050 PRINT* (3) Confidena Interval for the Mean 
1060 PRINT* (4) Confidena Interval for the Variance 
1070 JT-640 S KS-4 S GOSUB 760 
1080 ON 10 GOSUB 1100, 2000, 3000, 3300 
1090 RETURN 


HYPOTHESES CONCERNING THE MEAN 

1100 REM See HYPOTHESES CONCERNING MEAN **» 

1110 CLS 

1120 PRINT 364,-HYPQTHESISt NULL ALTERNATIVE* t PRINT 

1130 PRINT* (1) mu a specified value mu */* specified value 

1140 PRINT* <2) mu O specified value mu > specified value 

1150 PRINT* (3) mu >■ specified value mu < specified value 

1160 JP-640 I ES>3 ! GOSUB 760 

1170 CLS : PRINT*Is the variance known?* 

1180 JP-64 ! T1S**Y“ ! T2S»*N* 1 GOSUB 710 
1190 IF Z0«*n“ OR ZQ«*N* GOTO 1250 

1200 PRINT‘Do you want to enter VARIANCE or STD. DEVIATION?* 

1210 JP-192 5 T1S»*V* i T2S**D* 1 GOSUB 7101 ZV*ZQ 

1220 PRINT’Enter ■{ i IF ZV«‘V* OR ZV-*v‘ PRINT ‘VARIANCE-1 ELSE PRINT’STD. DEVIATION*! 
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1230 INPUT AV 

1240 IF ZV»“D“ DR ZV-"d" AV-AV*AV 

1250 CLS ! INPUT“Enter the value that the mean Is hypothesized to be'iMZ 
12AO IF ZQ»"N* OR ZQ-*n* GOTO 1400ELSE ON 10 GOSUB 1280, 1480, 1630 
1270 RETURN 


mumBTHTOl 

1280 REM set NULL HYPOTHESIS 1 
1290 TZ-( SQRtIN) • (XB - MZ) >/SQR(AV) 

1300 PRINT : PRINTTest Statistic «*}TT 

1310 NX*ABS(TZ> i VAR-1 1 MU-0 112-1 1 GOSUB 10000 

1320 PRINT-P-Value - *S2*QOX 

1330 PRINTS INPUT-Enter Desired Alpha Level*SAL 

1340 IF ALCO OR AU1 THEN PRINT**** ERROR **• ENTER A VALUE BETWEEN 0 ft l*SGOTO 1330 
1330 PRINT 1 PRINT’For Alpha -‘SALS 

1360 IF ALC(2*QOX) THEN PRINT* accept”! ELSE PRINT* reject*} 

1370 PRINT* the hypothesis that the true mean >*SM2 

1380 PRINT*-vs- the alternative hypothesis that true mean -/-*}MZ 

1390 GOTO 1450 

1400 OK 10 GOSUB 1410, 1590, 1750S RETURN 
1410 TZ-( SQRtIN) • (XB - MZ) >/SQR(Vl) 

1420 PRINT S PRINT'Test Statistic - *}TZ 
1430 TX-ABS1TZ) S DF-IN-1 S12-1 1 GOSUB 15000 
1440 GOTO 1320 

1430 PRINT Q960,"Press SPACE BAR when ready to return to main menu*}' 

1440 IF INKET* O * * GOTO 1460 
1470 RETURN 


HULL HYPOT H E SI S 2 

1480 REM •»* NULL HYPOTHESIS 2 »•* 

1490 TZ- (SQRtIN) • (XB - MZ))/SQR(AV> 

1300 PRINT} PRINTTest Statistic -*}TZ 

1310 HX-TZ I VAR-1 t HU-0 S12-1 S GOSUB 10000 

1320 FHMT*F-Value »*}QOX 

1530 PRINTIINPUT’Enter desired Alpha Level*}AL. 

1340 IF ALCO OR AL>1 THEN PRINT'*** ERROR *** ENTER A VALUE BE T WEEN 0 6 f'.GOTO 1530 
1530 FRINTJPRINT'For Alpha -*}AU 

1360 IP AL < QOX T HEN PRINT'acrept”! ft.sf PHINTVeiect*} 

1370 PRINT" the hypothesis that the true mean <-*!MZ 5 PRINT *-vs- the alternative hypothesis that the true me 
an is >*!MZ 
1300 GOTO 1450 

1390 TZ-( SQRtIN) • (XB-MZ) )/SQR(Vl) 

1600 PRINT 1 PRINTTest Statistic -*!TZ 
1610 TX-TZ } DF-GN-U ! 12-1 I GOSUB 13000 
1620 GOTO 1520 


HULL HYPOTHESIS 3 

1630REM *•# NULL HYPOTHESIS 3 •*♦ 

1640 TZ-(SQR(IN) * (XB-MZ))/SQR«AV) 

1630 PRINT l PRINT'Test Statistic -*}TZ 

1660 KX-TZ 1 VAR-t S MU-0 112-t ! GOSUB 10000 

1670 PRINT“P-Valu* -*}FOX 

1680 PRINT : lNFUT*Enter desired Alpha Level"!AL 

1690 IF ALCO OR AL>1 THEN PRINT"*** ERROR ♦** ENTER A VALUE BETWEEN 0 ft l'JGOTO 1680 
1700 PRINT 1 PRINT"For Alpha »*}AL| 

1710 IF AL < FOX THEN PRINT'accept*} ELSE PRINT "reject"} 
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1720 PRINT* the hypothesis that the true mean >*;MZ 

1730 PRINT "-vs- the alternative hypothesis that true mean <"!MZ 

1740 GOTO 1450 

1750 TZ-( SQRtlN) ♦ (XB-MZ) VSQR(Vl) 

1740 PRINT [ PRINT’Test Statistic =»"!TZ 
1770 TX=TZ ! DF-IN-1 ! 12*1 J GOSUB 15000 
1780 GOTO 1670 


HYPOTHESES CONCERNING THE VARIANCE 


2000 

2010 

2020 

2030 

2040 

2050 

2060 

2070 

2080 

2090 

2100 

2110 

2120 

2130 

2140 

2150 

2160 

2170 

2180 

2190 

2200 

2210 

2220 

2230 

2240 

2250 

2260 


REM ••• HYPOTHESES CONCERNING VARIANCE *** 
CLS 

PRINT 044,‘HYPOTHESIS! 


NULL 

PRINT* (1) Var * specified value 
PRINT" or Std, Dev. -dito- 
PRINT* 12) Var O specified value 
PRINT" or Std, Dev, -dito- 
PRINT* (3) Var >- specified value 
PRINT" or Std, Dev. -dito- 
JP-896 ! KS-3 ! GOSUB 760 


ALTERNATIVE* 1 PRINT 
Var ■/■ speafied value 
Std. Dev. -dito-" ! PRINT 
Var > specified value 
Std. Dev. -dito-" J PRINT 
Var < specified value 
Std. Dev. -dito-’! PRINT 


CLS I PRINT’Do you want to enter VARIANCE or STD. DEVIATION?" 

JF-64 ! T1*-‘V* i T2»**D" i GOSUB 710! ZV-ZQ 

PRINT"Enter the value the *5! IF ZV**V* OR ZV**v" PRINT-VARIANCE"! ELSE PRINT*STD. DEVIATION*! 
PRINT" is hypothesiied to be"!! INPUT AV 
IF ZV**D" OR ZV-*d" AV-AV*AV 


PRINT 0192,"Is the mean known?" 

JP-236 ! Tl**"Y* ! T2I-"N" ! GOSUB 710 
IF ZQ-*n* OR Z0*“N" GOTO 2190 
INPUT’Enter Mean’IMZ 

IF ZQ--N" OR ZQ-"n" OR 14-1 TZ-<ON-l) * VI)/AV ! GOTO 2210 

TZ-0 ! FOR I«l TO IN 1 TZ-TZ ♦ (Y(I) - MZ)I2 ! NEXT I! TZ-TZ/AV 

PRINT ! PRINT "Test Statistic »*!TZ 

IF ZQ-*N" DF-IN-1 ELSE DF-IN 

CX-TZ 112-1 I GOSUB 12000 

PRINT “P-Value •"! 

ON IO GOSUB 2270, 2410, 2520 
RETURN 


mLHTrgTHES&l 

2270 REM •** NULL HYPOTHESIS 1 **# 

2280 IF 2*FOX <-l PRINT 2*FOX ELSE PRINT 2*GOX 
2290 PRINT ! INPUT’Enter Desired Alpha Level’JAL 

2300 IF AL<0 OR AU1 THEN PRINT"*** ERROR #+* ENTER A VALUE B E T WEEN 0 6 1‘IGOTO 2290 
2310 PRINT ! PRINT"For Alpha -*!AL 
2320 IF 2*FOX > 1 GOTO 2350 

2330 IF AL < 2*FOX PRINT'Acrept"! ELSE PRINT"Reject"5 
2340 GOTO 2360 

2350 IF AL C 2*QOX PRINT"Accept*! ELSE PRINT’Reject*! 

2360 PRINT" the hypothesis that “!! IF ZV-"V" OR ZV-"v" PRINT -VARIANCE"! ELSE PRINT "STD. DEVIATION"! 
2370 PRINT * -"!!IF ZV-*V" OR ZV*"v" PRINT AV ELSE PRINT AVI.5 

2380 PRINT"-vs- the alternative that"!! IF ZV«“V* OR ZV-*v* PRINT-VARIANCE*! ELSE PRINT "STD. DEVIATIO 
N"! 

2390 PRINT" -/-"!UF ZV-*V" OR ZV»"v* PRINT AV PRINT AVI J 
2400 GOTO 1450 
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MULL HYPOTHESIS 2 


2410 REM • ** NULL HTPQTHESIS 2 *** 

2420 PRINT 00 X 

2430 PRINT : INPUT"Enter Desired Alpha Level'JAL 

2440 IF ALCO OR AOl THEN PRINT"*** ERROR ENTER A VALUE BETWEEN 0 & I'JGOTO 2430 

2450 PRINT ! PRINT “For Alpha =“}AL 

2460 IF ALCOOX PRINT "Accept"} ELSE PRINT "Reject"; 

2470 PRINT" the hypothesis that the "Si IF ZV=“V“ OR ZV-“v" PRINT'VARIANCE"} ELSE PRINT*STD. DEVIATION 

"I 

2480 PRINT* C-*{! IF 2V«"V" OR ZV»"v* PRINT AV ELSE PRINT AVt.5 

2490 PRINT"-vs- the alternative that the *}! IF ZV-"V* OR ZV»“v* PRINT'VARIANCE"! ELSE PRINT “STD. DEVIA 
TION"} 

2500 PRINT" >"!! IF ZV-"V" OR ZV-"v" PRINT AV ELSE PRINT AVt A 
2510 GOTO 1450 


HVU. HYPOTHEC ? 

2520 REM ••• HULL HYPOTHESIS 3 
2530 PRINT FOX 

2540 PRINT i INPUT’Enter Desired Alpha Level'SAL 

2550 IF AL<0 OR AOl THEN PRINT**** ERROR *** ENTER A VALUE BETWEEN 0 & 1‘iGOTO 2540 

2560 PRINT t PRINT"For Alpha *")AL 

2570 IF AL<FOX PRINT "Accept"! ELSE PRINT "Reject"! 

2580 PRINT" the hypothesis that the "}! IF ZV-"V" OR ZV-"v" PRINT'VARIANCE"! ELSE PRINT “STD. DEVIATION 
"I 

2590 PRINT" >"Jl IF ZV-"V" OR ZV-“v" PRINT AV ELSE PRINT AVC.5 

2A00 PRINT"-vs- the alternative that the "}! IF ZV-“V" OR ZV-"v" PRINT "VARIANCE"! ELSE PRINTSTD. DEVIA 
TION"! 

2A10 PRINT * <"!! IF ZV-"V" OR 2V»“v“ PRINT AV ELSE PRINT AVt.5 
2420 GOTO 1450 




ICE IKTSRm FQR THE MEAN 


3000REM ••• CONFIDENCE INTERVAL FOR THE MEAN ••• 

3010 CLS 

3020 PRINT*Is Variance Known?" 1 JP-64 ! T1*»"Y" 5 T2*-*N" | GOSUB 710 
3030 IF ZQ-"N* OR ZQ-"n" GOTO 3190 

3040 PRINT 80. BL* ! PRINT BL» ! PRINT 80,"Ho you wont to enter VARIANCE or STD. DEVIATION?* ! T1»-*V" 

1 T2«*"D" ! JF-44 I GOSUB 7101 ZV-ZQ 

3050 PRINT Enter *ii IF ZV-"V" OR ZV-"v" PRINT'VARIANCE"! ELSE PRINT "STD. DEVIATION"! 

3060 INPUT AV 

3070 IF ZV-«d" OR ZV-*D" AV»AV*AV 

3080 PRINT i INPUT "Enter Desired Confidence Level (e.g. .95HAL ! IF AOl THEN AL-AL/100 
3090 IP ALCO OR AOl THEN PRINT"*** ERROR *** ENTER A VALUE BETWEEN 0 & 1‘iGOTO 3080 
3100' M0*0 ! VR-t } PR-(i+AL)/2 ! 12-1 ! GOSUB 17000 
3110 CLS I PRINT 884. AL*100!"% CONFIDENCE LIMITS" 

3120 PRINT ! PRINT "TWO SIDED CONFIDENCE LIMITS!" 

3130 PRINT" "{SPRINT USING"#***#.***“!XB-IXQ • SQRlAV/INWiFRINT* < True Mean <"}!PRINT USING"* 

M**.****;XB+(XQ • SORIAV/IN1) 

3140 PR-AL ! GOSUB 17000 

3150 PRINT S FRINTONE SIDED CONFIDENCE LIMITS!" 

3160 PRINT Upper -"{SPRINT USING*«****.***"!XB * IXQ • SQRtAV/INl) 

3170 PRINT Lower -“’.SPRINT USING*#****.***“}XB - (XQ * SGR1AV/IN1) 

3180 GOTO 1450 

3190 PRINT ! INPUT "Enter Desired Confidence Level (e.g. .95)"!AL ! IF AL>1 THEN AL-AL/100 
3200 IF ALCO OR AOl THEN PRINT"*** ERROR *** ENTER A VALUE BETWEEN 0 & l".GOTO 3190 
3210 DF-GN-l) S PR-U+AL1/2 112-1 S GOSUB 19000 
3220 CLS i PRINT 884, AL*100!"% CONFIDENCE LIMITS' 

3230 PRINT ! PRINT “TWO SIDED CONFIDENCE LIMITS!" 

3240 PRINT "{SPRINT USING‘M*«*.**«“;XB-<XQ « SQR(V/aN-l)»}5PHINT" C True Mean <"!!PRINT USING" 
••«**.***"!XB*(Xa • SQR(V/(IN-im 
3250 PR-AL ! GOSUB 19000 
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3260 PRINT ! PRINT-ONE SIDED CONFIDENCE LIMITS!" 

3270 PRINT* Upper »*!!PRINT USING**##**.#*#*!XB ♦ IXQ • SQR(V/(IN-1»> 
3200 PRINT* Lower -“JiFRINT USING“*#***.#***JXB - (XQ # SQR(V/ON-l)» 
3290 GOTO 1450 


CONFIDENCE INTERVAL FOR THE VARIANCE 

3300 REM • * » CONFIDENCE INTERVAL FOR VARIANCE ••• 

3310 CLS 

3320 FRINT*Is the Mean Known?* 1 JP-64 ! T1S-*Y* l T2*«*N" 1 GOSUB 710 
3330 IF 2G»*N* OR Z0»“n“ GOTO 3540 
3340 INFUT'Enter Mean'IXB 

3350 INPUT*Enter Desired Confidence Level (e.g. .95)*!AL 1 IF A1D1 THEN AL-AL/100 

3360 IF ALCO OR ADI THEN PRINT* **♦ ERROR •« ENTER A VALUE BETWEEN 0 ft l’SGOTO 3350 

3370 IF 14-1 W2-IN*V ! GOTO 3390 

3380 W2-0 ! FOR 1-1 TO IN ! W2-W2 ♦ ( T(I) - XB )C2 1 NEXT I 
3390 CLS 1 PRINT Q84, AL»100{‘% CONFIDENCE LIMITS* 

3400 PRINT 0320,"TWO SIDED CONFIDENCE LIMITS!* 

3410 DF-IN ! PR-(I+AL)/2 ! 12-1 ! GOSUB 18000 

3420 PRINT0384,* * JIPRINT USING******.»**"JW2/XG!!PHINT" C VARIANCE <" 

3430 PRINT0448** *JiPRIMT USING"#*##e,###“5(W2/XQ)C JJiFRINT* < STD. DEVIATION <“ 

3440 PR-(l-AL>/2 ! GOSUB 18000 

3450 PRINT 0416,**t!PRINT USING******.***’JW2/XQ1 
3460 PRINT 0483,**!!FRINT USING’*****,***“}(W2/XQ>C.5 
3470 PRINT IPRINT’ONE SIDED CONFIDENCE LIMITS!* 

3480 PRINT* Variance Std. Dev.* 

3490 PR-l-AL ! GOSUB 18000 

3500 PRINT * Upper -‘JSPRINT USING******.****JW2/XQ}S PRINT USING*************.***“i(W2/Xa)I.5 
3510 PR-AL ! GOSUB 18000 

3520 PRINT * Lower -‘JiPRINT USING******.***’;W2/XQ|! PRINT USING*************.****SCW2/XQ)l,5 
3530 GOTO 1450 

3540 PRINT I INPUT*Enter deeired confidence level <e.g. .951*!AL ! IF ADI THEN AL-AL/100 
3550 IF AL<0 OR AOl THEN PfelNT"*** ERROR *« ENTER A VALUE B E T WEEN 0 ft 1'iGOTO 3540 
3560 CLS ! PRINT 084, AL*100J*% CONFIDENCE LIMITS’ 

3570 PRINT 0320,“TWO SIDED CONFIDENCE LIMITS!* 

3580 DF-(IN-l) 1 PR-(l+AL)/2 112-1 ! GOSUB 18000 

3590 PRINT0384,* *i!PRINT USING* *#*##. t*e“;(IN*V)/XQ!!PRINT* < VARIANCE C* 

3600 PHINT0448,’ *|!PRINT USING******,*#**J(aN*V)/XQ)I^JiFRINT“ < STD. DEVIATION <*! 

3610 PR-U-AL1/2 i GOSUB 18000 

3620 PRINT Q416,**;:PRINT USING******,****}(IN*V)/XQ} 

3630 PRINT 0483,**J!PRINT USING******. ****i<ON*V)/XQ)t J 
3640 PRINTiPRINr ONE SIDED CONFIDENCE LIMITS!* 

3650 PRINT* Variance Std. Dev.* 

3660 PR-l-AL t GOSUB 18000 

3670 PRINT- Upper ■*J!FHINT USING******.***‘;aN*V)/Xaj!PRINT USING*************.****}<(IN*V)/X 
QKA 

3680 PR-AL ! GOSUB 18000 

3690 PRINT* Lower -*J!PRINT USING"**•♦*.****!ON*VJ/XQIURINT US3NG’*********#**.***"J(ON*V)/X 
QIC a 

3700 GOTO 1450 


MEHV OP TESTS FAR EXTONENTIAl. DATA 

4000REM *•* MENU OF TESTS FOR EXPONENTIAL DATA *«• 

4010 CLS ! PRINT ’DATA IS ASSUMED TO BE EXPONENTIAL 
4020 PRINT i PHINT'MENU!* 

4030 PRINT* (11 Test Hypothesis Concerning Mean (or l/mean) 

4040 PRINT* (2) Confidence Interval for Mean (or 1/me an) 

4050 PRINT S PRINT*NOTE! bath of the above options allow either the Mean or the* 

4060 PRINT* exponential parameter Lambda (1/mean) to be used as the* I PRINT * argument 
4070 JP-896 ! KS-2 ! GOSUB 760 
4080 ON 10 GOSUB 4100, 5000 
4090 RETURN 
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HYPOTHESES CONCERNING THE MEAN 


4100 REM * * * HYPOTHESES CONCERNING MEAN * * • 

4110 CLS ! PRINT TAB(28)*MENUt* 

4120 PRINT 9128,“Do you want to use tha Mean or Lambda (1/mean) as the argument?* 

4130 ,TP*1»2 i T1*»*M* 1 T2*»*L* ! GOSUB 710! ZA*ZQ 

4140 IF ZA-*1* OR ZA»“L" THEN A*»*lambda* ELSE A»«“mean* 

4150 PRINT 9128, CHRX233)! PRINT CHRK2SS) 

4HO PRINT Q64,"HYPOTHESIS!* 

4170 PRINT 0128,* NULL ALTERNATIVE* 

4180 PRINT"! 1) "JASi* > specified value*!: PRINT TAB(33) AS)* «/• specified value* 
4190 PRINT*<2) *!A*i* O specified value*!! PRINT TABI3S) AS]* > specified value* 
4200 PRINT"(3) *5A»!" >» specified value*!! PRINT TAB(3S> AS!" C specified value* 
4210 JP-896 ! KS*3 ! GOSUB 760 

4220 CLS ! PRINT "Enter the hypothesized value of *;A«{!INPUT MZ 
4230 IF ZA»*1* OR ZA»*L* THEN CX»2*IN*XB*MZ ELSE CX-2*IN#X9/MZ 
4240 DF»2*IN 

4250 PRINT i PRINT 'Test Statistic »‘!CX 
4260 GOSUB 12000 
4270 ON IO GOSUB 4290, 4410, 4530 
4280 RETURN 


NULL HYPOTHESIS 1 

4290 REM *** NULL HYPOTHESIS 1 *** 

4300 IF 2*FOX <- 1 THEN PV-2*FOX ELSE PV«2*QOX 
4310 PRINT’P-Value -*!PV 

4320 PRINT ! INPUT'Enter Desired Alpha Level'IAL 

4330 IF AL<0 OR AL>1 THEN PRINT**** ERROR ••• ENTER A VALUE BETWEEN 0 & l*!GOTO4320 

4340 PRINT ! PRINT’For Alpha -*!AL 

4350 IF ALCPV PRINT*Accept**, ELSE PRINT "Reject*! 

4360 PRINT* the hypothesis that the true *!A*!* -*!MZ 

4370 PRINT*-vs- the alternative that the true *!A»I* -/«*!MZ 

4380 PRINT 9960,'-PRESS SPACE BAR TO RETURN TO MAIN MENU*! 

4390 IF MEET* O * * GOTO 4390 
4400 RETURN 

KBUHTTOIHRISZ 

4410REM ••* HULL HYPOTHESIS 2 «*♦ 

4420 PRIHT’P-Value •*!! IF ZA«*1* OR ZA«*L* PRINT FOX ELSE PRINT QOX 
4430 PRINT ! INPUT'Enter Desired Alpha Level*!AL 

4440 IF ALC0 OR AL>1 THEN PRINT**** ERROR ♦*♦ ENTER A VALUE BETWEEN 0 & 1NGOTO 4430 

4450 PRIHT-For Alpha =**!AL 

4460 IF 1A-*1* OR ZA«*L* GOTO 4490 

4470 IF AL<QOX PRINT * Accept*! ELSE PRINT 'Reject"! 

4480 GOTO 4500 

4490 IF ALCFOX PRINT ‘Accept*? ELSE PRINT 'Reject*J 
4500 PRINT * tha hypothesis that the true *?A»?" <-*!MI 
4510 PRINT “-vs- the alternative that the true *IA«!* >*!MZ 
4520 GOTO 4380 


HULL HYPOTHESIS 3 

4530 REM NULL HYPOTHESIS 3 **# 

4540 PRINT"P-Value ■*!! IF ZA»*1‘ OR ZA«*L* PRINT QOX ELSE PRINT FOX 
4550 PRINT t INPUT'Enter Desired Alpha Level’IAL 

4560 I? AL<0 OR AL31 THEN PRINT**** ERROR h» ENTER A VALUE BETWEEN 0 £ l*!GOTO 4550 
4570 PRINT "For Alpha -*!AL 
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4580 IF IA**1" OR IA«“L" GOTO 4610 

4590 IF AL<FOX PRINT "Adapt*! ELSE PRINT " Reject* i 

4600 GOTO 4620 

4610 IF ALIQOX PRINT "Adapt"! ELSE PRINT "Reject”! 

4620 PRINT " the hypothes that the true ‘!A>;‘ >‘!MZ 
4630 PRINT " -vs- the alternative that the true ”!A»!" <"!M2 
4640 GOTO 4380 


CONFIDENCE IN TERVAL FOR THE EXPONENTIAL MEAN 

5000 REM ••• CONFIDENCE LEVEL FOR EXPONENTIAL MEAN 
5010 CLS ! PRINT'CONFIDENCE INTERVAL FOR EXPONENTIAL MEAN* 

5020 PRINT ! INPUT-Entur Desired Confidence Level ( e.g. .95 HAL ! IF AL>1 THEN AL-AL/100 
5030 IF AL<0 OR AOl THEN PRINT"*** ERROR •** ENTER A VALUE BETWEEN 0 6 l“!GOTO 5020 
5040 PR-U-AL1/2 ! DF>2*IN ! 12*1 ! GOSUB 18000 
5050 L*X0/(2*IN*XB) 

5060 PR*<l*AL)/2 i GOSUB 18000 
5070 U*XQ/(2*IN*XB) 

5080 CLS ! PRINT 084, AL*100|"\ CONFIDENCE LIMITS" 

5090 PRINT I FRINT"TWO SIDED CONFIDENCE LIMITS!" 

5100 PRINT USING"###»#*.**#“!1/U!!PRINT" < True Mean <*!!PRINT USING"####*#.*##"!1/L 
5110 PRINT USING"*##**#.***“!L!!PRINT" < Lambda <"!!PRINT USING“**###*.**#*!U 
5120 PRINT 1 PRINT-ONE SIDED CONFIDENCE LIMITS!" 

5130 PR»AL I GOSUB 18000! U-XQ/(2*IN*XB) 

5140 PR-l-AL I GOSUB 180001 L*XQ/(2*IN*XS> 

5150 PRINT "MEANS Upper **!!/L 
5160 PRINT " Lower *“il/U 

5170 PRINT "LAMBDA! Upper **|U 
5180 PRINT " Lower *"!L 

5190 GOTO 4380 


6000 REM *•• MENU OP TESTS FOR DATA A38UMED TO BE BERNOULLI 
6010 CLS 1 PRINT-DATA IS ASSUMED TO BE BERNOULLI 
6020 PRINT 1 PRINT-MENU!" 

6030 PRINT" U) Teet Hypothesis Concerning Proportion <P) 

6040 PRINT" G> Confidence Interval for P 
6090 UP-448 1 KS*2 I GOSUB 760 
6060 ON 10 GOSUB 6100, 7000 
6070 RETURN 




.CatiC.EMIKQ.THS PROTORUBM 


CONCERNING PROPORTION «•* 
MULL ALTERNATIVE- ! PRINT 

P -/■ specified value 
P > specified value 
P < specified value 


6100 REM BTPOTHESES 

6110 CLS ! PRINT 964,‘HYPOTHESIS! 

6120 PRINT" (1) P* specified value 
6130 PRINT” (2) P O specified value 
6140 PRINT* (3) P >* specified value 
6150 UP-640 ! XS*3 ! GOSUB 760 
6160 CLS 1 INPUT*Enter Hypothesized Value of P"!FB 

6170 IF FB<0 OR PBM THEN PRINT"*** ERROR*** P MUST BE A VALUE B E T WEEN 0 AND 1" ! FOR 1*1 TO 
2000 1 NEXT ! GOTO 6160 

6180 IF IN > 20 THEN YS=(IN * XB - IN * FBI / SQRUN * PB • il-PB)) ELSE YS-IN * XB 

6190 PRINT 1 PRINT'Test Statistic *"!YS5! IF IN < 20 PRINT" (Using Binomial CDF)" ELSE PRINT" (Using N 

ormal CDF)" 

6200 ON 10 GOSUB 6220, 6410, 6550 
6210 RETURN 
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MULL HYPOTHESIS 1 

6220 REM **• NULL HYPOTHESIS 1 *** 

6230 IF IN >« 20 GOTO 6300 

6240 N-IN : 12*1 : K-YS ! GOSUB 11000 

6250 PV»2*PR 

6260 IF PV C 1 GOTO 6290 

6270 'K-YS-l ! GOSUB 12000 

6280 FV-2*(1-PR) 

6290 GOTO 6320 

6300 MU-0 1 VAR-1 ! NX-ABS<YS> ! 12-1 $ GOSUB 10000 
6310 PV»2«QOX 

6320 PRINT i PRINT" P-Value -*;PV 

6330 PRINT i INPUT"Enter Desired Alpha Level" !AL 

6340 IF ALCO OR ADI THEN PRINT"*** ERROR ••« ENTER A VALUE BETWEEN Oil"! GOTO 6330 
6350 IF AL C PV PRINT "Accept"! ELSE PRINT "Reject"! 

6360 PRINT" the hypothesis that the true value of P »*!PB 

6370 PRINT"-vs- the alternative that the true value of P »/»“}FH 

6380 PRINT ! PRINT"*** PRESS SPACE BAR TO RETURN TO MAIN MENU* 

6390 IF INKEY* O " " GOTO 6390 
6400 RETURN 


NULL HYPOTHESIS 2 

6410 REM ••• HULL HYPOTHESIS 2 ••• 

6420 IF IN >- 20 GOTO 6460 

6430 N-IN 1 K-YS-l ! 12-1 S GOSUB 11000 

6440 PV-l-PR 

6450 GOTO 6480 

6460 MU-0 ! VAR-1 ! NX-YS !12-1 ! GOSUB 10000 
6470 PV-QOX 

6480 PRINT 1 PRINT"F-Value -*!PV 

6490 PRINT 1 INPUT "Enter Desired Alpha Level’SAL 

6500 IF ALCO OR AL>1 THEN PRINT"*** ERROR *** ENTER A VALUE BETWEEN 0 ft I’iGQTO 6490 
6510 IF AL < PV PRINT "Accept"! ELSE PRINT "Reject"! 

6520 PRINT" the hypothesis that the true value of P <-"5FB 
6530 PRINT“-vs- the alternative that the true value of P >*JFB 
6540 GOTO 6380 


NULL HYPOTHESIS 3 

6550 REM ••* NULL HYPOTHESIS 3 **» 

6560 IF IN >- 20 GOTO 6600 
6570 N-IN 1E-YS i 12-1 i GOSUB 11000 
6380 PV-PR 
6590 GOTO 6620 

6600 MU-0 t VAR-1 ! NX-YS ! 12-1 1 GOSUB 10000 
6610 PV-FOX 

6620 PRINT ! PRINT"P-Value -*!FV 

6630 PRINT ! INPUT"Enter Desired Alpha Level'lAL 

6640 IF ALCO OR ADI THEN PRINT"*** ERROR *** ENTER A VALUE BETWEEN 0 8 1"! GOTO 6630 
6650 IF AL < PV PRINT "Accept"! ELSE PRINT "Reject"! 

6660 PRINT" The hypothesis that the true value of P >-"!FB 
6670 PRINT"-vs- the alternative that the true value of P <*!PB 
6680 GOTO 6380 
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CONFIDENCE TN TERVAL FOR PROPORTION 


7000 HEM *** CONFIDENCE INTERVAL FOR PROPORTION • ** 

7010 CLS 

7020 INPUT'Enter the desired confidence level (e.g. .9S)"!AL 1 IF ALM AL*AL/100 

7030 IF AL<0 OR AL>1 THEN PRINT"#*# ERROR *** ENTER A VALUE BETWEEN 0 & 1‘JGOTO 7020 

7040 PR*U*AL)/2 

7050 DN*2#IN - 2*IN*XB ♦ 2 ! DD»2*IN*XB ! 12*1 1 GOSUB 20000 
7060 CLS J PRINT 034, AL#I00S“% CONFIDENCE LIMITS" 

7070 PRINT ! PRINT'TWO SIDED CONFIDENCE LIMITS!" 

7080 L-UN • XB) / «IN*XB> * (IN - IN#XB * 1) • XQ) 

7090 PRINT" "!L}"< True Value of P <"S 

7100 DN*2#IN#XB * 2 ! DD*2*IN - 2*IN#XB ! GOSUB 20000 

7110 U*((IN # XB +1) * XQ) / ((IN - IN*XB> ♦ UN » XB * 1) • XQ) 

7120 PRINT U 

7130 PRINT ! PRINT'ONE SIDED CONFIDENCE LIMITS!" 

7140 PR*AL ! GOSUB 20000 

7150 PRINT" Upper *"!«IN # XB + 1) • XQ) / (UN - IN*XB> ♦ UN#XB ♦ 1) * XQ) 

7160 DN*2#IN - 2*IN*XB ! DD=2#IN»XB ! GOSUB 20000 

7170 PRINT" Lower -"iUN • XB) / UN # XB + UN - IN#XB ♦ 1) * XQ) 

7180 GOTO 6380 


NORMAL DISTRIBUTION 

10000 REM ••• NORMAL DISTRIBUTION #•* 

10010 SD-SQRCVAR) 

10020 XN*(NX-MU)/SD!AX-ABS(XN> 

10030 XI».3989423#EXP((-AX*AX>/2) 

10040 NT*1/(1+.2316419#AX) 

10030 OOX»XZ#((.3193815«NT>-{.3S63638*NTt2W1.7Bi478*NT£3>-<1.821236tNTt4)+(1.330274»NT[3)> 
10060 IF XN<0 GOTO 10090 . 

10070 FOX*lrQOX 
10060 RETURN 
10090 FOX-QOX 
10100 QOX*l-FQX 
10110 RETURN 


amotm aaBagnoa 

11000 REM ••• BINOMIAL DISTRIBUTION •«* 
11010 PR*0 

11020 IF K<0 GOTO 11080 

11030 IF K > N THEN PR-1 ! GOTO 11080 

11040 FOR BB-0 TO K 

11050 TB-N ! GOSUB 11090 

11060 PR*PR+ BC * PBIBB # (l-PB)t(N-BB) 

11070 NEXT BB 
11080 RETURN 

11090 BC*0 I B-0 ! BT*0 I SB*BB S QB»TB-BB 

11100 FOR BI*1 TO TB 

11110 IF SB>«QB GOTO 11170 

11120 XF SB >* 81 THEM B-B * LOG(BI) 

11130 IF QB > BI GOTO 11210 

11140 IF QB - BI THEN BC—B ! GOTO 11210 

11150 BC*BC * LOGCBI) 

11160 GOTO 11210 

11170 IF QB > BI THEM B*B * LOG(BI) 

11180 IF SB > BI THEN GOTO 11210 
11190 IF SB - BI THEN BC—S I GOTO 11210 
11200 BC-BC * LOGIBD 
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U210 NEXT BI 
11220 BC=EXP(BC) 
11230 RETURN 


CHIiSQy m PIStfflB UTW N 

12000 REM CHI-SQUARE DISTRIBUTION 

12010 IT DF>30 GOTO 12220 
12020 DP-1 

12030 FOR I»INT(DF> TO 2 STEP -2 
12040 DP»DP*I 
12050 NEXT I 

12060 NP-C X C dNTI1DF+1 1/21MEXF1-CX/21/DP 

12070 IF INT(DF/2)°DF/2 GOTO 12100 

12080 PC-SQR(2/C X/3.1415931 

12090 GOTO 12110 

12100 PC-1 

12110 LC-1 

12120 MC-1 

12130 D-DF 

12140 D-D+2 

12150 MC-MCtCX/D 

12160 IF MC<.000001 GOTO 12190 

12170 LC-LC+MC 

12180 GOTO 12140 

12190 FOX-PC*NP*LC 

12200 QOX-l-FOX 

12210 RETURN 

12220 MX-UCX/DF1C(1/3) - U-(2/(9»DFm>/SQRC2/(9«DF)> 

12230 MU-0 1 VAR-1 112-1 
12240 GOSUB 10000 
12250 RETURN 


t, ^DISTRIBUTION 

15000 REM **• t-DISTRIBUTION 
15010 IF TX<013-1 ELSE 13-0 
15020 IF DF>25 GOTO 15450 
15030 THETA-ATN(ABS(TX)/SQR(DF» 

15040 IF DFOl GOTO 15090 
15050 AT-»6366198*TH£TA 
15060 FOX-1AT+1I/2 
15070 QOX-l-FOX 
15080 GOTO 15400 

15090 IF IHTCDF/2HDF/2 GOTO 15220 

15100 IF DF-3 GOTO 15370 

15110 TS-0 i TP-1 ! TC -C OSCTHET A) 

15120 FOR 1-1 TO INTHDF-31/2) 

15130 TP-TP*2*I/(l+2*I) 

15140 TC«TC*<COS(THETA)£2> 

15150 TS-TS+(TP*TC1 

15160 NEXT I 

15170 TS-TS+COSCTHETA) 

15180 AT-.6366198*{THETA*ISIN(THET AMT51) 

15190 FQX-lAT+11/2 

15200 QOX-l-FOX 

15210 GOTO 15400 

15220 IF 0FO2 GOTO 13270 

15230 AT-SIN(THETA) 

15240 FOX-<AT+l)/2 
15250 QOX-l-FOX 
15260 GOTO 13400 
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15270 TS*0 : TP*1 : TC*1 
15280 FOR 1*1 TO INT((DF-2)/2) 

15290 TP*TF»((2*1)-1)/<2*I) 

15300 TC*TC*(CQS(THETA)C2) 

15310 TS=TS+(TP*TCi 
15320 NEXT I 

15330 AT a SINfTHETA)*<l+TS) 

15340 FQX*<AT+ll/2 
15350 OOX»1-FOX 
15360 GOTO 15400 

15370 AT*.6366198«<THETA*<SIN<THETA>«COSCTHETA») 

15380 FOX*(AT+l)/2 

15390 GOX*l-FOX 

15400 PX-TX 

15410 IF 13*0 RETURN 

15420 QOX*FOX 

15430 FOX-l-QOX 

15440 RETURN 

15450 NX-I TX • U - ,25/DF) 1/SQR(1 ♦ ,5*TXC2/DF) 
15460 MU*0 ! VAR*1 112*1 1 GOSUB 10000 
15470 PX*TX 
15480 RETURN 


r - 

16000 REM F-DISTRIBUTION *** 

16010 IF FX*0 THEN FOX*0 t QOX*l ! RETURN 

16020 IF DIO40 AND DD740 GOTO 16610 

16030 Xl-DN ! K2-OD 

16040 A»DN*FX/(DN*FX+DD> !,A1*1-A 

16050 IF AK1E-37 THEN A1-1E-37 

16060 Dl-DNO ! D2*DD*J 1 D3-D1+D2-1 

16070 R*0 1 S1*0 ! S2*0 5 DL«U C*.2S 

16080 XM-1 1 XX* 1 t PI-3.141593 ! KN*DD 

16090 KM«INTCD2>*2 

16100 IF KMOKN GOTO 16280 

16110 XN-D2-1 

16120 IF KN*0 GOTO 16210 

16130 FOR 1*1 TO KH 

16140 Sl-DL+SUR 

16150 D2-D2-1 

16160 D3-D3-1 

16170 TEM-A1/D2 

16180 R*D3*TEM 

16190 S2*(R+TEM)«S2 

16200 NEXT I 

16210 S1*DL+S1*R 

16220 DL*0 

16230 T—l 

16240 D3»-l 

16250 S2*A*S2 

16260 C-C+.5 

16270 GOTO 16460 

16280 KN*INT(D2) 

16290 IF KN*0 GOTO 16380 
16300 FOR 1*1 TO KN 
16310 Sl-DL+SUR 
16320 D2-D2-1 
16330 D3*D3-1 
16340 TEM-A1/D2 
16350 R*D3»TEM 


91 



16360 S2=(R+TEM)*S2 
16370 NEXT I 
16380 S1»XK*S1 
16390 S2-XK«S2 
16400 ART-SQR(Al) 

16410 XM-XM»ART 
16420 T»(XM-ART)/Al 
16430 D3—.3 
16440 XK-2/FI 
16450 r«C #2 

16460 IF 0.873 GOTO 16330 

16470 D2-D1 

16480 D3-D2+D3 

16490 52=51 

16300 St-0 

16310 Al-A 

16320 IF AK1E-37 THEN A1-1E-37 
16330 KN-DN 
16540 GOTO 16090 

16330 IF C<lil23 THEN DL-4/PI * ATN(T> 

16360 FOX-XM*(S2-S 1HDL 
16570 IF FOXCO THEN FOX-O 
16580 IF FOXOl THEN FOX-1 
16390 QOX-l-FOX 
16600 R E T U RN 

16610 NX-(((FXC(l/3»*(l-(2/(9*DD)})>-{l-(2/(9*DN))))/SQR({2/(9*DN))+(FXt(2/3»*2/(9*DD)) 
16620 MU-0 ! VAR-1 
16630 GOSUB 10000 
16640 RETURN 


INVERSE NORMAL 

17000REM *'*« INVERSE NORMAL DISTRIBUTION *•* 

17010 IF TtOa THEN RL-SQRO.QG(l/«l-PR>C2)» ELSE RL-5QR(LOG( 1/CPRC2))) 

17020 XQ-RL-<(2J13517+(^02833#RL)+(.010328#RLt2))/(l+(1.432788»RL>+<.189269*RLt2W.001308*RLt3m 
17030 IF FR>-J THEN XQ-MO(SQR(VR>*XQ) ELSE XQ-MO-(SOH(VR)*XQ> 

17040 RETURN 


INVERSE CHI- SQUARE 

18000 REM #♦* INVERSE CHI-SQUARE DISTRIBUTION • •• 
18010 MO-O t VR-1 112-1 ? PZ-FR 
18020 IP DFM GOTO 18070 
18030 PR-«5*(1-FZ) 

18040 GOSUB 17000 
18050 XQ-XQI2 I FR-FI 
18060 RETURN 
18070 IF DF>2 GOTO 18100 
18080 XQ—2tLQG(l-FW 
18090 RETURN 
18100 GOSUB 17000 
18110 DQ-XQ 

18120 IF DF>(2+IKT(4tABS(DQ))) GOTO 18220 

18130 X9-1.0000886 - .2237368/DF - .01513904/DFI2 

18140 X9-X9*(DFt- J • DQ • (.4713941 ♦ ,02607083/DF - .008986007/DFC2)) 

18130 X9-X9 + (1/DF * DQI2 * (.0001348028 + .01128184/DF + .02277679/DFC2J1 
18160 X9-X9 ♦ (DFI-1.3 * DQC3 * (-.008533069 - .01153761/DF - ,01323293/DFt2» 
18170 X9-X9 ♦ (DFC-2 * DQ£4 * (.00312558 ♦ .003169654/DF - .006930336/DF12» 
18180 X9-X9 + (DFI-2J * DOC3 * (-.0008426812 ♦ 4W253001/DF ♦ .001060438/DFI2)) 
18190 X9-X9 ♦ (DFt-3 • DQI6 # (.00009780499 - ,001450117/DF ♦ ,001565326/DFC2» 
18200 XQ-DF«X9C3 
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18210 RETURN 

18220 X9-1 - 2/(9«DF) ♦ (4#DQC4 ♦ 16»DQC2 -28)/(1213»DFC2> 

18230 X9-X9 ♦ <8»D0t6 ♦ 720»DQC4 ♦ 3126O0C2 ♦ 2904>/(229635*DFC3) 

18240 XB-DQ/3 ♦ (-DQC3 ♦ 3*DQ)/(162*DF) 

18250 X8-X8 - (3*DQC5 ♦ 40»DQC3 ♦ 43*DQ)/(3832*DF£2) 

18260 X8-X8 ♦ I301-D0C7 - 1517»DQI5 -32769*DQC3 - 79349»DQ)/(7873200*DFC3> 
18270 X8- SQR(2/DF) • X8 
18280 XQ-DF * <(X8+X9)C3> 

18290 RETURN 
18300 CX-XQ 
18310 0ELTA*>01 
18320 C05UB 12000 

18330 IF ABSIFR-FOXX.OOOOl THEN XQ-CX t RETURN 

18340 IF FR < FOX THEN CX-CX-DELTA 

18330 IF PR < FOX THEN DELTA 3 DELTA/10 

18360 IF PR < FOX GOTO 18320 

18370 CX-CX+DELTA 

18380 GOTO 18320 


mVEB§E.--,t 

19000 REM «»» INVERSE t-DISTRIBUTION 
19010 M0*0 t VR-l 112»l 
19020 GOSUB 17000 
19030 XC-XQ 

19040 XCW.23 • (XCC3 ♦ XCH/DF 

19050 XQ-XQ -*• (11/96) • (S*XCC3 + 16*XCC3 ♦ 3*XC»/DFC2 

19060 XQ-XQ ♦ (0/384) • <3»XCC7 ♦ 19*XCI3 ♦ 17«XCC3 -15*XC»/DF£3 

19070 XQ-XQ ♦ (0/92160) * (79»XCC9 ♦ 776-XCI7 ♦ 1482«XCC5 -1920*XCC3 -943*XC))/DFC4 

19080 XQ- XC ♦ XQ 

19090 IF DF>3 RETURN 

19100 DELTA-1 

19110 TX-XQ 

19120 GOSUB 15000 

19130 IF ABS(FH-FOX) < .00001 THEN XQ-TX t RETURN 

19140 IF PR < FOX THEN TX-TX-DELTA 

19150 IF PR < FOX THEN DELTA-DELTA/10 

19160 IF PR < FOX THEN GOTO 19120 

19170 TX-TX+DELTA 

19180 GOTO 19120 


BTCEBSSrF 

20000 REM INVERSE F-DISTRIBUTIOM *« • 

20010 IF DN < 20 OR DO < 20 GOTO 20090 
20020 MO-O 1 VR-l 112-1 
20030 GOSUB 17000 

20040 H - 2 • ((l/GDD-t) ♦ 1/0DK-1»C-1 > 

20050 LM - (XQ12 - 3>/6 

20060 HX-<XQ • SQR(H * LM»/H - ((1/(DN-1) - l/OO-l)) * (LM ♦ 3/6 - 2/<3»H))) 
20070 XQ - EXPC2»HX> 

20080 RETURN 

20090 IF DNOl AND DDOl OR DN-1 AND DD-1 GOTO 20190 
20100 PH-PR 1 IF DN-1 THEN PR-(l+PR)/2 ELSE PR-l-PR/2 
20110 IF DN«1 THEN DF-00 ELSE DF-ON 
20120 GOSUB 19000 

20130 XQ»XQ*XQ 1 IF DD-1 THEN XQ-l/XQ 
20140 PR-PH 

20130 IF PK>-.99 THEN XQ-.994-XQ ! GOTO 20180 

20160 IF FH>-.973 AND DN>3 THEN XQ-.993«XQ 1 GOTO 20180 
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20170 IF PR>-.95 AND DN>3 THEN XQ».996*XQ 
20180 RETURN 

20190 IF FR>.999 AND DD-1 THEN XQ-9E+36 ! RETURN 
20200 IF PR>.997 AND DD-1 THEN DELTA-500000 : GOTO 20270 
20210 IF PR>.995 AND DD-l THEN DELTA-15000 : GOTO 20270 
20220 IF PR>.99 AND DD-1 THEN DELTA-5000 : GOTO 20270 

20230 IF PR>.973 AND DD=1 OR PR>.999 AND DD-2 THEN DELTA-1000 i GOTO 20270 

20240 IF PR>i9 AND DD-1 OR PR>.99 AND DDC-2 OR PRQ-.999 AND DD<-3 THEN DELTA -100 : GOTO 20270 

20230 IF PR>*.99 AND DD<-3 THEN DELTA-50 i GOTO 20270 

20260 IF PR>-.975 AND DOC-2 THEN DELTA-50 ELSE DELTA-10 

20270 LL-0 

20280 FX-DELTA 

20290 GOSUB 16000 

20300 IF ABSIPR-FOXK-.OOOOl THEN XQ-FX ! RETURN 

20310 IF PRCFOX THEN DELTA-DELTA/2 ! FX-DELTA : GOTO 20290 

20320 LL-FX 5 FX-FX+DELTA i UL-FX ! GOSUB 16000 

20330 IF ABS<PR-FOXK-.00001 THEN XQ-FX 1 RETURN 

20340 DELTA-(UL-LL)/2 

20330 FX-LL+DELTA 

20360 GOSUB 16000 

20370 IF ABStPR-FOXK-.OOOOl THEN XQ-FX ! RETURN 
20380 IF PR>FOX THEN LL-FX ELSE UL-FX 
20390 GOTO 20340 


PROGRAM SELECTION MENU 


22000 

22010 

22020 

22030 

22040 

22030 

22060 

22070 

22080 

22090 

22100 

22110 

22120 

22130 

22140 

22150 

22160 

22170 


REM • PROGRAM SELECTION MENU ••• 

CLS 

PRINT TAB(20)*PROGRAM SELECTION MENU* 

PRINT @128«*#1 Data I/O, Transformations, Prab. & Inv. Prob. Distributiona-'FRINT 
PRINT**2 Sample Stat's, 7-Num-Sum, Plots*!PRINT 

PRINT**3 Parametric Tests of Hypotheses and Confidence Intervals*!PRINT 
PRINT"*4 Robust Estimates and Confidence Intervals*SFRINT 
PRINT**3 Multiple Linear Regression'IFRINT 
PRINT**6 Analysts of Variance* 

JP-896 ! KS=6 ! GOSUB 760 

IF 10-3 GOTO 10 

IF 10-4 THEN ZT-*2"'.GOTO 22130 

IF 10-1 OR IO-2 THEN ZT-*1* ELSE IT-"3* 

IF 10-1 OR 10-3 THEN 10-1 ELSE 10-2 
CLS 

PRINT-INSERT TAPE #*TO* INTO ESF-60 AND PRESS SPACE BAR* 

IF INKETSO* * GOTO 22160 
QLOAD 10 


END OF PROGRAM *3 



PROGRAM =*4 


MENU 

10 LBS=CHRS(123> ! RBS=CHRS(125) ! BLS=CHR»(250> 

20 REM *•* M A I N M E N U 
30 IF 14=1 GOTO 130 

40 CLS 5 PRINT TAB (20) “*«* MENU i PRINT 

50 PRINT"(1) Robust-t Confidence Intervals for the Population Mean* 5 PRINT 

A0 PRINT“<2) Confidence Interval for the Median* ! PRINT 

70 PRINT*I3) Jackknife (Point Estimates & Confidence Intervals)* 5 PRINT 

60 PRINT"(4) Return to Program Selection Menu 

90 KP-768 ! XS=4 

100 GOSUB ISO 

110 ON 10 GOSUB 500, 2000, 3000, 10000 
120 GOTO 20 
130 CLS 

140 FRINT*None of the subroutines in Program *4 are applicable since only sample statistics were entered. Pre 
ss the SPACE BAR to access the program selection menu. 

150 IF INKEYIO* * GOTO 150 
160 CLS 

170 GOTO 10000 

ISO PRINT QKP+5,“OPTION DESIRED “lLBSiCHRi(176)iRB»J 
190 ZI-INKEYS 5 IF ZI=*“* GOTCJ 190ELSE IO=VALOI) 

200 PRINT 0KP+21,a i FOR 1=1 TO 100 ! NEXT 

210 IF I0<1 OR IQ>KS THEN PRINT (JKF, BLS 1 PRINT QKP,**** NOT A VALID OPTIOM ***“ 1 FOR 1=1 TO 1000 I 
NEXT i PRINT SKP, BL* S GOTO 180 
220 RETURN 


BOattSI-r-t 

500 REM •** ROBUST-t SUBROUTINES *** 

510 CLS 

520 PRINT TABl22)*ROBUST-t TECHNIQUES" J PRINT 

530 PRINT" The two robust techniques listed below are useful in* 

540 PRINT “calculating confidence intervals for the population mean* 

550 PRINT*when known or suspected outliers violate the assumption of* 
560 PRINT’near normality!* 

570 PRINT 

580 PRINT" (1) Winsorized-t Technique 
590 PRINT" (21 Trimmed-t Technique 
600 KP-896 ! JCS-2 1 GOSUB 180 
610 CLS 

620 IF 10*1 ZR=* Winsorixing * ELSE ZR»"Trimroing* 

630 PRINT”Enter Desired Uvel of “JZR5* (10 max)*;: INPUT JG 
640 IF JG=0 THEN WM-XB ! HS-SQR(Vl) ! GOTO 810 
650 FOR 1=1 TO JG 
660 Wa>=Y(JG+l) 

670 WW(I)=Y(IN-JG) 

680 NEXT I 
690 SM»0 
700 S2«0 

710 FOR 1=1 TO JG 
720 SM-SM+WairWWai 
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730 S2*S2+W<IH2 ♦ WW(I)C2 
740 NEXT I 

750 FOR I-JG+1 TO IN-JG 
760 SM«SM*Y(I> 

770 S2*S2*Y(I)C2 
780 NEXT I 
790 WM*SM/1N 

800 WS*SQR«S2-<SMC2)/IN)/(IN-1» 
810 ON 10 GOSUB 1000, 1500 
820 RETURN 


WINSORIZED - t 

1000 REM • WINSORIZED-t ♦ ** 

1010 CLS 

1020 PRINT TAB(25)‘WINSORIZED - t* 

1030 PRINT 0128,"At Uvtl *-,JG 

1040 PRINT* Winsorixed Point Estimate for the Mean »*}WM 

1050 PRINT 1 INPUT*Enter Desired Confidence Level (e»g< ,93)*SAL 

1060 IF AL>*1 THEN AL-AL/100 

1070 PR-U+ALJ/2 

1080 DF-IN-2*JG-1 

1090 GOSUB 8000 

1100 PRINT 1 PRINT AL*100{*\ CONFIDENCE INTERVALS* 

1110 PRINT S PRINT TABUS) WM - (IN-1)/{DF) ♦ XQ * WS/SORON)}* < True Mean <“{ WM * CIN-l)/tDF) * XQ • H 
S/SQRON) 

1120 PRINT 0940,"PRESS SPACE BAR TO RETURN TO MAIN MENU*} 

1130 IF INKEYt O * " GOTO 1130 
1140 RETURN 


riMIUja*; 


1300 REM •** TRIMMED-t 
1310 CLS 1 PRINT TAB(23)*TRIMMED-t* 

1320 PRINT 0128,*At Level -“SJG 

1330 PRTNT’TruTvned Point Estimate for the Mean »"! 

1340 SM«0 

1530 FOR WG-t-1 TO IN-JG 
1360 SM-SM+Yd) 

1370 NEXT I 

1580 TM-SM/ON-2*JG> 

1390 RW*MtIN-1 )/GN-2» JG-1 »»WSC2) 

1600 PRINT TM 

1610 PRINT S INPUT"Enter Desired Confidence Level <e.g. .93)‘UL 

1620 IF AL)*1 THEN AL-AL/100 

1630 DF»IN-2*JG-1 

1640 PR*<l+AL>/2 

1630 GOSUB 8000 

1660 PRINT S PRINT TABUS) TM - XQ«SQR<RW/(IN-2*JG))I*< True Mean <*STM ♦ XQ*SQHtRH/ttN- 2 tJG» 
1670 PRINT 0960,‘PRESS SPACE BAR TO RETURN TO MAIN MENU*} 

1680 IF INKEYS O * * GOTO 1680 
1690 RETURN 


CONFIDENCE INTERVAL FOR THE MEDIAN 

2000 REM ••• CONFIDENCE INTERVALS FOR THE MEDIAN 
2010 CLS 

2020 PRINT 013,‘CONFIDENCE INTERVAL FOR THE MEDIAN- 

2030 PRINT 0128,* Calculation of a confidence interval for the median is* 

2040 PRINT'based on the binomial distribution which is, of course, a dis-* 
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2090 PRINT'crete distribution. An interval will be calculated such that* 

2060 PRINT'the confidence level you soeafy will be met or exceeded.* 

2070 PRINT ! INPUT”Enter Desired Confidence Level le.g. .95)*IAL 

2080 IF AL>»1 THEN AL=AU100 

2090 CLS : PRINT Q471,'CALCULATING* 

2100 P-.5 
2110 PR-0 

2120 IF IN<=30 BB-0 ! GOTO 2180 
2130 IF IN7-150 BB-55 ! GOTO 2170 
2140 IF IN>100 BB=35 1 GOTO 2170 
2190 IF IN>60 BB-20 ! GOTO 2170 
2160 IF 110-40 THEN BB-12 ELSE BB-S 
2170 MU-.5 ! VAR-.25/IN 1 GOTO 2210 
2180 GOSUB 2900 
2190 PH - PR 

2200 PR-PR + BC * PtIN ! GOTO 2220 

2210 PH-PR 1 NX-(BB+.5)/IN ! GOSUB 9000S PR-FOX 

2220 IF PR < (1 - AL1/2 THEN 88-BB+l ! IF IN>30 GOTO 2210ELSE GOTO 2180 
2230 PR-1-2* PH 
2240 CLS 

2290 IF PR-1 THEN PR-.995 

2260 PRINT 910,""{SPRINT USING*###.**{PR*100{:FHINT-X CONFIDENCE INTERVAL FOR THE MEDIAN* 
2270 IF BB-0 THEN BB-1 

2280 PRINT 9148, Y(BB);*C Median <"{ YUN-BB+l) 

2290 FOR BB-R1 TOSl-l 

2300 PRINT 9960,-PRESS SPACE BAR TO RETURN TO MENU*; 

2310 IF INKEY* O * * GOTO 2310 
2320 RETURN 

2500HEM *** BINOMIAL COEFFICIENT 

2910 BC-0 ! B-0 t SB-BB ! QB-IN-BB 

2920 FOR BI-1 TO IN 

2930 IF SB3-QB GOTO 2990 

2940 IF SB>*BI THEN B-B+LOGiBI) 

2950 IF QB>BI GOTO 2630 

2960 IF QB-BI THEN BC—B ! GOTO 2630 

2570 BC-BC ♦ LOG(BI) 

2980 GOTO 2630 

2990 IF QBJ-BI THEN B-B+LOG(BI) 

2600 IF 5B>BI THEN GOTO 2630 
2610 IF SB-BI THEN BC—B I GOTO 2630 
2620 BC-BC+LOG(BI) 

2630 NEXT BI 
2640 BC - EXPfflC) 

2690 RETURN 


JACKKNIFE ESTIMATES & CO NFIDENCE INTERVALS 

3000 REM **# JACKKNIFE SUBROUTINE **• 

3010 CLS : PRINT TAB (22) -JACKKNIFE* 

3020 PRINT! PRINT*The jackknife technique is useful for reducing bias in the estimates of parameters and i 
n obtaining robust interval estimates. 

3030 PRINT 

3040 FRINT*HOTE 1’ Estimate & Conf. Int. far Variance are based on Jack- knife of log (base 10) of s* 1 PRINT 
3090 PRINT*NOTE 2: Jackknife Est, & Conf, Int. far the Mean are identical to classical calculations and are md 
uded here for convenience only, 

3060 PRINT 

3070 PRINT* *;LB»:*l-mBS5- Mean- 
3080 PRINT* *;lB* 5*2*;RB*;* Variance* 

3090 KP-960 ! KS-2 S GOSUB 180 

3100 CLSIIF 10-1 THEN PE-XB 5IH-1 : SD-SQR(Vl); GOTO 3220 
3120 PM-LOG (SQRCV1 ))/LOG (10) 

3130 PE-0 ! S2-0 
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3140 FOR 1*1 TO IN STEP IH 
3130 GOSUB 3450 

3160 SV*(IN/IH)*PM - (IN/IH-1)*L01 
3170 S2-S2 ♦ SV»SV 
3180 PE-FE + SV 
3190 NEXT I 

3200 SD-SQRl (S2 - (PEC2)/(IN/IH)t/(IN/IH-l) ) 

3210 FE*FE*IH/IN 
3213 CLS 

3220 PRINT 00 iC HRS (253) ! FRINT’Fotnt estimate of *1 
3230 IF IQ»1 ZF*“Mean“ ELSE ZP*‘Variance* 

3240 PRINT ZP5* *“! 

3230 IF 10*1 PRINT PE ELSE PRINT 10KPE+PE) 

3260 IF 10*2 PRINT"Point estimate of Std. Deviation **!10CPE 

3270 PRINT ! INPUT-Enter the desired Confidence Level <e.g. .93)*’AL 

3200 IF AL>1 THEN AL-AL/100 

3290 PR*1 - (1-AL1/2 ! IF PR>1 PR-PR/100 

3300 DF»IN/IH-1 ! 12*1 : GOSUB 8000 

3310 L*PE - XQ * SD/SQR(IN/IH) ! IF 10*2 L»10t(L*L) 

3320 U*PE ♦ XQ * SO/SQRGN/IH) i IF 10*2 U*10C(U*U) 

3330 PRINT ! PRINT AL*tOO;-\ Confidence Interval!" 

3340 PRINT* *!L!* < True Value of "!ZP!" <*!U 

3350 IF 10*2 PRINT" “!SQR(L)i" C True Value of Std. Deviation <"!SQR<U> 

3360 PRINT ! PRINT i PRINT ! PRINT'*** PRESS SPACE BAR TO RETURN TO MAIN MENU* 

3370 IF INKEYS O * * GOTO 3370 

3380 RETURN 

3400 SP*0 ! FOR 1-1 TO 5 

3405 IF IN/I*INT(IN/I) THEN PG(I)-1 ELSE PGID-O 
3410 IF PG(I>»1 THEN SP=SP*1 
3413 NEXT I 

3420 CLS t IF SP>1 THEN 3432 ELSE PRINT * 5ince*;iN5*(the sample size) is a prime number it is only" 

3423 PRINT"possible to jackknife in groups of 1. For large sample sizes > computation time is excessive**! 

3430 FRINTJGOTO 3433 

3432 PRINT* You may jackknife in groups as listed below* Computation time Is a function of both the sample 
size and group size and becomes excessive with large samples and small groups*' 

3433 PRINT:PRINT*e.g. Sample Size Group Size Approx. Time“!PHINT* 100 1 23 min.'TIF 

SFM PRINT* 100 2 0 min.'SPRINT* 50 1 4 min." 

3434 PRINTiPRINT'You may! *!LB*!*Q*!RB*i"uit* 

3433 TL*1'PRINT* or*!FRINT*Jackkmfe in groups of!* 

3436 FOR 1-1 TO 5 

3438 IF PG(JL)-1 THEN PRINT* *rLB»;JL!RBS! ELSE JL-JL+1 ! IF JL<*3 THEN 3438 

3442 IF PGIJD-l THEN JL-JL+1! NEXT I 

3444 PRINT 0965.*Indicate Option *!LBS;CHR*<176)!RBS; 

3446 OS-INKEYS ! IF OS**“ OR OS-CHRSU3) THEN 3446 ELSE PRINT Q 982,0»;!FOR 1*1 TO lOOHEXT 

3447 IF OS-*Q* OR OS**q“ THEN 40 ELSE FOR IH*1 TO 5 

3448 IF VAL(01)*IH AND FGOH1-1 THEN RETURN 

3449 NEXT IH ! GOTO 3444 
3430 L01-0 ! SM*0 

3460 FOR J*i TO IK 

3470 IF Ja OR X>*I+IH THEN SM«SM+RY(J) t L01*L0t+RT(JK2 
3400 NEXT J 

3490 L 01 *(LQ 1 - (SMC2)/(IN-IH))/aN-IH-l> 

3300 LOl*LOG(SQR(LOl»/LOG(10) 

3310 RETURN 

3530 FOR 1*1 TO IH I RYID-YCD ! NEXT ! IF IH-1 RETURN 
3333 FOR I” IN TO 1 STEP -1 
3360 IU*RND(IN) 

3563 RT-RYO)! RYCD-RYdU)! RYCIU)-RT 
3370 NEXT I 
3573 RETURN 
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NORMAL DISTRIBUTION 


5000 REM * • • NORMAL DISTRIBOTIOK *•# 

5010 SD-SQR(VAR) 

5020 xn=<nx-mu>/sd:ax»abs<xn> 

5030 XZ-.3989423»EXP((-AX*AX)/2) 

5040 NT* 1/(1 ♦•2316419*AX) 

5O3OQOX-XZ*«.319381S*NT)-<.3363638*NTt2)+(1.701478*NT[3>-(1.321236*NT[4>*<1.33O274*HTt3» 

5060 IF XNCO GOTO 5090 

5070 FOX-l-OOX 

5080 RETURN 

5090 FOX=QOX 

5100 QOX-l-FOX 

5110 RETURN 


\.=aswssim. 

6000 REM t-DISTRIBUTION *•« 

6010 IF TX<0 13*1 ELSE 13-0 
6020 IF DF>25 GOTO 6450 
6030 THETA-ATN(ABS(TX)/SQR(DF» 

6040 IF DFOt GOTO 6090 
6050 AT-.6366198*THETA 
6060 FOX-(AT+l)/2 
6070 QOX-l-FOX 
6080 GOTO 6400 

6090 IF INT(DF/2>»DF/2 GOTO 6220 

6100 IF DF-3 GOTO 6370 

6110 TS-0 ! TP-l 1 TC=C05(THETA) 

6120 FOR 1-1 TO INT((DF-3>/2> 

6130 TP-TF*2*I/(1+2«I) 

6140 TC-TC»(COS(THETA)C2) 

6150 TS-TS+(TP*TC> 

6160 NEXT I 

6170 TS-TS*C OS (THETA) 

6180 AT-.6364198*lTHETA*lSIN(TH£TA>*TS)> 

6190 FOX-(AT*I)/2 
6200 QOX-l-FOX 
6210 GOTO 6400 
6220 IF DF02 GOTO 6270 
6230 AT-SIM (THETA) 

6240 FOX-(AT*l)/2 
6230 QOX-l-FOX . 

6260 GOTO 6400 

6270 TS-0 S TP-1 1 TC-1 

6280 FOR 1-1 TO INT«DF-2)/2> 

6290 TP-TP*((2*IM)/(2*I) 

6300 TC-TC*(COS(THETA)t2> 

6310 TS-TS+(TP*TC) 

6320 NEXT I 

6330 AT-SIN(THETA>*(1+TS) 

6340 FOX»(AT*l)/2 
6350 QOX-l-FOX 
6360 GOTO 6400 

6370 AT-,6366198*(THETA+(SIN(THETA)*C OS (THETA))) 

6380 FOX*(AT*l)/2 

6390 QOX-l-FOX 

6400 PX-TX 

6410 IF 13-0 RETURN 

6420 QOX-FOX 

6430 FOX-t-QOX 

6440 RETURN 

6430 NX-( TX * (1 - >2S/DF) )/SQR(l ♦ .3*TXI2/DF> 
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mao mu=o : var* i: 12*1 : gosub sooo 

M70 PX-TX 
MBO RETURN 


IMVEE S S NO RM AL 

7000 REM *** INVERSE NORMAL DISTRIBUTION **♦ 

7010 IF PR>J THEN RL=SQR<LOG(l/«l-FR)t2)» ELSE HL*SOR(LOG( 1/IPRI2IIJ 

7020 XG=RL-<(2.51S517«-<.802853»RL>t-<.010328»RLC2»/<l'M1.432788»RLH-M89269*RLt2)+<.001308*RLt3)» 
7030 IF FR>*.5 THEN XG=MO+<SQR(VR)*XQ> ELSE XO*MO-(SQR(VR>*XQ> 

7040 RETURN 


fflVCT5E -t 

6000 REM *** INVERSE t - D I S T RI B U T I O N 
8010 MO-O ! VR-1 ! I2-V 
8020 GOSUB 7000 
8030 XOXQ 

8040 XG-1.23 * IXCI3 ♦ XCO/DF 

8030 X0*XQ ♦ (tl/96) * (5*XCI3 ♦ 16#XCt3 ♦ 3*XC»/DFI2 

8060 XQ*XQ ♦ (tl/384) * I3*XCC7 > 19*XCC5 ♦ 17«XCI3 -15*XC»/DFC3 

8070 XG-XQ + Ul/92160) # (79*XCC9 ♦ 776«XCC7 ♦ 1482*XCI5 -1920*XCt3 -94S«XC»/DFt4 

8080 XQ- XC + XQ 

8090 IF DF>3 RETURN 

8100 DELTA* 1 

8110 TX*XQ 

8120 GOSUB 6000 • 

8130 IF ABS(FR-FOX) < .00001 THEN XO-TX ! RETURN 

8140 IF PR < FOX THEN TX*TX-DELTA 

8130 IF PR < FOX THEN DELTA*DELTA/10 

8160 IF PR < FOX THEN GOTO 8120 

8170 TX-TX+DELTA 

8180 GOTO 8120 


PROGRAM SFl.ECnOK MENU 


10000 

10010 

10020 

10030 

10040 

10030 

10060 

10070 

10080 

10090 

10100 

10110 

10120 

10130 

10140 

10130 

10160 

10170 


REM **♦ PROGRAM SELECTION MENU 
CL5SFRINT TAB<20PPROGRAM SELECTION MENU- 

PRINT (3128,*61 Data I/O* Transformations, Frob. & Inv. Prob. Di»tntiution**lFRINT 
PRINT-62 Samp Is Stat's, 7-Num-Som, Plots*:FRINT 

PRINT"#3 Parametric Tests of Hypotheses and Confidence Intervals*:PRINT 

PRINT*64 Robust Estimates and Confidence Intervals-, PRINT 

PRINT* 63 Multiple Linear Regression" IPRINT 

PRINT-66 Analysis of Variance- 

KF-896 ! JCS*6 ! GOSUB 180 

IF 10*4 GOTO 20 

IF IQ-3 THEN 2T--2- 1 GOTO 10120 

IF 10*1 OR 10-2 THEN 2T«"1* ELSE ZT--3- 

CLS 

IF 10*1 OR 10-3 OR 10*3 THEN 10*1 ELSE 10-2 
CLS 

PRINT-INSERT TAPE 6*3Tr INTO ESF-80 AND PRESS SPACE BAH* 

IF INXETSO- * GOTO 10160 
OLOADt 


END OF PROGRAM 64 
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PROGRAM 


MENU 
5 CLEAR 

10 DIM MMU00.10), VCUOO.IO), VTVU0.100), RV(IOO), RHtlOO), EUOO), HI100), ID(20) 

20 DEFSTR 2 ! DEFINT I-K 

23 LBS-CHRSU23) S RBt*CHR*(123> ! BLS=CHRS(232> 

30 REM • ** PROGRAM DISCREPTION #•• 

40CLS 

30 PRINT’ • MULTIPLE LINEAR REGRESSION#* 

60 PRINT ! PRINT - This program fits a linear model to the data by the method of Ordinary Least Squares (OLS 
)• Parameters and statistics cal-culated include the equation coefficients, standardized coeffi- dents. *1 
70 PRINT ‘partial F statistics for the coefficients, correlation matrix. Hat Matrix (diagonal elements) and AHOV 
A table. 

80 PRINT 

90 PRINT* MODEL! T * BO ♦ B1 # XI ♦ B2 • X2 ♦ ... «• Bn # Xn 
100 PRINT ! PRINT* WHERE! Y * Response Variable 
110 PRINT- XI ... Xn - Carrier Variables 

120 PRINT* BO * Constant < Y Intercept ) 

130 PRINT* B1 ... Bn * Coefficients of Carrier Variables 

140 PRINT S PRINT*# PRESS SPACE BAR TO CONTINUE*! 

150 IF INKEYS O * * GOTO 150 
160 CLS 

170 FOR 1*10 TO 39 ! SET(I,3)! NEXT 
173 PRINT Q10,-PROGRAM LIMITATIONS*! 

180 PRINT 0>28.*(a) No more than 9 carrier variables 

190 PRINT @236,*(b) No more than 100 observations 

200 PRINT 0463,-MENU!*! FOR 1*30 TO 46! SET(I,24>:NEXT 

210 PRINT ©376,*<1) DatalVO 

220 PRINT ! PRINT*(2) Program Selection Menu* 

230 KF-896 i KS*2 ! GOSUB 300 
240 ON IO GOTO 400, 10000 


OPTION SELECTION 

300 REM #*• OPTION SELECTION 

310 PRINT 0KP*S,’OPTION DESIRED *!LBS!CHRS(176)iHBS! 

320 U-INKEYS ! IF ZI»** GOTO 320ELSE IO-VALOI) 

330 PRINT 0KP*21, II! ! FOR 1*1 TO 100 ! NEXT 

340 IF I0<1 OR IO>KS THEN PRINT 0KP, BLSJ ! PRINT 0KP,*### NOT A VALID OPTION *♦#*!! FOR >1 TO 1000 
1 NEXT ! PRINT 0KP, BL*!! GOTO 310 
330 RETURN 

360 PRINT 0KP#5,LB*!* *!LSJ* or *!RS!RBSi* *;LBs;CHRSU76)!RBs; 

370 ZQ-INICET* ! IF ZQ«** GOTO 370ELSE PRINT 0KP#17, ZQ!! FOR 1*1 TO 100 ! NEXT 
380 IF ZQOLS AND ZQOCHRS(ASC(LS)-32) AND ZQORS AND ZQOCHR»(ASCCRS>-32) THEN 360 
390 RETURN 


400 REM *•• IMFUT/OUTPUT OPTIONS ### 
410 CLS 

420 PRINT TABt20)*INPUT/OUTP U T OPTIONS* 

430 PRINT 0128,‘(1) Input Data via Keyboard 
440 PRINT*(2) Input Data via Tape 
430 PRINT-G) Save Data to Tape 
460 KP-896 ! KS-3 ! GOSUB 300 
470 OK IO GOTO 1000, 300, 700 



INPUT DATA VIA TAPE 


500 REM *»• INPUT DATA VIA TAPE »*♦ 

510 CLS l PHINT-Inaert DATA TAPE into ESF-80 and press SPACE BAR* 
520 IF INKEY*0“ * GOTO 520 

530 PRINT ©132,“Data File Number -;LS«:CHR*(176>;RB*J 
540 ZI-INKEY* 5 IF ZI*~ THEN 540! ELSE FI=VAL<ZI) 

550 PRINT 0150, ZI|! FOR 1*1 TO 100 ! NEXT 

560 IF no THEN 530 

570 ©CLOSE 

580 ©OPEN FI 

570 ©INPUT CV.IN.YM.Y2 

600 FOR 1*1 TO IN 

610 0INFUT RV(I) 

620 FOR J«1 TO CV 
623 IF IOl THEN 630 
628 0INPUT MU(J),X2(J> 

630 0INPUT MM(I,J) 

640 NEXT J 
650 NEXT I 
660 ©CLOSE 
665 TA*-*T* 

670 CLS ! PRINT-CALCULATING CORRELATION MATRIX...* J GOTO1230 


9 A YE PATA TQ TAPS 

700 REM *** SAVE DATA TO TAPE *•* 

710 CLS ! PRINT-Insert DATA TAPE into ESF-80 and press SPACE BAR- 
720 IF INKEY60- • GOTO 720 

730 PRINT 0132,-File Number to Save Data Under "JLBt!CHR*(176)JRB*; 
740 H-INXEY* ! IF ZI-“ THEN 740ELSE FI*VALOI) 

750 PRINT 0164, Zi;t FOR 1-1 TO 100 ! NEXT 

760 IF FIG THEN 730 

770 ©CLOSE 

780 ©OPEN FI 

770 ©PRINT CV,IN,TM,Y2 

800 FOR I-1 TO IN 

810 ©PRINT RV(I) 

820 FOR J-l TO CV 
825 IF IOl THEN 830 
828 ©PRINT MU(J),X2(J1 
830 ©PRINT MM(I<J) 

840 NEXT J 
850 NEXT I 
860 ©CLOSE 

870 PRINT ©512,'OPTIONS:* 

880 ©CLOSE 

890 PRINT* (1) Continue with Regression 
900 PRINT* (2) Return to Program Menu 
910 PRINT* (3) Go to Program Selection Menu 
920 KP-896 5 XS-3 ! GOSUB 300 
930 ON 10 GOTO 1500, 160. 10000 


INPUT DATA VIA KEYBOARD 

1000 REM INPUT DATA VIA EETBOARD *** 

1010 CLS , 

1020 INPUT "Number of Carrier Variables to be entered-JCV 

1030 IF CV>9 PRINT-*** ERROR **♦ NO MORE THAl^.9 CARRIERS ALLOWED"!GOTO 1020 
1040 IF CVG PRINT-EHROR *** MUST HAVE AT LEAST 1 CARRIER VARIABLE':GOTO 1020 
1043 FOR 1*1 TO CV 1 MU(I)*0 ! XS(I>*0 ! X2(IM> ! SX(I)*0 ! NEXT 1 
1050 INPUT*Enter the Number of Observations'ON 

1060 IF INOCV+1 THEN PRINT**** ERROR *** MIN. NUMBER OF OBSERVATIONS -*iCV+2 ! GOTO 1050 
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1070 CLS 

1083 YS-0 ! Y2-0 ! SY-0 : YM-0 

1090 FOR 1*1 TO IH 

1100 PRINT-OBSERVATION *“!l 

1120 FRINT"Enter the observed value for the response variable, Y’JSINFUT RV(I) 

1122 YM=YM+RV(I> 5 SY=RV(I)-Y5 : YS-YS+SY/I5 T2*Y2rSY«(RV<I>-YS> 

1130 FOR 7*1 TO CV 

1140 PRINT-Enter the Value for carrier variable, X*;MIDt<STES<7),2,2>;:lNFUT MM(ImT) 

1142 MUt7)*MU<7)+MM(I,7) i SX(J)*MM(I.7>-XSt7) 1 XS<J>-XSC7H'SX<7>/I ! X217)*X2(7H5X<7>*<MM(I»r>-XS(J» 
1150 NEXT 7 
1160 PRINT 
1170 NEXT I 


CALCULATE AND DISPLAY CORRELATION MATRIX 

1180 CLS : PRINT-CALCULATING CORRELATION MATRIX...* 

1190 YM-YM/IN 
1200 FOR 7-1 TO CV 
1210 MUt7)-MU<J)/IN 
1220 NEXT 7 
1223 TAS»*N* 

1230 FOR I>1 TO IN 

1240 IF TASO’Y* THEN RV(I)»(RV(I)-YM>/SQRCr2) 

1250 FOR 7*1 TO CV 

1260 VCa J>-(MMa J)-MU(7)I/SQR<X2(7» 

1280 NEXT 7 
1290 NEXT I 

1293 FOR 1-1 TO IN ! VCCI J)*RVa»! NEXT I 
1300 IV-CV+1 ! RETS*“Y" i GOSUB 2000 
1310 RET**-N- ! EPS*lE-3 

1320 CLS J PRINT-CORRELATION MATRIX! - ! PRINT 
1323 IF CV<-5 THEN PRINT* - ! ELSE PRINT* *! 

1330 FOR I-1 TO CV SPRINT* X*!RIGHT*(STRSdU>;! IF CV<-5 THEN PRINT * "J ELSE PRINT* *! 
1335 NEXT ! PRINT* Y* 

1340 FOR 1-1 TO IV 
1330 FOR 7-1 TO IV 

1355 IF CV73 THEN PS-**.###* ELSE PS-**.####* 

1360 IF 7*1 AND K*CV THEN PRINT *X*!RIGHT«(STRfa),l>!* *5 ELSE IF 7-1 PRINT *T *5 

1362 IF ABSfWaj)-lKEPS AMD CV>5 THEN PRINT * 1 *J! GOTO 1380 

1364 IF ABSCWajHlXEPS THEN PRINT* 1 *5 ! GOTO 1380 

1370 PRINT USING Pliwajl!! IF CV<-3 THEN PRINT* *1 ELSE PRINT* *1 

1380 NEXT 7 

1390 PRINT 

1400 NEXT I 

1410 PRINT 3969,''Press SPACE BAR to Continue*! 

1420 IF INXETSO* * THEN 1420 


SELECT SUBSET OF CARRIER VARIABLES 

1500REM *** SELECT SUBSET OF CARRIERS#*# 

1310 IF CV*1 THEN 1580 

1520 CLS • PRINT* The data array entered mntains*;CV;*carrier variables. You* 

1530 PRINT*may regress on aU*!CV!*carriers or any subset.* 

1540 PRINT 3197,LBSrA'iRBSJ'U or *;LB*rS"!RBS;*ubset *SLB*!CHR«176)!RB*! 

1550 HCt-INKEY* ! IF HC»=** THEM GOTO 1550ELSE PRINT 3221 ,HCS ! FOR 1-1 TO 100 ! NEXT 

1560 IF HCtO*A* AND HCSO'a* AND HCSO*S* AND HCSO's* THEN 1540 

1370 IF HC»**S* OR HC*-*s‘ THEN 1600 

1580 FOR 1-1 TO CV ! U*(I)**Y* ! NEXT ! GOTO 1810 

1600 REM CHOOSE SUBSET OF CARRIERS 

1620 CLS i IF DO-5 THEN 7D*S ELSE 7D-IN 

1630 FRINT*To facilitate recognition, the first*J7D!”values for each carrier* ! PRINT*are displayed below. When 
prompted indicate whether you want* 

1640 PRINT"to ■!LBs;*R*iRB»!*etain or *;LBt;"D*!RBS;"eletB «ach carrier for tire regression.* 

1650 KP-256 
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1660 FOR 1=1 TO CV 

1670 PRINT 0KP*3,“X“;MID*(STRS(I).2,2>5* “} 

1660 FOR J=1 TO JD 

1690 PRINT MMUMNIIF JOJD PRINT V5 

1700 NEXT J 

1710 KP=KP*64 

1720 NEXT I 

1730 KP-256 

1740 FOR 1*1 TO CV 

1750 PRINT 0KP,LB»;CHR»U76>;RB*; 

1760 RTS-INKEYl ! IF RT*=“* GOTO 1760ELSE PRINT @KP*l,RTt; 

1770 IF RT*0“R“ AND RTSO'r* AND RTSO"D“ AND RT*0“d“ THEN 1730 
1780 IF RT**“R‘ OR RTS»‘r* THEN US(I)»“Y' ELSE USdl-'N* 

1790 KP“KP*64 
1000 NEXT I 
1810 U«(0M*“Y* 

1820 IV-0 

1830 FOR 1=1 TO CV 

1840 IF U*<I)-*Y* THEN IV-IV-H 1 IWOVM 
1850 NEXT I 
I860 J-0 

1870 PRINT 0896,“Regress on! *1 
1880 FOR I-l TO CV 

1890 IF U»(I)«*Y“ THEN IP(J+1>-I: J-J+l 
1900 NEXT I 

1910 IF IV=0 THEN PRINT 0896,*•«« ERROR *« YOU HAVE DELETED ALL CARRIERS'JFOR I-l TO 2500 ! NEXT 
1 GOTO 1300 
1920 FOR I» 1 TO IN 
1930 FOR JT-1 TO IV 

1940 IF J-IV AND IC»IV THEN PRINT *X“:MIDS(STRS(IWttl),2,2>;*, *5 
1950 VC(I,J>»<MMU,IP(J»-MU(IP<J»)/SQR(X2<IF<J)» 

1960 NEXT J 
1970 NEXT I 

1980 PRINT 0960,'Calculating...'} 


CALCULATION OF OL5 ESTIMATES 


CAICVLATS V X 

2000 REM * * * CALCULATION OF LEAST SQUARES ESTIMATOR 
2010 REM BETA.hat - ( X' X l.inverse X' T 
2030 REM * * * CALCULATE X'X * * * 

2040 FOR >1 TO IV 
2050 FOR J«1 TO IV 

2060 vva,j>»o: vid^n-o 

2070 NEXT J 

2080 NEXT I 

2090 FOR I-l TO IV 

2100 FOR J-l TO IV 

2110 SUM-0 

2120 FOR IK-1 TO IN 

2130 SUM-SUM ♦ VCCIK,I)*VC(IXrT> 

2140 NEXT IK 
2150 VV(I,J>SUM 
2160 NEXT J 
2170 NEXT I 

2180 IF R£T*-*T* THEN RETURN 


CALCULATE (X'XHnverse 

2300 HEM * * * CALCULATE (X'X hnverse • * ♦ 
7310 FOR I-l TO IV 1 VKI.I1-1 i NEXT 
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3240 SUM»SUM+VTV<I,IK)*RV<IK) 
3250 NEXT IK 
3260 B(I) * SUM 
3270 NEXT I 


TRANSFORM BETA COEFFICIE NTS TO REGR ESSION EQUATION COEFFICIENTS 

3300 REM # * • TRANSFORM BETA COEFFICIENTS TO REGRESSION EQUATION COEFFICIENTS 
3305 B(0)»YM ! SC(0)»YM 
3310 FOR 1=1 TO IV 
3320 sca>=sa> 

3330 Ba>=SC(I)»SQR(Y2>/SQR(X2HF<I») 

3340 B(0>»B(0)-BtI)*MU<IPa» 

3350 NEXT I 


CALCULATE ELEMENTS OF ANQVA TABLE 

3500 REM **• CALCULATE ELEMENTS OF ANOVA TABLE 

3510 RSS=»0 : ESSO i E2S*0 

3520 TS5=Y2 

3570 FOR 1=1 TO IN 

3580 SUM-O : S2MO 

3590 FOR IK«1 TO IV 

3600 SUM«SUM+MM(I,IP<IX»*BaX> 

3605 S2M«S2M+VCa,IK)*SCIIX) 

3610 NEXT IK 
3620 SUM*SUM+B<0) 

3630 RH(I)*SUM 

3640ESS«ESS+(RVa)«SQR(Y2HYM-SUM)*(RVa)*SQR{T2H-YM-SUM) 

3645 E2S«E2S+(RVa>-S2M>«<Rva>-S2M) 

3650 ECI>=<RV<I>«SQR<Y2>+YM)-SUM 
3660 NEXT I 
3665 RSS-TSS-ESS 
3670 MRS-RSS/aV) 

3680 M£S=ESS/GN-IV-1) 

3685 M2S»E25/(IN-IV-1) 

3690 F-MRS/MES 

3700 R2-RSS/TSS 

3710 FOR 1-1 TO IV 

3720 T(I>-SCa>/(SQR(M2S * Via,I»> 

3770 SEHJ-SQRIM2S * VM4» • SQH(Y2)/SQR(X2aFH)» 

3780 NEXT I 


DISPLAY ANOVA TABLE. REGRESSION EQUATION. OTHER STATISTICS 

4000 REM PRINT OLS SOLUTION, ANOVA, & OTHER STATISTICS 

4010 CLS 

4020 PRINT 026,* A NOVA* 

4030 PRINT* Source DF Sum of Squires Mein Squire F-Ritio 

4040 FOR 1=5 TO 15! SETH,6)1 NEXT! FOR 1=26 TO 29S SETH,6>: NEXT! FOR I«40 TO 671 SETa,6>: NEXT! FOR 1-78 
TO 99! SETH,6)1 NEXT! FOR 1=106 TO 1201 SETtI,6)i NEXT 

4050 PRINT 0192,*RegnMBion*[lPRINT 02O4,IV;:PRINT 0213,RSS[! PRINT @231,MRS5IPHINT 0243,* *[HF F7999.99 

9 THEN PRINT USING*##.#**##CttC*:F ELSE PRINT USING****#.#*#*[F 

4060 PRINT 8258»*Error" [SPRINT 0268,IN-UVM SPRINT 0277,ESS;:PRINTQ295,MES 

4070 PRINT 0322»“Total*,!FRINT 0332,IN-1 [SPRINT 0341.TSS 

4080 PRINT : PRINT "The Regression Eqmtion i»l* 

4090 PRINT* Y * *J 
4100 FOR I»0 TO IV 
4110 IF I O 0 GOTO 4140 

4120 IF B(0)>0 THEN PRINT * *S ELSE PRINT *- *} 

4130 GOTO 4150 

4140 IF BH) 7= 0 THEN PRINT* ♦ *} ELSE PRINT* - *) 

4150 PV=A3S(Ba» 
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4160 IF PV>9999.9999 OR FVC .0001 THEN PUs»*e.###*CCtC* ELSE PU*=*##*#«.eeee* 

4170 PRINT USING PU*;PVI 

4180 IF 100 PRINT* X*;MID*<STR*dWd»,2,2>;* ", ELSE PRINT * ", 

4190 IF 1=2 OR 1=5 OR 1=8 THEN PRINT" * 5 PRINT" ", 

4200 NEXT I 

4205 IF IVC=5 THEN PRINT 3704,"“; ELSE PRINT (?768,"*; 

4210 PRINT *R-Square! "SPRINT USING" ».**#**;R2;;PRINT* * I* Adjusted R-Square! *j:PRINT USING*# 

.•••••;i-U-H2>*aN-l)/(IN-IV-l) 

4220 PRINT’Std. Error:*::PRINT U5ING"#ee#e.ee#e*;SQR(MES> 

4225 PRINT"Std. Error as % of Response MeanJ’I'.PRINT USING**#ee*e#.##e";iSQR(MES>/YM>*100;:FRINT*\*! 
4230 PRINT G?60,"PRESS SPACE BAR TO DISPLAY ADDITIONAL STATISTICS*; 

4240 IF INKEYSO* * GOTO 41'0 

4250 REM ♦ ♦ * PRINT COEFF's, t-STAT's. STD. ERROR of COEFF's, & STANDARDIZED COEFF's * * * 

4260 CLS 

4270 PRINT*Var. Coeffiaent Partial F Std. Error Beta Coeff" 

4280 PRINT*- - - - -* 

4290 Ul*=*##*##*#***.#*##* J U2»=“#.*#*#lCtt* 

4300 FOR 1=1 TO IV 

4310 PRINT“X*iMID*<STR*dW(IH,2,2); 

4320 IF A B5tB<I)»999999 OR ABS(B(I»<.0001 THEN PRINT * ‘JSPRINT USING U2*iBd)i:GOTO 4340 

4330 PRINT USING U1*;B(I>; 

4340 IF ABS(TG)«TQ)»999999 OR ABS<T(I)*T(I»<.0001 THEN PRINT" ‘"PRINT USING U2*rni)«Ta>}: GOTO 43 
60 

4350 PRINT USING Ul*rr<I)«TtI>! * 

4360 IF ABSiSEd >>>999999 OR ABSISEtDK.OOOl THEN PRINT" *!iPRINT USING U2SS SEdKiGOTO 4380 
4370 PRINT USING Ul*! SEd)! 

4380 IF ABS(SCd)»999999 OR ABS(SC(I)X.0001 THEN PRINT" *;!PRINT USING U2*J SC (I) : GOTO 4400 
4390 PRINT USING Ul*; SC(I> 

4400 NEXT I 

4410 PRINT G896,“Da you want to see the Diagonal Elements of the Hat Matrix?*; 

4420 L*=*y* ! R*»‘n* : KP*960 ! GOSUB 360 

4430 IF ZQ="N* OR ZQ=*n" THEN PRINT 8896, BL*i! PRINT 3960, BLSJIGOTO 4970 


PISPLAT BIAfiPNAl QF MATO 

4500 REM DISPLAY DIAGONAL ELEMENTS OF HAT MATRIX **« 

4510 CLS 

4520 PRINT * The diagonal elements of the Hat Matrix ( H = XtX'XJinv. X' )give a measure of the influence of th 
e i th observation of theresponse variable on the corresponding fitted value.*! 

4530 PRINT* A 'large'value of h(i,i> is indicative of an outlier. This, of course, does not imply that the data po 
int is necessarily a 'bad' pointAut should alert you to inquire into its nature." 

4540 PRINT : PRINT* A reasonable rule of thumb for what is considered a large H(14> is any value larger than 
twice the rank of the hat ma- trix divided by the number of observations.* 

4550 PRINT 3704, *. i.e. Consider HU,i> >*J2*IV/IN;* as large* 

4560 IF LH»=*N* GOTO 4590 

4570 PRINT 0834,’NOTES A blinking cursor will highlight large' H(i,i)'s* 

4580 FOR JG-1 TO 10 ! PRINT G832,CHR*(143>;: FOR JH=t TO 25 ! NEXT JHJ PRINT 3832.* *!5 FOR JH-1 TO 25 ! 
NEXT JH !NEXT JG ! PRINT 3832,**; 

4590 PRINT 0?60,*Press SPACE BAR to display the H<i,i> 's"{ 

4600 IF INKEY* O * * GOTO 4600 
4610 CLS 

4620 FRINT’Diagonal Elements of Hat Matrix:* 

4630 FOR 1=1 TO IN STEP 3 

4640 IF INTO/37) O 1/37 AND INTd/73) O 1/73 GOTO 4700 
4650 IF LH*“*Y* GOSUB 4810 

4660 PRINT 3?60,*Press SPACE BAR to display remaining elements*, 

4670 IF INKEY* O * * GOTO 4670 
4660 CLS 

4690 PRlNT*Diagonal Elements of Hat Matrix (continued)!* 

4700 FOR J=0 TO 2 
4710 IF I*J > IN GOTO 4790 
4720 IF J«0 GOTO 4740 

4730 IF J»1 PRINT TAB(21>! ELSE PRINT TAB(42)J 
4740 J*=MID*(STR*d+J>,2,3> 

4730 frint*hc;j*;v;j*{*> «*;Hd+J>! 
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4760 NEXT J 
4770 PRINT 
4780 NEXT I 

4790 IF LH*»’Y“ GOSUB 4810 
4800 GOTO 4980 
4810 FOR JX-1 TO JL 

4820 IF I<-37 AND IDIJXK37 THEN GOSUB 48705 GOTO 4850 

4830 IF I>37 AND IO 73 AND ID(JX»=*37 AND IDUXK73 THEN GOSUB 4870! GOTO 4850 

4840 IF 1373 AND ID(JX)3-73 THEN GOSUB 4870 

4850 NEXT JX 

4860 RETURN 

4870 FOR JZ=l TO 10 

4880 IF K-37 THEN IF-IDCJX) 

4890 IF D37 AND I<*73 THEN IP*ID(JX>-36 
4900 IF 1373 THEN IP-ID<JX>-72 

4910 IF IP/3-INT1IP/3) THEN PP-<IP/3>*64 ELSE PP-aNTOP/3Hl)*64 
4920 IF IP/3-INTGF/3) THEN PP-PP+42 5 GOTO 4940 
4930 IF (IP*l>/3 » INTUIPM>/3> THEN PP-PP+21 

4940 PRINT GPP.CHR4U43)55 FOR JY-1 TO 25 5 NEXT JT 5 PRINT GPP,*H*J 
4950 NEXT JZ 
4960 RETURN 

4970 PRINT G896,*Do you want to SB* the residuals?'? 

4980 PRINT G896,*Do you want to ses the residuals?'; 

4990 KF-960 5 L*-*y* ! R*»*n* 5 GOSUB 360 
5000 IF ZQ-*Y* OR IQ»*y* THEN 6010 

5010 PRINT @896■ BLS 5 PRINT G896,"Do you want to see a plot of the residuals?*; 

5020 GOSUB 360 

5030 IF ZQ-*Y* OR ZQ-'y' THEN 6510 
5040 GOTO 6910 


DISPLAY RESIDUALS 

6000 REM DISPLAY RESIDUALS **♦ 


6010 CLS 

6020 PRINT* T Observed Y Estimate Residual* 

6030 FOR 1-1 TO IN 

6040 PRINT USING***##**#.*##*;RV(I)*SQR(Y2)+'YM;5PRINT* 'J5PRINT USING* 


(IHJPRINT* *;:PRINT USING*4#**#**##.e#**;E(D 
6050 IF 1/12 O INTO/12) GOTO 6100 
6060 IF I-IN GOTO 6100 

6070 PRINT5PRINT“PR£SS SPACE BAR TO DISPLAY REMAINDER OF RESIDUALS* 
6080 IF INKEYS O * * GOTO 6080 

6090 CLS5PRINT GO,* Y Observed Y Estimate Residual* 

6100 NEXT I 

6110 PRINT 0896,"Da you want to see a plot of the residuals?*; 

6120 KP-960 5 L»-‘y* 5 RS»‘n* ! GOSUB 360 
6130 IF ZQ-*Y* OR ZQ-*y* THEN 6510 
6140 GOTO 6910 


DISPLAY PLOT OF. RSSTPVAL? 

6500 REM *** PLOT OF RESIDUALS *** 

6510 CLS SPRINT 90,'THE HORIZONTAL AXIS REPRESENTS Y ESTIMATE 

6520 PRINT 9192,'THE VERTICAL AXIS REPRESENTS RESIDUAL OR ERROR OF ESTIMATE 

6530 TS«1£*25 5 YB-1E-25 5 RL-YS 5 RB-YB 

6540 FOR l-l TO IN 

6550 IF HHtIKYS THEN YS-RH(I) 

6560 IF RH<I)3YB THEN^B-RHtt) 

6570 IF EdKRL THEN RL-E(I) 

6580 IF E(I)>RB THEN RB-EG) 

6590 NEXT I 

6600 PRINT G64,* RANGE OF Y ESTIMATE -*;YS;*TO*5YB 
6610 PRINT Q256,“ RANGE OF RESIDUALS »*;RL;*TO*;RB 

6620 PRINT G768,“«** AFTER VIEWING THE PLOT, PRESS SPACE DAR TO CONTINUE* 


>*JRH 
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6430 HU=tia/(YEJ-YS> ! VU=48/(R8-KL) 

6640 IF YS>*0 THEN VA-10 ELSE VA=10+ABS<YS«HU> 

4630 IF RL>=0 TIEN HA=47 ELSE HA=47-ABS(RL*VU> 

6440 PRINT 0896,"PRESS SPACE BAR TO DISPLAY PLOT- 
6470 IF INKEYS > • * GOTO 6670 
6480 CLS 

6490 FOR 1=0 TO 47 
6700 SET(VAtl) 

6710 NEXT I 

6720 FOR I« 10 TO 127 

6730 SET a,HA) 

6740 NEXT I 

6730 Kl=INT<VA/2>-5 ! K2=960+INT(VA/2>-3 : K3-INT1HA/3+1.3) i IF HA/3 - INT(HA/3)>-2/3 THEN K3-K3-1 
4740 IF K3M3 THEN K3-13 

4770 PRINT 0Kl,-j: IF ABS<INTIRB+.5)»=t THEN PRINT STR*<INT(RB+.3)>; ELSE PRINT USING*9.994*IRB; 
6780 PRINT SK2,-;: IF ABS(INT(RL+.3)»=t THEN PRINT 5TR$(INT(RL+.3»; ELSE PRINT USING" •••♦♦"IRL; 

6790 IF VA<»14 AND YS<0 THEN K4=K3*64 ELSE SC4*K3»64+6 

6800 PRINT 0K4,-;:iF ABS<INT<YS+.5»)=1 THEN PRINT STR*<INT<YS+.5»; ELSE PRINT USING*9.8994*;TSS 
6810 IF ABS<INT<YB+.5»>=1 THEN PRINT 0K3»64+63-LEN(STRS(INT(YB+.3)»,STRt(INT(YB+.3)); ELSE PRINT 010 

la.-jiFRiNT using*9.949*;yb; 

6820 FOR I»1 TO IN 

6830 O - 10+(RHU)-YS)*HU * D » 47 - t (E<D-RL>*VU ) 

6840 IF 0)127 THEN 0=127 

6830 IF D<0 THEN D»0 

6840 IF 0=INT<VA) THEN 0=0+1 

4870 IF D-INT(HA) THEN RESETtO-l,D)l IF 00127 THEN RESET(0+1,D) 

6880 SET <0,D) 

4890 NEXT I 

6900 IF INKEYS O * * GOTO 6900 

6910 PRINT 0832, BLS I PRINT BLS IPRINT BLS; 

6920 PRINT 0832,“(1) RETURN TO ANOVA TABLE*; 

6930 IF CV)1 THEN PRINT 0896,*<2> RERUN REGRESSION WITH A DIFFERENT SUBSET OF CARRIERS*} 

6940 IF CV)1 THEN PRINT 0960,'’OH ELSE PRINT 0960,*(2)*} 

6930 PRINT 0964,-RETURN TO PROGRAM MENU*} 

6940 KP-768 ! IF CV)1 KS-3 ELSE KS=2 
6970 GOSUB 300 

6980 IF 10=1 THEN 4010ELSE IF 10=2 AND CV«1 THEN 160ELSE IF 10-2 THEN 1500ELSE 160 


PBQflRAM SRSCTIOM MSNV 

10000 HEM *•* PROGRAM SELECTION MENU 
10010 CLS 

10020 PRINT TAB(20)*PROGRAM SELECTION MENU* 

10030 PRINT 0128,*91 Data I/O, Transformatons, Prob, & Invars* Prob, Distributions* 

10040 PRINT**2 Sample Stat's, 7-Num-Sum, Plots’iPRINT 

10050 PRINT**3 Parametric Tests of Hypotheses and Confidence Intervals*SPRINT 
10040 PRINT*#4 Robust Estimates and Confidence Intervals-tPRINT 
10070 PRINT*#3 Multiple Linear Regression" IPRINT 
10080 PRINT*94 Analysis of Variance" 

10090 KP*896 : KS«4 5 GOSUB 300 

10100 IF IOC1 OR 10)6 PRINT 0896,CHRS(25O)!:FHINT 0896,**** NOT A VALID OPTION WjaFOR 1=1 TO 1000J 

NEXT ! PRINT 0896X HRS (230IJIGOTO10090 

10110 IF 10=5 GOTO 30 

10120 IF IO-4 THEN 2T=*3* ELSE 2T="1* 

10130 CLS 

10140 IF 1006 THEN PRINT "From program 93 acxess to programs 91 thru 94 is through the Control Module in 
program 91,* ! PRINT 

10130 PRINTTNSERT TAPE 9*!ZT}* INTO ESF-80 AND PRESS SPACE BAR* 

10160 IF INKETSO* * GOTO 10160 
10170 IF 10=6 THEN 10=2 ELSE 10=1 

10180 IF 10=1 THEN PRINT ! PRINT ’Program 91 is now being loaded...* 

10190 0LOAD 10 


END OF PROGRAM 95 
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PROGRAM 


MEMJ 


3 CLEAR 
10 CLEAR 200 
20CLS 

30 DEFINT I-K ( DEFSTR L 

33 LB«*CHR«(123) : MB**CHR*<174) i RB»«CHRS(123) *. BLMCHR«(232) 
40 DIM PSUQO), C(10,20) 


30 PRINT TAB(20)(* Analysis of Variance* 

40 PRINT TAB(26>!’(ANOVA)“ 

70 PHINT(PRINT"Menu(“ 

80 PRINT* (1) One-Way Classification 

90 PRINT* (equal or unequal no. of observations for each factor)* 
100 PRINT* (21 Two-Way Classification 
110 PRINT* (with or without replication) 

120 PRINT* (3) Three-Way Classification 

130 PRINT* (4) Return to Program Selection Menu* 

130 LW«CHR*(233)( LP»* A###*#.###" ( LD**#** 

140 KP-894 ( KS*4 ! GOSUB 1000 ( WS-IO ! XZ*0 
170 CLS 5 IF WS*4 THEN 20000 

180 PRINT ‘Will Data be Entered via *!LBS!*K*;RB*;*eyboard or *; 
**t* 1 GOSUB 1200 

190 IF Q***K* OR Q*-*k* THEN n*l ELSE 11*2 I GOSUB 4000 
200 CLS 

210 ON WS GOSUB 3000,10000.13000,20000 


LB*(*T*5RB*;*ape ?*!: KP-128 1 L*«*k* 1 R* 


PRINT ANQVA TABLE 
213 CLS 

220 IF WS-1 LH»*(One-Way Classification)*(GOTO 240 
230 IF WS-2 LH*‘(Two-Way Classification)*:GOTO 240 
240 IF WS-3 LH**(Three-Way Classification)* 

240 PRINT TAB(15)l*ANOVA “JLH 

270 PRINT* SOURCE D.F. SUM OF SQUARES MEAN SQUARE F 
290 PRINT STRING*(42,*-*) ( ON WS GOSUB 300,400.800 

283 CLS ! PRINT'Do you want to save the data to tape 1 *-;: KP*37 ( L*-*y* : R*-'n* 1 GOSUB 1200 : LQ-Q* 

290 IF LQO*y* AND LQO’Y* AND LQO'n* AND LQO'N" THEN GOTO 283 
293 IF LQ**y* OR LQ»*Y* GOSUB 7000 
297 CLS (GOTO 50 • 

300 PRINT LCrrABtl3)5CDFrTAB(19)5:PRIMT USING LPTTRTSS(PRINT TABQWCPRINT USING LPJMCJIFHINT TAB(3 
01NPRINT USING LPSF 
310 PRINT 

320 PRINT*R»sidual"!TAB(13)ffiDFrrAB(19)i:PRINT USING**##***.****;ESS;:FRINT TAB<34)!(FRIHT USING'*# 

m**.»#*-;me 

330 PRINT 

340 PRINT STRING*<42,*-*) 

330 FRINT*Total*fTAB(13)!TDFJTAB(19)!:PRINT USING* #*#***.#***?rTS 
340 PRINT 0940,"Press SPACE BAR to continue*! 

370 IF INKEY* O * * GOTO 370 
380 CLS : RETURN 

400 IF I EM PRINT *Interaction*;TAB<13);:PRINT USING LDJADFJIPRINT TAB(19);iFRINT USING LFJTISSNPRINT 
TAB(34)i (PRINT USING LP(MA((PRINT TAB(30)i(PRINT USING LPI FA (PRINT 
410 IF IE*l PRINT 9192,■*! ELSE PRINT Q320,**! 
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420 PRINT LR;TA3<t3K:FRINT USING LDiRDFUPRINT TABU?»;:PRINT USING LP!TRSS!!PRINT TAB(36);;PRIHT u 
sing lp;mr;:print taS(50>;:priht using lpifr:print 
430 IF IE*1 PRINT 0320,""! ELSE PRINT 0443,"“} 

440 PRINT LCrrAB(13i;:PRINT USING LD!CDF!!PRINT TABU9);!FRINT USING LPJTCSS!SPRINT TAB(36>;:PRINT u 

sing lp;mc;:print tab<50>;:print using lp;fc:frint 

450 IF IE*1 PRINT 0448,""! ELSE PRINT 0576,“'! 

460 PRINT 'Residuar;TAB<13)!!PRINT USING LDiEDFJSPRINT TAB(19)!!PRINT USING LPJESSJIPRINT TAB<36>i!FR 

INT USING LPJME 

470 PRINT STRING»<62,*-*» 

480 IF XZ*1 GOTO 670 

490 PRINT"Total’!TAB<13>!:PRINT USING LDJTDFJiPRINT TABU9KSPRINT USING LP.TTS 
500 IF IE*t GOTO 670 

510 FOR 1*1 TO 25 ! PRINT 0253,CHR«<139> J FOR J*1 TO 25 i NEXT J 1 PRINT 0253,* * ! FOR J*1 TO 25 i NEXT 
J 1 NEXT I 

520 PRINT 0832,“REMINDER! a Urge F-statistlc «> rejection* 

530 FOR 1*1 TO 1500 1 NEXT I 
540 PRINT 0832, LW ! PRINT LW 

550 PRINT 0832,'Do you accept or reject the hypothesis that there is no inter-*!PRINT*actian between the row a 
nd column factors?"’,!KP*936 ! L»**a" 1 RS-'r* ! GOSUB 1200 ! LOGS 
570 PRINT 0832,LW 1 PRINT LW 
580 IF LO*a* OR LQ-'A* GOTO 620 

590 PRINT 0768,'Since you reject the hypothesis of no interaction the linear*!PRINT“model is violated and the re 
mainder of the ANOVA is invalid"!PRINT**** Press space bar to continue*; 

600 IF INKEY* O * * GOTO 600 

610 CLSiRETURN 

620 PRINT 0832, LW ! PRINT LW 

630 PRINT 0192, LW ! PRINT 0192,'Hypothesis that there is no interaction is accepted* 

640 PRINT 0832,“Do you want t'o pool interaction fr residual?*! LS**y* ! R»**n" ! KP*873 ! GOSUB 1200 1 LQ*GS 
640 IF LQ - V OR LQ **Y* GOTO 700 

670 PRINT 0768, LW ! PRINT LW ! PRINT LW 1 PRINT 0768,LRi* t "!LC1* hypotheses can be assessed using" 1 PR 
INT*the F-ratios shown above,"! PRINT* < Press space bar to continue)*! 

480 IF INKEY* O * " GOTO 680 
690 CLS t RETURN 

700 ESS*TISS ♦ ESS ! EDF«ADF*EDF ! ME-ESS/EDF ! FC*MC/ME ! FR-MR/ME ! XI* 1 ! GOTO 410 

800 PRINT LHJTAB<13){!PHINT USING LDJRDFJIFRINT TAB<19)!!PRINT USING LPrTRSSSlFRINT TAB<36)!!FRINT U 

SING LPIMR)!PRINT TABt50)!!PRINT USING LP! Ft 

810 PRINT LC!TAB<13)!!PRINT USING LD!CDF!!PRINT TAB(19)!!PRINT USING LPJTCSSNPRINT TABQ6WPRINT U 

sing lp;mc;:print tabisokjprint using lpi F2 

820 PRINT LG!TAB(13);:PRIKT USING LDIGDFJJPRINT TAB(19)|!PRINT USING LPJTGSSilPRINT TAB(36);iFRINT U 

sing lp;mg;sprint tab<so);:print using lp; fs- 

830 PRINT ! PHINT*Interactions!' 

840 PRINT LEFT*(LR,5)!*-*!LEFT*(LC,5>!TAB<13)ja s RINT USING LDJDDFJJPRINT TAB(19)J!PHINT USING LP!DOS8 
51PRINT TAB(36);:PRINT USING LP!HD!!PRINT TAB(50)i!PRINT USING LP! F4 

850 PRINT LEFT*(LR,5)f-*!LEFT*(LG^)!TAB(l3);:PRINT USING LD!BDF!!PRINT TABU9)J!PRINT USING LPJDBSS 
ItFRINT TAB(36)i!PRINT USING LPJMBiJPRINT TABt50)!!PRINT USING LPJF5 

860 PRINT LEFT*(LC^)!*-*!LEFT*(LG^)!TAB(t3);:PHINT USING LD!ADF!!PRINT TAB<19I!1PRINT USING LPiAISS 
{SPRUIT TAB(36)!!PRINT USING LP!MA!!PRINT TAJB(50)!!PRINT USING LP! F6 

870 PRINT ! PRINT"Residual*!TABU3)i!PRINT USING LDiEDF!!PRINT TAB(19>!:PRINT USING LPiESSJIPRINT TAB 
(3A)!.'FRINT USING LP! ME 
880 PRINT STRING»<63,*-*> 

890 PRINT "Total*!TAB(13)!!PRINT USING LD«TDF!SPRINT TAB(19)i!PRINT USING LP!TSS 
900 PRINT 0960,“»«• Press SPACE BAR to return to continue*! 

910 IF INKET*0* - GOTO 910 
920 CLS 
930 RETURN 


OPTION 5KSCTON 

1000 PRINT 0KP+5,'OPTION DESIRED *!LB*!MB*!RB*i 
1010 I*»INKEY* ! IF I***“ THEN 1010 ELSE IO«VALH*> 

1020 PRINT QKP+21, I*t ! FOR 1*1 TO 100 ! NEXT 

1030 IF 10<1 OR I07KS THEN PRINT 0KP, BL*{! PRINT 0KP,**** NOT A VALID OPTION *•*■}! FOR 1*1 TO 100 
0 ! NEXT ! PRINT 0KP, BL*!i GOTO 1000 
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1040 RETURN 

1200 PRINT 0KP+3(LB*5“ *iL*J* or "IRS!* "IRB*!* *:LB*iMB*;RB*! 

1210 QS-INKEYS ! IF Q*="“ THEN 1210 ELSE PRINT 0KMO,Q*;: FOR I-t TO 100 : NEXT 
1220 IF Q«OL* AND Q*OCHR*(ASC(L*)-32) AND QIORI AND QSOCHR»(ASC(R*>-32> THEN 1200 
1230 RETURN 


ONE-WAY CLASSIFICATION 

3000 REM **♦ ONE-WAY CLASSIFICATION 
3003 IF 11*2 THEN ON QL GOTO 5110, 5390 

5010 PRINT TAB(27)!“ANQVA" ! PRINT TAB(17)S‘(One-Way Classification)' 

3020 PRINT 0128,"Enter a reference name for the factor (i>a. Operator"(INPUT"(enter up to 12 character*)"JLC 
3030 PRINT 0128.LW 5 PRINT LW 

3040 PRINT 0128,“Enter the number of "!LC’,* levels (up to lOI'JSINPUT IC 

3050 PRINT 0128,LW : PRINT 0128,'Are there an equal number of experiments / observations' t PRINT"for each 
level?'! 

3032 KP-203 S L*-'y“ t R*="n' : GOSUB 1200 : LQ-Q* 

3033 IF LQ»"y" OR LO-'Y" THEN QL-1 ELSE QL-2 
3060 PRINT 0128.LW : PRINT LW 

3080 IF LQ-*n" OR LQ-"N* GOTO 3390 
3090 PRINT 0128,LW : PRINT LW 

3100 PRINT 0128,"Enter number of observations for each 'JLC,' (up to 20)'GNFUT IE 
3110 S«0 SOS-O 

3113 CLS ! IF 11-2 THEN PRINT 00,"Calculating" 

5120 FOR >1 TO IC 

3130 IF II-l THEN PRINT LC!I 

5140 CS(I>-0 

3130 FOR J«1 TO IE 

3160 IF H-l THEN PRINT' Enter value of observation'JJJONPUT Cd,J> 

3170 csa)-cstiHca.j> 

3180 S-S+CCI.JK2 

3190 os-os+caj) 

3200 NEXT J 

3210 IF n -1 THEN PRINT 

3220 NEXT I 

3230 CLS : PRINT 00,"Calculating' 

3240 OS«OPI2/<IC»IE) 

3230 TRTS-0 
3260 FOR 1-1 TO IC 
3270 TRTS-TRTS-*CS<IK2 
3280 NEXT I 

3290 TRTS«(TRTS/IE)-OS 
3300 TTS-S-OS 
3310 ESS-TTS-THTS 
3320 CDF-lC-l 
3330 EDF»IC»(IE-1) 

3340 TDF*IC«EE-1 
5330 MC-TRTS/CDF 
3360 ME-ESS/EDF 
3370 F-MC/ME 
3380 CLS 5 RETURN 

3390 CLS : IF H»2 THEN PRINT 00,'Calculating' 

3400 os-o ; s«o : cob-o : jn-o 

3410 FOR I-l TO IC 

5420 IF n«l THEN PRINT'Enter the number of observations for 'JLCJIS'tup to 20)'HNPUT IXG) 

3430 csai-o : catiw): jn«jn+ix<d 

3440 FOR J-l TOIXtl)* 

3430 IF H-l THEN FRINT'For "ILCJi;" enter observation'JJJJINFUT C(I,J> 

3460 cs<n«csn>cti,j> 

3470 os-os+caj) 

3480 S«S*C(I,J)t2 

3490 NEXT J 
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5500 PRINT 
5510 NEXT I 

5520 CLS ! PRINT GO,"Calculating" 
5525 TRTS*0 
5530 FOR 1*1 TO IC 
5540 TRTS=TRTS*CS<I)C2/IX(I> 
5550 NEXT J 

5560 TRTS-TRTS-05C2/JN 
5570 TTS*S-Q5C2/JN 
5580 ESS*TTS-TRTS 
5590 CDF*IC-1 
5600 EDF-JN-IC 
5610 TDF*CDF+EDF 
5620 MC*TRTS/CDF 
5630 ME-ESS/EDF 
5640 F»MC/ME 
5650 CLS 1 RETURN 


LOAD DATA FROM TAPE 

6000 HEM »*• LOAD DATA FROM TAPE **# 

6010 CLS 

6020 PRINT-NOTE! You must specify a file number containing data which was" 

6030 PRINT*recorded in the format for *! 

6040 IF WS*1 PRINT-one-way classification.- ! GOTO 6100 

6050 IF WS*2 THEM PRINT*two-way classification." ELSE PRINT "three-way classification* 

6100 PRINT ! PRINT"Insert Data Tape into ESF-80 and press SPACE BAR" 

6110 IF IN KEY SO" • GOTO 6110 

6120 FRINTJINPUT "Enter FUe • of Data to be Loaded'tfF 

6130 GCLOSE 

6140 QOPEN JT 

6150 SS-WS 

6160 GIN PUT WS 

6170 IF WS-SS GOTO 6300 

6180 CLS I PRINT-ERROR... the data file you are attempting to load is for * 

6190 LLlll-'one" 1 LL(2)»"two" ! LLGP>*three* 

6200 PRINT LL<WS>S*-way classification instead of "»LL<SS)i'-way classification!" 

6210 PRINT 1 PRINT "Da you want to continue to load data and run a "!LL(WS)!"-way" 

6220 PRINT"dassification ANOVA or da you want to start over?" 

6230 PRINT Q896,-;aNPUT" (1) Continue (2) Start Over 11 or 2 )"!IO 

6240 IF 1001 AND 1002 PRINT 0896, LW ! PRINT Q896,"»** NOT A VALID OPTION **«"! FOR 1*1 TO 1000! NE 

XT 1 GOTO 6230 

6250 IF 10-2 HUN 

6300 ON HS GOTO 6310, 6500, 6700 

6310 GIN PUT LC 

6320 GINPUT QL 

6330 GINPUT IC 

6340 IF QL*2 GOTO 6400 

6350 GINPUT IE 

6360 FOR I-1 TO IC 

6370 FOR J-l TO IE 

6380 ginput ca,j> 

6385 NEXT J 
6390 NEXT I 
6395 GCLOSE 
6397 RETURN 
6400 FOR I-1 TO IC 
6410 GINPUT IX (I) 

6420 FOR J»1 TO IXU1 
6430 GINPUT Ca ,J) 

6440 NEXT J 
6450 NEXT I 
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6460 GCL05E 
6470 RETURN 
6500 ©INPUT LH.LC 
6510 ©INPUT IRiIC.IE 
6520 FOR 1*1 TO IR 
6530 FOR J*1 TO IC 
6540 FOR K*1 TO IE 
6550 ©INPUT XdtJ.K) 
6560 NEXT E 
6570 NEXT J 
6580 NEXT I 
6590 ©CLOSE 
6600 RETURN 
6700 ©INPUT LR.LC.LG 
6710 ©INPUT IXtIR.IC'IG 
6720 FOR 1*1 TO IR 
6730 FOR J*1 TO IC 
6740 FOR K*1 TO IG 
6750 ©INPUT XU.J.K) 
6760 NEXT K 
6770 NEXT J 
6780 NEXT I 
6790 ©CLOSE 
6800 RETURN 


SAVE DATA TO TAPE 

7000 REM •*# SAVE DATA TO TAPE 
7005 CLS 

7010 PRINT"Data will be recorded in accordance with the format of* 

7020 PRINT"ANOVA “JLH 

7030 PRINT i PRINT"Insert Data Tape into ESF-80 and press SPACE BAR* 
7040 IF INKEY* O * " GOTO 7040 

7050 PRINT I INPUT*Enter File • under which you want data reoorded’JJF 
7055 ©CLOSE 
7060 ©OPEN JT 
7070 ©PRINT HS 

7080 ON WS GOTO 7100, 7400,7700 

7100 ©PRINT LC 

7110 ©PRINT QL 

7120 ©PRINT IC 

7130 IF QL*2 GOTO 7250 

7140 ©PRINT IE 

7150 FOR 1*1 TO IC 

7160 FOR J-1TOIE 

7170 ©PRINT C(I,J) 

7180 NEXT J 
7190 NEXT I 
7200 ©CLOSE 
7210 RETURN 
7250 FOR 1*1 TO IC 
7255 ©PRINT IXtil 
7260 FOR J»1 TO IXtl) 

7270 ©PRINT CdJ) 

7280 NEXT J 
7290 NEXT I 
7300 ©CLOSE 
7310 RETURN 
7400 ©PRINT LRiLC 
7410 ©PRINT IRdCiIE 
7420 FOR 1*1 TO IR 
7430 FOR J*1 TO IC 
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7440 FOR K-l TO IE 

7450 ©PRINT XdiJiK) 

7460 NEXT K 

7470 NEXT J 

7480 NEXT I 

7490 ©CLOSE 

7500 RETURN 

7700 ©PRINT LRiLCiLG 

7710 ©PRINT IX.IR,IC,IG 

7720 FOR I-t TO IR 

7730 FOR J«t TO IC 

7740 FOR K-l TO IG 

7750 QPRINT XU,J,K) 

7760 NEXT iC 

7770 NEXT J 

7780 NEXT I 

7790 QCLOSE 

7800 RETURN 




10000 REM •«« TWO-WAY CLASSIFICATION 
10005 IF II»2 THEN GOTO 10070 

10010 PRINT TAB(27irANOVA“tPRINT TAB<17)i“ (Two-Way Classification)* 

10020 PRINT ©192rEntsr a reference nans for the row factor (up to 12 characters)" (INPUT LR 

10030 PRINT ©192, LW 5 PRINT.LW 5 PRINT ©192,"Enter the number of *!LR1* lavels-JdNFUT IR 

10040 PRINT Q192, LW ( PRINT ©192,"Enter a reference name for the colunw factor (up to 12 char.'sl'iINFUT LC 

10050 PRINT ©192, LW ! PRINT LW 1 PRINT ©192,‘Enter the number of "iLCi* levels'WNPUT IC 

10053 PRINT ©192, LW 

10060 PRINT Q192,*Entar tha nunbar of replication* in each call (up to 10)*{(INFUT IE 

10070 CLS I IP IE O 1 PRINT "Tha modal is a C ELSE GOTO 10130 

10080 PRINT* (1) Filed Effects Modal 

10090 PRINT* <2> Random Effects Modal 

10100 PRINT* (3) Rows Fixed, Columns Random 

10110 PRINT* (4) Columns Fixed, Rows Random 

10120 KP-896 ( KS»4 ( GOSUB 1000 I JM-IO 

10130 CLS I IF IE-1 THEN JM-1 

10133 IF n-2 PRINT ©0,'Calculating* 

10140 S*0 ! GS-0 

10150 FOR 1-1 TO IR 1 RS(1>*0 

10160 IF II-l THEN PRINT'For *JLHSI 

10170 FOR J-l TO IC 

10180 IF n»l THEN PRINT* For *JLCtf 

10190 U-<I-1)*IC ♦ J J PSOJ>0 

10200 FOR K-l TO IE 

10210 IF II-l THEN PRINT* enter replication ♦'iKCINFUT XttJ.K) 

10220 Y-X(I,J,K) 

10230 S*SrY£2 5 OSOS+Y i RS(I]»RS(IHT 1 FStUJ-PSdJMY 

10240 NEXT K 

10250 NEXT J 

10253 IF II-l THEN CLS 

10260 NEXT I 

10270 CLS i FRINT*Calailating* 

10280 OS-OS[2/(1R*IC*IE) 

10290 FOR J-l TO IC 
10300 CS(J)-0 
10310 FOR I-l TO IR 
10320 FOR K-l TO IE 
10330 CSU>CS(JHXCI,J,1C) 

10340 NEXT K 
10330 NEXT I 
10360 NEXT J 
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10370 TRSS=0 ! TCSS=0 ! TISS=0 

10380 FOR 1*1 TO IR ! TRSS*TRSS+RS(I>t2 ! NEXT I ! TRSS=TRSS/(IC*IE»-OS ! RDF=IR-1 I MR-7RSS/RDF 
10390 FOP. J*1 TO IC t TCSS=TCSS»CS(J)C2 ! NEXT J i TCSS=TCSS/(IR«IE)-OS ! CDF*IC-1 ! MC*TCSS/CDF 
10400 IF IEM ADF*0 i TISS*0 ! GOTO 10420 

10410 FOR K*1 TO IC»IR ! TISS*TISS«-PS(1C>C2 ! NEXT K ! TISS-TISS/IE-OS-TRSS-TCSS i ADF*RDF«CDF 1 MA-TI 
SS/ADF 

10420 TTS*S-OS ! TDF*IR«IC*tE-l 

10430 ESS*TTS-TRSS-TCSS-TISS 1 EDF*TDF-ADF-CDF-RDF ! ME*ESS/EDF 
10440 IF IE* I MA*ME 

10450 ON JM GOTO 10460, 104S0, 10500, 10520 
10460 FR*MR/ME ! FC*MC/ME ! FA*MA/M£ 

10470 GOTO 10530 

10460 FR*MR/MA ! FC»MC/MA I FA-MA/ME 
10490 GOTO 10530 

10500 FR-MR/MA J FC*MC/ME ! FA*MA/ME 

10510 GOTO 10530 \ 

10520 FR*MR/ME 1 FC-MC/MA i FA-MA/ME 
10530 CLS 1 RETURN 


THESE - WAY CMWiaCA TO M 

15000 REM •** THREE HAY CLASSIFICATION *** 

15002 IF 11*2 THEN GOTO 15085 
15005 IX*0 

15010 PRINT TAB127F ANOV A* t PRINT TAB(16r(Three-Hay Classification)' 

15020 PRINT 3192,‘Enter a reference name for the row factor (up to 12 characters)‘lINPUT LH 
15030 PRINT 0192, LW SPRINT LW 

15040 PRINT Q192,‘Enter tha number of *JLRi* levels (up to lOriUNPUT IR 
15045 IF IR>IX THEN IX*IR 

15050 PRINT 0192J.H i PRINT 0192,‘Enter a reference name for the column factor (up to 12 char,‘s)*:iNPUT LC 
15060 PRINT Q192, LH *. PRINT LH I PRINT 0192,‘Enter the number of “,LC|* levels (up to lOPJHNPUT IC 
15065 IF IOIX THEN IX-IC 

15070 PRINT 0192, LH I PRINT LH 1 PRINT 0192,'Enter a reference name for the group factor (up to 12 char.'s) * 
(INPUT LG 

15080 PRINT 0192, LH ( PRINT LH l PRINT 0192,‘Enter the number of ‘JLG| # levels (up to lOPJlINFUT IQ 
15082 IF IGMX THEN IX-IG 

15085 CLS ! S*0 ! OS*0 l IF n«2 THEN FRINTCalculating* 

15090 FOR 1*1 TO IX 

15095 RSUI-0 1 CS(IM> 1 GSCIW) 

15100 FOR 7*1 TO IX 

15105 RC(I,J)*0 1 RGU,J>*0 1 caa JM) 

15110 NEXT 7 

15115 NEXT I 

15t30 FOR 1-1 TO IG 

15140 IF 11*1 THEN PRINT’For ‘JLG1I 

15150 FOR 7*1 TO IR 

15160 IF II»1 THEN PRINT* For ‘ILRJJ 

15180 FOR K-1TOIC 

15190 IF II»1 THEN PRINT* Enter observation for ‘(LCSOHNPUT X(7,K,I) 

15200 Y-X(7,1C,I) 

15220 S»S*Yt2 i OS-OS+Y 

15230 R5(J)*RS(7HY 1 CS(K)*CS(K)*Y l GS(I)*GS<IHY 

15240 RC(I,K)*RC(7,K>-»Y ( RG(7,I1-RG(7,I>+Y 5 CG(K.I)*CG(K,IHY 

15250 NEXT K 

15260 NEXT 7 

15265 IF II-I THEN CLS 

15270 NEXT I 

15280 CLS ! PRINT‘Calculating* 

15290 OS*OS12/(IR*IC*XG) 

15300 TRSS*0 1 TCSS=0 J TGSS-0 

15310 FOR 1*1 TO IR ! TRSS*TRSS+RS(I)t2 1 NEXT 15 TRSS*TRSS/(IC»IGI - OS 1 RDF-IR-1 ! MR-TRSS/RDF 
15320 FOR 7*1 TO IC i TC5S*TCSS*CS(J>C2 1 NEXT 7 1 TCSS*TCSS/(IR»IG> - OS ! CDF»IC-1 ! MC*TCSS/CDF 
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15330 FOR K*t TO 10 : TGSS*TGSS-*GS<K)C2 i NEXT K : TGSS-TGSS/<IR*IC) - OS 1 GDF-IG-l : MG-TGSS/GDF 
15340 DOSS=0 S DBSS*0 1 AISS-0 
15380 FOR 1*1 TO IX 
15390 FOR J*1 TO IX 

15400 DOSS*DOS5*RC(I,J>C2 ! DBSS=DBSS*RG(I,J>C2 ! AISS-AISSrCG(I,J)C2 
15410 NEXT J 
15420 NEXT I 

15430 DOSS*DOSS/IG - OS - TRSS - TCSS 1 DDF*(IR-l>*aC-l) ! MD-DOSS/DDF 
15440 DBSS=DBSS/IC - OS - TRSS - TGSS 1 SDF»UR-1>*UG-1) 1 MB*DBSS/BDF 
15450 AISS*AISS/IR - OS - TC5S - TGSS i ADF*<IC-1)«(IG-1> 1 MA*AISS/ADF 
15455 TSS*S-OS 5 TDF*IR«IC*IG-1 

15460 ESS=TSS-TRSS-TCSS-TGSS-DOSS-DBSS-AISS 1 EDF=TDF-RDF-CDF-GDF-DDF-BDF-ADF ! ME*ESS/EDF 

15470 Fl-MR/ME 

15480 F2»MC/ME 

15490 F3*MG/ME 

15500 F4*MD/ME 

15510 F5-MB/ME 

15520 F6*MA/M£ 

15530 RETURN 


PROGRAM SELECTION MENU 

20000 REM ••• PROGRAM SELECTION MENU 
20010 CLS 

20020 PRINT TAB<20>-PROGRAM SELECTION MENU- 

20030 PRINT 0128,*41 Data I/O, Transformations, Prob. & Inv. Prog. Distrlbutions'lPRIHT 
20040 PRINT*#2 Sample Slat's, 7-Num-Sum, Flots-IPRINT 

20050 PRINT"*3 Parametric Tests of Hypotheses and Confidence Intervals'lPRINT 

20060 PRINT-44 Robust Estimates and Confidents Intervals’lPRINT 

20070 PRINT-45 Multiple Linear Regression-IPRINT 

20080 PRINT’46 Analysis of Variance- 

20090 KP-896 1 KS=6 1 GOSUB 1000 

20110 IF 106 GOTO 5 

20120 IF 10*5 THEN LT--3* ELSE LT--1* 

20125 CLS 

20130 IF 1005 THEN PKINT'Frorn program 46 access to programs #1 thru 44 Is thorugh the Control Module in 
program 41.* ! PRINT 
20150 10*1 

20170 PRINT-INSERT TAPE 4 *1LT;* INTO ESF-90 AND PRESS SPACE BAR- 
20180 IF INKEY* O * * GOTO 20180 

20185 IF LTO-3* THEN PRINT-Program 41 not* being loaded...- 
20190 QLOAD 10 


END OF PROGRAM 46 
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